PharmaCalc

Pharmaceutical Calculation Suite

← All calculator guides Open calculator →
📅

Stability / Shelf Life (ICH Q1E)

Linear regression analysis to establish product expiration dating. Determines when 95% CI crosses acceptance limit.

ICH Q1E Long-Term Accelerated
🏭

Industry Use

Establishing Expiration Dates

Stability testing evaluates the degradation of an active pharmaceutical ingredient (API) and the overall product quality over time under defined storage conditions. ICH Q1E guidance provides the statistical methodology for analyzing stability data and establishing the shelf life (expiration dating period) that ensures the product meets specification throughout its labeled storage period.

Stability Testing Conditions (ICH Q1A/Q1B):

Regulatory Requirements:

ICH Q1E 95% CI Method 3–5 Year Shelf Life
🔢

Calculation Explanation

Linear Regression and Confidence Interval

ICH Q1E recommends linear regression of assay results (% label claim) versus time (months). The slope indicates degradation rate. The shelf life is the time at which the lower 95% confidence limit of the regression line intersects the acceptance limit.

Linear Regression Equation:

Linear Fit
ŷ = a + b·t

Where ŷ = predicted assay (% LC), a = intercept (assay at t=0), b = slope (%/month), t = time (months).

Lower 95% Confidence Interval:

95% CI Lower Bound
Lower CI = ŷ − t(α,n−2) × s × √[1/n + (t − t̄)²/Σ(tᵢ − t̄)²]

Where:

Shelf Life Determination:

Shelf Life Calculation
Set Lower CI = Acceptance Limit (e.g., 90%)
Solve for t (months)

The calculated shelf life is typically rounded down to a standard interval (12, 24, 36, 48, 60 months) for commercial labeling.

Linear Regression 95% CI Acceptance Limit
📋

Variable Definitions

Stability Regression Parameters

Symbol Parameter Name Units Description
tᵢ Time Points months Sampling time points for stability testing (e.g., 0, 3, 6, 9, 12, 18, 24 months).
yᵢ Assay Results % Label Claim API content at each time point, expressed as %LC. Initial (t=0) typically ≈100%; decreases over time due to degradation.
a Intercept % LC Predicted assay at time zero (t=0). Should be close to 100% for properly manufactured product.
b Slope %/month Rate of assay change per month. Negative value indicates degradation. Steeper (more negative) = faster degradation.
Coefficient of Determination 0–1 Measure of fit quality. R² > 0.99 indicates excellent linear fit; R² > 0.97 is typically acceptable.
Shelf Life Expiration Period months Time when lower 95% CI equals acceptance limit. Rounded to standard interval (12, 24, 36, 48, 60 months).
Acceptance Limit Lower Specification Limit % LC Minimum acceptable assay (e.g., 90% LC). Product must meet this throughout labeled shelf life.
👨‍🏫

Step-by-Step Tutorial

Worked Example with Long-Term Stability Data

Scenario: A pharmaceutical company is establishing the shelf life for a tablet product containing aspirin. Long-term stability data has been collected at 25°C/60% RH at 7 time points over 24 months. The acceptance limit is 90% LC. Linear regression and 95% confidence interval analysis are used to determine the expiration date.

Long-Term Stability Data (25°C/60% RH)
Month Assay (%LC)
0 100.2
3 99.8
6 99.1
9 98.6
12 98.0
18 97.1
24 96.3
  1. Calculate Regression Statistics:
    Using linear regression (least squares method):
    n = 7 time points
    t̄ = (0+3+6+9+12+18+24)/7 = 72/7 = 10.29 months
    ȳ = (100.2+99.8+99.1+98.6+98.0+97.1+96.3)/7 = 689.1/7 = 98.44%
  2. Calculate Slope (b) and Intercept (a):
    Using standard regression formulas (detailed algebra omitted for brevity):
    b ≈ −0.174 %/month (negative slope = degradation)
    a ≈ 100.4 % (intercept at t=0)
  3. Calculate Fit Quality (R²):
    R² ≈ 0.9985 (excellent fit; nearly perfect linearity)
  4. Calculate Standard Error and Residual SD:
    Residual sum of squares (RSS) ≈ 0.85
    s (residual SD) = √[RSS/(n−2)] = √[0.85/5] ≈ 0.41%
  5. Determine Shelf Life (Lower 95% CI = 90%):
    At t = shelf life, the predicted assay is:
    ŷ = 100.4 + (−0.174) × t

    The lower 95% CI at that time point must equal 90%:
    90% = ŷ − 2.571 × 0.41 × √[1/7 + (t − 10.29)²/Σ(tᵢ − t̄)²]

    Solving iteratively (or using regression software):
    t (shelf life) ≈ 59.8 months (raw calculation)
  6. Apply Commercial Rounding:
    Raw shelf life: 59.8 months
    Rounded down to standard interval: 48 months (4 years)
✓ SHELF LIFE = 48 MONTHS (4 YEARS)
Linear Fit: R² = 0.9985
Slope: −0.174 %/month | Intercept: 100.4%
Expiration Dating Period Established

Interpretation: The stability data shows a linear, predictable decline in assay at approximately 0.174% per month at 25°C/60% RH. With this degradation rate, the product's lower 95% confidence limit reaches the 90% acceptance limit at approximately 60 months. For commercial labeling and a conservative margin, a 48-month (4-year) shelf life is assigned. The product will meet the 90% LC specification throughout this period when stored at room temperature.

↪ Open Stability Calculator