Jump to content
  • Spotfire Statistica® Stability Analysis


    Stability analysis is the study of how drug product potency degrades over time. The primary statistical quantity of interest is the expiration date, or shelf life. Typically, a drug product is manufactured in batches. When estimating the shelf life of a medication, it is necessary to evaluate how the batches differ with respect to the potency degradation of the drug product over time.

    Stability analysis is the study of how drug product potency degrades over time. The primary statistical quantity of interest is the expiration date or shelf life. Typically, a drug product is manufactured in batches. When estimating the shelf life of a medication, it is necessary to evaluate how the batches differ with respect to the potency degradation of the drug product over time. Specifically, three different model types are considered.

    stability_node.thumb.png.41ddba9901c0fa4dc7b1e38690bfaf4b.png

    Model 1: Common Regression across All Batches

    This model is appropriate when the intercepts and slopes do not differ across batches and the data therefore can be pooled in order to estimate a common intercept and common slope.

    stability-1_0.png.08649c30cb56ebf7806003ace8043744.png

    Model 2: Separate Intercepts and Common Slope

    This model is appropriate when the intercepts differ across batches, but the slopes are common across batches.

    stability-2_0.png.29e9f6503b6be723b56e68be7ff8ec57.png

    Model 3: Separate Intercepts and Separate Slopes.

    This model is appropriate when the intercepts and slopes differ across batches. Additionally, the Mean Square Error (MSE) is not pooled across different batches, that is, homogeneity of variance across batches is not assumed.

    stability-3_0.png.10551967ae7496a06c31b5c457d5ee8a.png

    Parameter Estimation

    Parameters are estimated by ordinary least squares. This can be accomplished by solving the normal equations:

    stability-4_1.thumb.png.01346031d785f145d3927e7edaecf729.png

    t-Statistic

    The t-statistic associated with the null hypothesis that the ith parameter is equal to 0 is given by:

    stability-5_0.png.a922682df6710890274558a3b67ddc6f.png

    Confidence Interval for a Single Parameter

    A two-sided 95% confidence interval for the ith parameter is given by:

    stability-6_0.png.b26288504a2e33541cd1ac12e54be817.png

    Predicted Values

    The predicted value or estimated mean response for a given covariate vector Xi is given by:

    stability-7_0.png.f45a87ed7bd07c8d149d9b3b1dad3387.png

    Confidence Interval

    The two-sided 95% confidence interval for the mean response for a given covariate vector Xi

    is given by:

    stability-7_0.png.25adccfcf79a60e5fb66ddbe9bef59ad.png

    For a one-sided interval, then the 95th quantile of t distribution is used.

    Prediction Interval

    The two-sided 95% prediction interval for a new observation Xi is given by:

    stability-8.png.65326bcff669d62f0b2beae0306ff77d.png

    For a one-sided interval, then the 95th quantile of the t distribution is used.

    Model Selection

    Statistica selects the model by first evaluating if the slopes are equal across different batches.

    If the slopes' p-value is significant (by default, less than 0.25), the separate intercepts and separate slopes model is used to estimate shelf life.

    If the slopes' p-value is not significant, Statistica determines if the intercepts are equal across different batches. If the intercepts' p-value is significant (by default, less than 0.25), the separate intercepts, common slope model is used. If neither the slopes' p-value nor the intercepts' p-value is significant (by default, less than 0.25) the common intercept, common slope model is used.

    The p-values are computed by using the type 1 sums of squares decomposition of the full linear model with effects placed into the model in the following order: time, batch, batch by time interaction. The slopes' p-value is associated with the batch by time interaction and

    the intercepts' p-value is associated with the batch effect.

    Shelf Life Estimation

    Shelf life estimation is performed by determining where the 95% lower and/or upper confidence or prediction interval intersects with the user-defined specification limit(s). If either Model 2 or Model 3 has been selected, shelf life is estimated separately for each batch. Shelf life is then defined as the shelf life associated with the worst-case batch, that is, the smallest shelf life estimate. By default, Statistica will base shelf life computations on a two-sided 95% confidence interval.

    Residual Analysis

    A user can optionally perform a residual analysis to determine the selected model's adequacy.

    A residual is defined as the difference between the observed and predicted value, that is:

    stability-10.png.2f2e3ad4337e71578424f2c614b5155b.png

    Additionally, studentized residuals are computed. The formula for the studentized residuals is given by:

    stability-11.png.0fce5ceee58f7ffb995e3c63e7a862f8.png

    normal probability plot of the (studentized) residuals is constructed as follows. First, the (studentized) residuals are rank ordered. From these ranks, Statistica computes z values (i.e., standardized values of the normal distribution) based on the assumption that the data come from a normal distribution:

    stability-12.png.16a85503adf1f7b34efb5a84a59f4f74.png

    These z values are plotted on the y-axis in the plot. If the observed (studentized) residuals (plotted on the x-axis) are normally distributed, all values should fall onto a straight line. If the residuals are not normally distributed, they will deviate from the line. Outliers can also become evident in this plot. If there is a general lack of fit and the data seem to form a clear pattern (e.g., an S shape) around the line, the variable may have to be transformed in some way (e.g., a log transformation to "pull?in" the tail of the distribution, etc.)

     


    User Feedback

    Recommended Comments

    There are no comments to display.


×
×
  • Create New...