PharmaCalc

Pharmaceutical Calculation Suite

← All calculator guides Open calculator →

Stability / Shelf Life (ICH Q1E)

Fits assay against time, finds where the one-sided lower 95% confidence limit crosses the acceptance criterion, and applies the ICH Q1E cap on extrapolation beyond the observed data.

ICH Q1A(R2) ICH Q1E Linear Regression Shelf-Life Extrapolation
Get a shelf life from the ICH Q1E confidence limit, with optional α = 0.25 poolability analysis.
Free demo — no account needed. Every run produces a GMP PDF report.
Open the Stability & Shelf Life calculator →
Open the Stability & Shelf Life calculator →
🏭

Industry Use

How this calculation is applied in pharmaceutical manufacturing

This guide explains shelf-life estimation from long-term stability data using linear regression, following ICH Q1E, Evaluation for Stability Data. The method fits assay (% label claim) against time, computes the one-sided lower 95% confidence limit on the regression line, finds where that limit crosses the acceptance criterion, and then applies the Q1E limit on extrapolation beyond the observed data.

Two guidelines govern this work. ICH Q1A(R2) defines the study design — batches, storage conditions, time points, and the data package required at submission. ICH Q1E defines how to evaluate the data and how far a shelf life may extend beyond it.

Where this is used:

Storage conditions (drug product, general case, ICH Q1A(R2) §2.2.7.1):

StudyConditionMinimum data at submission
Long-term25°C ± 2°C / 60% RH ± 5% RH or 30°C ± 2°C / 65% RH ± 5% RH12 months
Intermediate30°C ± 2°C / 65% RH ± 5% RH6 months
Accelerated40°C ± 2°C / 75% RH ± 5% RH6 months

If 30°C/65% RH is chosen as the long-term condition, there is no intermediate condition. Q1A(R2) requires a minimum of 12 months of long-term data on at least three primary batches at submission — long-term data through the full proposed shelf life are not required; Q1E exists precisely because extrapolation beyond the observed data is permitted within defined limits.

Calculation Explanation

Mathematical basis and regulatory foundation
Regression model (ordinary least squares)
ŷ = a + b·t   ·   s = √( RSS / (n − 2) )

Assay y (% label claim) versus time t (months); s is the residual standard deviation.

One-sided lower 95% confidence limit
Lower bound(t) = a + b·t − t(0.95, n−2) · s · √[ 1/n + (t − t̄)² / Σ(tᵢ − t̄)² ]

For an attribute that decreases with time, Q1E §2.6 directs that “the lower one-sided 95 percent confidence limit should be compared to the acceptance criterion.” t(0.95, n−2) is the one-sided 95th-percentile Student's t (for df = 5, 2.015). Using the two-tailed α = 0.05 value (2.571 at df = 5) is a common error; it produces a 97.5% one-sided limit.

Shelf-life estimate and extrapolation cap
Solve Lower bound(t) = criterion  →  then cap at min( 2 × Tdata, Tdata + 12 )

Solve numerically for the crossing (e.g. at 90.0% LC), then apply the ICH Q1E §2.4.1.1 limit: for long-term data amenable to statistical analysis, the proposed shelf life “can be up to twice, but should not be more than 12 months beyond, the period covered by long-term data.”

Multiple batches: shelf life is established on a minimum of three batches. Before pooling into one regression, test poolability (equality of slopes and intercepts, e.g. ANCOVA); Q1E directs a significance level of 0.25 to compensate for the low power of the design. If batches cannot be pooled, the shelf life is based on the worst-performing batch. In the calculator the poolability analysis is opt-in: enable it to run the α = 0.25 tests; left off, multiple batches are evaluated individually and the worst-performing batch governs — conservative and Q1E-acceptable, but no poolability test is performed.

ICH Q1A(R2) ICH Q1E §2.6 ICH Q1E §2.3 / §2.4.1.1
📋

Variable Definitions

All input parameters and their meaning
Symbol Variable Units Description
y, t Assay result; storage time % LC; months Long-term stability results for the attribute being trended, at each scheduled pull.
a, b Intercept; slope % LC; %/month Ordinary least-squares estimates of the change over time.
s Residual standard deviation % LC √(RSS/(n−2)); the scatter of the data about the fitted line.
t(0.95, n−2) One-sided 95th-percentile Student's t dimensionless The one-sided value — 2.015 at df = 5 — not the two-tailed 2.571.
criterion Acceptance criterion % LC The specification limit the confidence bound is compared against (e.g. 90.0% LC for assay).
Tdata Long-term data coverage months The period covered by long-term data; the Q1E cap is min(2 × Tdata, Tdata + 12).
📖

Step-by-Step Tutorial

Worked example — single batch, arithmetic illustration only

Scenario: long-term data at 25°C/60% RH, acceptance criterion ≥ 90.0% LC. A real submission requires at least three batches and a poolability assessment; a single batch is shown here to illustrate the arithmetic only.

  1. Assemble the long-term data.
    Example Data
    t (months): 0, 3, 6, 9, 12, 18, 24
    y (% LC): 100.2, 99.8, 99.1, 98.6, 98.0, 97.1, 96.3
  2. Fit the regression.
    Calculation
    n = 7; t̄ = 10.2857; ȳ = 98.4429; Sₓₓ = 429.4286; Sₓᵧ = −71.4857
    Slope b = −71.4857 / 429.4286 = −0.16647 %/month
    Intercept a = 98.4429 + 0.16647 × 10.2857 = 100.155 % LC
    RSS = 0.0771 → s = √(0.0771/5) = 0.1242 · R² = 0.9936
  3. Take the one-sided t value: t(0.95, 5) = 2.015 — one-sided, per Q1E §2.6.
  4. Find where the lower bound crosses 90.0% LC.
    Calculation
    Solving Lower bound(t) = 90.0 numerically gives t ≈ 57.5 months.
    (The point estimate alone crosses at 61.0 months — the confidence bound must cross earlier.)
  5. Apply the extrapolation cap — do not skip this step.
    Calculation
    Data cover 24 months → ceiling = min(2 × 24, 24 + 12) = 36 months
    The statistical crossing point is not the shelf life. Q1E §2.4.1.1: the proposed shelf life “can be up to twice, but should not be more than 12 months beyond, the period covered by long-term data.” Q1E §2.3 adds that extrapolation assumes the same change pattern continues, and that an extrapolated shelf life must be verified by additional long-term data as it becomes available.
  6. Conclude.
    Interpretation
    Supportable shelf life: 36 months — not the 57.5-month statistical crossing.
✓ Supportable shelf life = 36 months (confidence-limit crossing 57.5 months, capped by ICH Q1E at 24 + 12)
⚠️

Interpretation, Assumptions and Limitations

What the calculation presumes — and what it does not cover
📚

References

Primary sources for this guide
Stop rebuilding this in a spreadsheet

PharmaCalc returns a shelf life from the ICH Q1E confidence limit, with optional α = 0.25 poolability analysis. It is computed server-side against the published method and written into a GMP PDF report with the inputs, formula chain, references, document control and signature pages — traceable to the software release that produced it.

Open the Stability & Shelf Life calculator — free demo
No account required to try it · Browse all 23 calculator guides
Open the Stability & Shelf Life calculator →