Comparison of shelf life estimates generated by asap prime tm with the king kung fung approach
1 / 16

Juan Chen 1 , Sabine Thielges 2 , William R. Porter 3 , Jyh-Ming Shoung 1 , Stan Altan 1 - PowerPoint PPT Presentation

  • Uploaded on

Comparison of Shelf Life Estimates Generated by ASAP prime TM with the King-Kung-Fung Approach. Juan Chen 1 , Sabine Thielges 2 , William R. Porter 3 , Jyh-Ming Shoung 1 , Stan Altan 1. 1 Nonclinical Statistics and Computing, Janssen R&D 2 BE Analytical Sciences and COES, Janssen R&D

I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.
Download Presentation

PowerPoint Slideshow about ' Juan Chen 1 , Sabine Thielges 2 , William R. Porter 3 , Jyh-Ming Shoung 1 , Stan Altan 1' - kolya

An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript
Comparison of shelf life estimates generated by asap prime tm with the king kung fung approach

Comparison of Shelf Life Estimates Generated by ASAPprimeTM with the King-Kung-Fung Approach

Juan Chen1, Sabine Thielges2, William R. Porter3, Jyh-Ming Shoung1, Stan Altan1

1Nonclinical Statistics and Computing, Janssen R&D

2BE Analytical Sciences and COES, Janssen R&D

3Peak Process Performance Partners LLC


  • Arrhenius Equation and Extensions to include Humidity Effects

    • Extended Arrhenius Equation

    • Extended King-Kung-Fung (KKF) Model

  • ASAP Approach

    • Zero order, 2-temperature Example

  • Case Study using Pseudo-data

    • DoE of Pseudo-data

    • Comparison of Outputs from ASAP and KKF Model

  • Conclusion

Arrhenius equation
Arrhenius Equation

Named for Svante Arrhenius (1903 Nobel Laureate in Chemistry) who established a relationship between temperature and the rates of chemical reaction:

Where kT= Degradation Rate

A = Non-thermal Constant

Ea= Activation Energy

R = Universal Gas Constant (1.987 cal/mol)

T = Absolute Temperature

Arrhenius equation with humidity term
Arrhenius Equation with Humidity Term

A humidity term with coefficient B is introduced to account for the effect of relative humidity on rate parameter.

activation energy

humidity sensitivity factor

degradation rate

Pre-exponential factor

gas constant (1.987cal/mol)

Extended king kung fung model
Extended King-Kung-Fung Model

King-Kung-Fung (KKF) model is widely used for analyzing accelerated stability data

Let T =298oK (25oC)

H = 60

  • Directly estimate Shelf Life (SL) at 25C/60%RH and its uncertainty

  • Parameter estimates are calculated based on the Arrhenius relationship conditional on an assumed zero order kinetic

  • Can be further extended to a nonlinear mixed model context and Bayesian calculation

Introduction to asap prime tm
Introduction to ASAPprimeTM

The ASAPprimeTMcomputerized system is a computer program that analyzes data from accelerated stability studies using a 2-week protocol accommodating both temperature and humidity effects through a extended Arrhenius model.

The program makes a number of claims:

  • Reliable estimates for temperature and relative humidity effects on degradation rates,

  • Accurate and precise shelf-life estimation,

  • Enable rational control strategies to assure product stability.

    These claims require careful statistical considerations of the modeling strategies proposed by the developers.

    Our objective is to evaluate the first two claims in relation to widely accepted statistical approaches and considerations.

Illustration of ASAPprimeTMApproach Shelf Life (SL) Estimation using Zero-order, 2-Temperature Example

Data SD Hierarchy:

1. Calculate from replicate data , if >LOD

2. User-defined SD (fixed) or RSD, if >LOD

3. Default 10%RSD, if >LOD

4. LOD

  • Calculate SD of SL at accelerated conditions = Mean SL – Extrema SL

  • Method of SD calculation is not consistent with the standard definition of a SD

  • No a priorivariance structure proposed for analytical variability in terms of a statistical model  The uncertainty in the SL estimates cannot be understood in relation to statistical principles; empirical comparisons only

Mean SL

Extrema SL

Error propagation through a MC simulation

Pairs of lnK at 50C and 60C form 49 regression lines across 1/T with slopes (=Ea/R) and intercepts (lnA)

  • Calculate lnK at 25C for each simulation: lnK = lnA – Ea/(RT)

  • Simulation was drawn from an undocumented distribution.

  • Arrhenius parameters (lnA and Ea) are determined from simulated degradation rates Statistical properties of the estimated lnA and Eaare not unknown.

ASAP Probability statement about lnK and SL at 25C

  • SL at 50C and 60C previously simulated from an unknown distribution Cannot verify lognormal distribution at 25C

  • Model fitting cannot be confirmed by standard statistical procedures.

Case Study using Pseudo-data

  • 4 x 4 factorial design of temperature and humidity

  • 9 sets of data simulated from combinations of D0, Ea, B values each at L, M, H, 73-day sampling design assuming a zero-order model

  • lnA back-calculated to obtain SL at 2 years at 25C/60%RH

  • Normally distributed random errors with mean 0 and SD 0.1 were added as analytical variability

  • KKF model carried out in SAS Procnlin(Initial values: SL = 0.27 year, Ea = 10, B = 0.001, D0 = 0)

  • KKF estimated residual errors range from 0.08 to 0.12, whereas the true value is 0.1

  • ASAP different user specified SD affect probability statement about SL

ASAP generally underestimated lnA.

ASAP generally underestimated Ea.

Standard errors of Ea are affected by user specified SD, and are generally larger than KKF estimates.

ASAPunderestimated B for data HHH and HLH, and overestimated B for LLL.

Standard errors of B are affected by user specified SD, and are generally larger than KKF estimates.


  • Uncertainty measure for ASAP estimated shelf life is derived from an “error propagation” calculation using either replicate error or a user defined quantity.

    • Statistical rationale for uncertainty limits is not clear.

  • Does not lead to a statistical confidence statement.

  • ASAP simulation of SL to predict room temperature SL:

    • Underlying distribution of SL at accelerated conditions is not documented.

    • The precision of Arrhenius model parameter estimates is influenced by user specified SD and cannot be validated statistically.

    • Overall model fitting is unclear and lacks documentation.

    • Manufacturing variability cannot be accommodated.

  • ad