From the Garrett Two-Stage Model to the Bayesian View: A Historical Survey of Stability Modeling 1954 to 2017
|
|
- Jason Howard
- 6 years ago
- Views:
Transcription
1 From the Garrett Two-Stage Model to the Bayesian View: A Historical Survey of Stability Modeling 1954 to 2017 William R. Porter, PhD Principal Scientist, PPPP-LLC@comcast.net Stan Altan, PhD Manufacturing, Toxicology and Applied Statistics, Janssen R&D SALTAN@its.jnj.com
2 Outline Early period : Accelerated stability assessment using Arrhenius models to estimate shelf life at room temperature. First Transitional period : Moving from accelerated studies to RT studies Real time room temperature period : 1987 FDA guidance and 2004 Q1E. Second transitional period 2010 current: Movement from real time studies back to accelerated studies with focus on degradants. Example real time and accelerated studies 2
3 Arrhenius, 1889 Reaction rate depends on temperature! Z physik Chemie 4: (1889) That is, k T k T T T 0 exp C kt 0 T0T1 1 exp C T T 1 3
4 Prehistory: Bigelow, 1921 Early attempt to relate time-tofail to storage temperature J Infec Dis 29: (1921) 4
5 Prehistory: Higuchi et al., 1950 Graphical analysis of data no mention of any statistical modeling. Note that t½ was used instead of rate k. 5
6 McBride & Villars, 1954 Anal Chem 26(5): (1954) First 2-stage regression model for Arrhenius kinetics Rates were calculated in Stage One Logarithms of rates were calculated for Arrhenius plot Weighted regression used for Stage Two 6
7 Garrett, 1954 (1) J Amer Pharm Assoc 43(9): (1954) The Ur paper for statistical analysis of accelerated stability assessment studies. McBride & Villars were ignored for years. 7
8 Garrett, 1954 (2) Details of the statistical modeling did not appear in Garrett s first paper. 8
9 Huyberechts et al., 1955 Don t publish in French if you want to be cited. Uses weighted 2 nd stage (McBride & Villars) 9
10 Garrett, 1956A J Amer Pharm Assoc 45(3): (1956) First paper by Garrett with details of how to fit Arrhenius model to data using two-stage unweighted regression. Tedious explanation. 10
11 Garrett, 1956B J Amer Pharm Assoc 47(7): (1956) The predictions work 1 st real time stability data plot with confidence band. 11
12 McLeod et al., 1958 H. A. McLeod, O. Pelletier and J. A. Campbell., "The Prediction of Expiration Dates for Multivitamin Preparations by Accelerated Storage Tests." Canadian Pharmaceutical Journal Scientific Section 91: (1958). The authors provide complete worked out numerical details with raw data and all statistical calculations using a 2-stage Garrett pseudo-zero-order kinetic model. Individual pseudo-zero-order rates and initial values (time zero potency) were obtained for each elevated storage temperature. The rates were converted to base 10 logarithms. The base 10 logarithms were fitted to an Arrhenius model of the form log 10 (rate) = intercept - slope reciprocal absolute temperature. NO WEIGHTING was used, and the authors cite Edward Garrett's paper. The fitted Arrhenius model was used to generate an estimate of the base 10 logarithm of the rate at some other temperature and the standard error of the base 10 logarithm of the rate was also calculated at that temperature. The rate at the new temperature was obtained by taking antilog of the base 10 logarithm of the projected rate at the new temperature. The standard error of the rate at the new temperature was calculated by taking the antilogs of the difference between the base 10 logarithm of the rate plus its standard error (in log units) and the base 10 logarithm of the rate minus its standard error (in log units) and then dividing the result by 2. The shelf life was calculated by taking the initial potency and the potency at the end-of-shelf-life as fixed constants, finding their difference, and dividing this difference by the estimate rate plus or minus its estimated standard error. The standard error of the shelf life was calculated by calculating the shelf life at rate plus standard error and at rate minus standard error, then finding the difference and dividing by two. Student's t distribution was used to generate a confidence interval estimate using T 2 d.f., where T is the number of different rates used in the second stage Arrhenius model. 12
13 Box & Lucas, 1959 Biometrika 46:77 90 (1959) Optimal design for Arrhenius model: Impractical Requires initial observation at 1/T 13
14 Toothill, 1961 Uses common intercept, isoconversion design to estimate rates in first stage. Adaptation of bioassay methods (Finney). Point & interval estimates for shelf life. 14
15 FDA cgmps, August,1969 Federal Register, 34(161): , August 22, 1969 No mention of statistical methods for evaluating shelf life assignment 15
16 Conference on Dating of Pharmaceuticals, October,1969 Comer JA. Processing of stability data for FDA submission. The Dating of Pharmaceuticals. Univ. Wisconsin pp , 1969 The computer system in use when this paper was written was the IBM The paper tape is read on the paper tape reader into the processing unit and onto magnetic tape. A cooperative effort of academic, governmental, and industrial organizations is needed to improve these systems until they become a common means of assisting the pharmaceutical scientist in the planning, interpretation, and submission of stability data. NOTE CONFIDENCE LIMITS USED TO DEFINE SHELF LIFE! 16
17 Carstensen & Su, 1971 Bull Parent Drug Assoc 25(6): ()1971) 1 st nonlinear regression fitting of Arrhenius model. Use pseudo-zero-order model if degradation is less than 10%. Use pseudo-first-order model if degradation exceeds 10%. 17
18 Carstensen & Nelson, Carstensen JT, Nelson E. Terminology regarding labeled and contained amounts in dosage forms. J Pharm Sci 65(2): (1976) Incorporated into 1987 FDA Stability Guidance 18
19 Davies & Budgett, 1980 (1) Distribution of Arrhenius parameter estimates: Prior distributions, posterior distributions. Lognormal distribution for Arrhenius frequency factor. Lognormal distribution for shelf life. 19
20 Davies & Budgett, 1980 (2) Prior distributions! 20
21 Davies & Budgett, 1980 (3) Prior distributions (cont.) Monte Carlo confirmation of posterior distributions 21
22 Davies & Budgett, 1980 (4) Posterior distributions 22
23 Statistics In the Pharmaceutical Industry, 1981 Editors: Charles Ralph Buncher, Jia-Yeong Tsay First edition covers ONLY accelerated stability studies. NO MENTION of real time studies. 23
24 King, Kung & Fung, 1984 J Pharm Sci 73(5): (1984) Direct point and interval estimation of shelf life at normal storage temperature from accelerated stability assessment. 24
25 FDA Stability Guideline, 1987 First guidance with detailed description of appropriate statistical methods for analysis and interpretation of real time stability data. Shelf life (expiration dating period) set using Carstensen & Nelson (1976). Fit pseudo-zero-order kinetic model to loss of potency as a function of time using linear regression. Construct one-sided lower confidence bound for the mean potency as a function of time. Shelf life is time when LCB reaches specification limit. Other approaches can be used but must be justified. 25
26 Statistics In the Pharmaceutical Industry, 2 nd Edition, 1991 Editors: Charles Ralph Buncher, Jia-Yeong Tsay Second edition covers BOTH accelerated stability studies AND real time studies. Only Frequentist methods described. 26
27 Shao & Chow, 1991 Release limits based on stability considerations; focus on potency loss. 27
28 Chow & Shao 1991 Confidence limits evaluated using constant normal prior by Monte Carlo simulation. 28
29 Shao & Chow, 1994 Confidence limits evaluated using constant normal prior by Monte Carlo simulation. 29
30 Su et al., 1994 Studied effect of varying kurtosis on estimates of ln(k) at various T conditions. ln(a) E a DID NOT look at prior distribution of raw data! 30
31 Shao & Chow, 2001 Three methods for estimating confidence limit approach to shelf life estimation using constant normal prior distribution evaluated using Monte Carlo simulation. 31
32 Statistics In the Pharmaceutical Industry, 3rd Edition, 2006 Editors: Charles Ralph Buncher, Jia-Yeong Tsay Third edition covers ONLY real time studies. No mention of accelerated studies. Only Frequentist methods described. 32
33 Chen, Zhong & Nie, 2008 Used uniform prior distributions for initial value and pseudo-zero-order degradation rates. 2 nd stage of 2-stage model. Did not consider estimation from raw data (1 st stage of 2-stage model). 33
34 Statistics in the Pharmaceutical Industry Edited by Charles Ralph Buncher, Jia-Yeong Tsay Second Edition 1991 RT and Accelerated First Edition 1981 Accelerated Studies Third Edition 2006 RT Studies 34
35 EXAMPLE Real Time Studies: Mixed Model Framework Mechanistic and empirical basis exists to forego pooling tests (constrained models by strength is recommended). Assume a fixed common package-temperaturespecific slope. Assume different batch-specific Intercepts. Main requirement is to estimate the parameters and account for incipient variation in such a way that control over the lot mean is assured. 35
36 Mixed Effects Model Model estimates: Process Mean at time of Manufacture Rate parameter Variance Structure Process (Batch-Batch) Analytical Variation Number of Parameters Number Form Fixed Effects Variance 1 n yijk ( 0 i ) Bj Tijk c +1 2 ijk n c +1 3/4 a 2 yijk ( 0 i ) B j i Tijk ijk Index: i=batch, j=condition, k=time. Model 4: Random Term in the Intercept, Model 5: Random Terms in Intercept and Slopes a if correlated Random Terms in Intercept and Slope Main objection small number of batches. 36
37 Case Study (Assay) 12 months stability study for the Assay of 3 Drug Product Batches at 2 storage conditions. Estimate Degradation Rate and Shelf Life for Specification limits, 90% to 110%. 37
38 Bayesian Approach Bayesian approach provides a flexible framework for incorporating scientific and expert judgment, exploiting past experience with similar products and processes and a more natural way to approach decision making. Posterior Distribution: X ~ f ( ), ~ prior( ) post( x) f f ( x f ( x ( x ) prior( ) ) prior( ) ) prior( ) d Posterior Distribution ~ Likelihood x Prior Distribution 38
39 Expiration Date Mixed Model 1 with random batch term on intercept, condition-specific slope: Y ijk ( ) X The expiration date is the solution to the equation: LSL ˆ ˆ T ED c Bayesian approach, consider: LSL T Z or T i SL j Z is independent of data and symmetric about 0: 2 Z ~ N(0, ). ijk ijk Var ˆ ( ˆ ˆ T SL ED ) ˆ 2 LSL Z 39
40 Prior Distributions and Simulation Parameters Prior Distributions: Expert opinion Process mean is likely between 99% and 101%: ~ N(100, 0.1) Lot to lot variance is likely between 0.1 and 0.5: 2 1 ~ (10, 2) Flat prior on the yearly degradation rates: 1, 2 ~ I(, ) Analytical variance is likely between 0.1 to 1.0: R/WinBUGS Simulation Parameters: 2 1 ~ (6, 3 chains, 500,000 iterations/chain, discard 1 st 100,000 simulated values in each chain. Retain every 100 th simulation draw, 27,000 simulated values for each parameter. 2) 40
41 Comparison of Parameter Estimates & Expiration Dates Frequentist Bayesian Parameters 95% Confidence Mean 95% Credible Estimate Interval (Median) Interval , , , , , , (0.20) 0.11, (0.26) 0.17, 0.42 Shelf Life Expiration Date Storage Condition 25C 30C 25C 30C Frequentist Bayesian 94 ( 89 ) 64 ( 62 ) Bayesian estimate of shelf life is based on the posterior mean (median) and the expiration date corresponds to the 5%-tile value. Specification limits, 90% to 110%. 41
42 Shelf Life Distribution Bayesian Approach 42
43 EXAMPLE Accelerated Stability Testing Product is subjected to stress conditions. Temperature and relative humidity are the most common stress factors. Purpose is to predict long term stability and shelf life. Arrhenius equation captures the kinetic relationship between rates and temperature. The usual fixed and mixed models used in real time studies ignore any relationship between rate and temperature. 43
44 Arrhenius Equation with Relative Humidity Term A humidity term with coefficient B is introduced to account for the effect of relative humidity on rate parameter. degradation rate ln(k T,H ) = ln A Pre-exponential factor activation energy humidity sensitivity factor E a R T + B h gas constant (8.314 x 10-3 kj mol -1 K -1 ) 44
45 Nonlinear Parameterization (Extended King-Kung-Fung Model) k T,H = Ae E a R T +B H k T, H A k Let T =298 o K (25 o C) H = ,60 e Ea 1 1 B( H 60) R 298 T A = k 298,60 e E a 298 R B 60 Assuming zero order kinetics, total degradation is: D t = D 0 + k T,H t Q D 0 k298,60 t SL D t = D 0 + Q D 0 t SL t e Ea R T +B H 60 +ε Estimate Shelf Life at 25C/60%RH and its uncertainty w.r.t spec = Q. Parameter estimates are calculated based on the Arrhenius relationship conditional on an assumed zero order kinetic model. 45
46 Linearized Expanded Arrhenius Model Two-stage approach (Garrett/McBride & Villars, 1954), modified to include relative humidity. Assume a zero order kinetic model. Stage 1 : fit a pseudo zero order kinetic model to the concentration measurements: D T t D 0 Stage 2 : Model the rate estimates according to Arrhenius relationship: ln kt, h ln A t Expressed as linear regression problem: ln, E R k a T T 1 T k T h 2 B h h 46
47 Case Study of Degradant under Accelerated Conditions Data Table day Temp RH Deg
48 Bayesian Considerations for Model Specification Priors: D 0 ~ Uniform(0, 4) lne A ~ Uniform(2, 7) B ~ Uniform(0, 0.20) lnt SL ~ Uniform(1, 8) σ ε 2 ~ Inverse-Gamma(1.41, 0.055) JAGS specifications: Burn-in of 100,000 iterations followed by 3,000,0000 iterations. Thinning rate of chains. Total of 30,000 simulated values retained for each parameter. Randomly generated initial values for each parameter from corresponding prior distribution. 48
49 JAGS model specification 49
50 Accelerated Degradant Study Results Parameter (se) King, Kung, Fung Garrett Two- Stage Bayesian Posterior Estimates Mean Median 5% D initial 0.27 (0.40) 0.36(0.50) E a (kj Mol -1 ) (3.9) (4.4) B (RH) 0.06(0.002) 0.06(0.002) t SL, 25C/60%RH (Exp 107) (Day) (Exp 72) (Exp 56) t SL, 30C/75%RH (Exp 23) (Day) (Exp 17) (Exp 14) MSE
51 Comparison of Approaches Linearized Arrhenius Model (Garrett): Simple and does not require specialized software. Requires GLS approach for appropriate weighting. Nonlinear Model (King-Kung-Fung): Computationally intensive. Computing convergence issues. 51
52 Summary Current regulatory guidelines for assessing stability and shelf life claims are being challenged. Poolability paradigm is a concern. A Mixed Effects model is a natural representation of a manufacturing process. Bayesian framework can incorporate process engineering and scientific judgment. Addresses poolability concerns. Accelerated study models are being reassessed and updated through linear and non-linear models. King, Kung and Fung model has been extended to include a humidity term. Technology is evolving rapidly which will require greater statistical collaborations. 52
53 Future 1. Degradation rates are parameters, not data. 2. Consider effect of informative prior distribution on: a. Measurement and within-batch sampling errors: i. Repeatability (within test day errors). ii. iii. Intermediate precision (over testing interval errors). Reproducibility (between lab errors). b. Batch variation (manufacturing errors). c. Environmental errors (temperature, relative humidity). 3. Most drug products fail due to excess degradation products, not loss of potency. 4. A drug product is not stable unless some other analyst working in some other lab testing some other sample from some other batch stored in some other environment says so. 53
54 References (1) Regulatory Guidance 1. International Conference on Harmonization ICH Harmonised Tripartite Guideline: ICH Q1A(R2) Stability Testing of New Drug Substances and Products. Revision 2 2. International Conference on Harmonization ICH Harmonised Tripartite Guideline: ICH Q1D Bracketing and Matrixing Designs for Stability Testing of New Drug Substances and Products. 3. International Conference on Harmonization ICH Harmonised Tripartite Guideline: ICH Q1E Evaluation of Stability Data if New Drug Substances and Products. Revision 2 4. International Conference on Harmonization ICH Harmonised Tripartite Guideline: ICH Q6B Specifications: Test Procedures and acceptance Criteria for Biotechnological /Biological Products. 5. International Conference on Harmonization ICH Harmonised Tripartite Guideline: ICH Q8(R2) Pharmaceutical Development. 6. International Conference on Harmonization ICH Harmonised Tripartite Guideline: ICH Q9 Quality Risk Management. Approaches to Pooling 1. Chen WJ, Tsong Y. Significance levels for stability pooling test: a simulation study. Journal of Biopharmaceutical Statistics 2003; 13: Chow SC, Shao J. Test for batch-to-batch variation in stability analysis. Statistics in Medicine 1989; 8: Liu W, Jamshidian M, Zhang Y, Bretz F, Han XL. Pooling batches in drug stability study by using constant-width simultaneous confidence bands. Statistics in Medicine 2007; 26: Ruberg SJ, Stegeman JW. Pooling data for stability studies: testing the equality of batch degradation slopes. Biometrics 1991; 47: Ruberg SJ, Hsu JC. Multiple comparison procedures for pooling batches in stability studies. Technometrics 1992; 34: Tsong Y, Chen WJ, Lin TY, Chen CW. Shelf life determination based on equivalence assessment. Journal of Biopharmaceutical Statistics 2003; 13:
55 References (2) Mixed and random Models 1. Chow SC, Shao J. Estimating drug shelf-life with random batches. Biometrics 1991; 47: Chen JJ, Hwang JS, Tsong Y. Estimation of the shelf-life of drug with mixed effects models. Journal of Biopharmaceutical Statistics 1995; 5: Shao J, Chow SC. Statistical inference in stability analysis. Biometrics 1994; 50: Michelle Quinlan, W. Stroup, J.D. Christopher, J. Schwenke. 2013a. On the Distribution of Batch Shelf Lives. Journal of Biopharmaceutical Statistics, 23:4, Michelle Quinlan, W. Stroup, J. Schwenke, J.D. Christopher. 2013b. Evaluating the Performance of the ICH Guidelines for Shelf Life Estimation, Journal of Biopharmaceutical Statistics, 23:4, Robert Capen, J.D. Christopher, P. Forenzo, C. Ireland, O. Liu, S. Lyapustina, J. O Neill, N. Patterson, M. Quinlan, D. Sandell, J. Schwenke, W. Stroup, T. Tougas On the Shelf Life of Pharmaceutical Products. AAPS PharmSciTech, 13: Norwood, T.E.: Statistical analysis of pharmaceutical stability data, Drug Development and Industrial Pharmacy, 12(4), (1986) 8. Kiermeier, A.; Jarett, R.G., Verbyla, A.P.: A new approach to estimating shelf-life, Pharmaceutical Statistics, 3, 3-11 (2004) Bayesian Modeling 1. Chen, J, Zhong,J, Nie,L. Bayesian hierarchical modeling of drug stability data. Statist. Med. 2008; 27:
56 Acknowledgements Thanks to (Janssen): Madelyn Drevets Areti Manola Jyh-Ming Shoung 56
57 Questions? 57
Issues in Non-Clinical Statistics
Issues in Non-Clinical Statistics Stan Altan Chemistry, Manufacturing & Control Statistical Applications Team Department of Non-Clinical Statistics 1 Outline Introduction Regulatory Considerations Impacting
More informationAN ALTERNATIVE APPROACH TO EVALUATION OF POOLABILITY FOR STABILITY STUDIES
Journal of Biopharmaceutical Statistics, 16: 1 14, 2006 Copyright Taylor & Francis, LLC ISSN: 1054-3406 print/1520-5711 online DOI: 10.1080/10543400500406421 AN ALTERNATIVE APPROACH TO EVALUATION OF POOLABILITY
More informationEFFECT OF THE UNCERTAINTY OF THE STABILITY DATA ON THE SHELF LIFE ESTIMATION OF PHARMACEUTICAL PRODUCTS
PERIODICA POLYTECHNICA SER. CHEM. ENG. VOL. 48, NO. 1, PP. 41 52 (2004) EFFECT OF THE UNCERTAINTY OF THE STABILITY DATA ON THE SHELF LIFE ESTIMATION OF PHARMACEUTICAL PRODUCTS Kinga KOMKA and Sándor KEMÉNY
More informationBest Linear Unbiased Prediction: an Illustration Based on, but Not Limited to, Shelf Life Estimation
Libraries Conference on Applied Statistics in Agriculture 015-7th Annual Conference Proceedings Best Linear Unbiased Prediction: an Illustration Based on, but Not Limited to, Shelf Life Estimation Maryna
More informationImportance Of Accelerated Stability Study
Importance Of Accelerated Stability Study Accelerated stability testing All medicinal products decompose with time. Paradoxically, when this decomposition is being assessed the skilled formulator becomes
More informationAccelerated Stability Assessment Program (ASAP), a QbD case study
Accelerated Stability Assessment Program (ASAP), a QbD case study Applicazione del QbD nella produzione dei medicinali Università degli Studi di Milano, 6 Maggio 2016 1 Approaches for Shelf Life determination
More informationASEAN GUIDELINES FOR VALIDATION OF ANALYTICAL PROCEDURES
ASEAN GUIDELINES FOR VALIDATION OF ANALYTICAL PROCEDURES Adopted from ICH Guidelines ICH Q2A: Validation of Analytical Methods: Definitions and Terminology, 27 October 1994. ICH Q2B: Validation of Analytical
More informationASAP implementation and use Sabine Thielges
ASAP implementation and use Sabine Thielges Science of Stability Conference, 12-14 October 2015 Introduction to ASAP principles Accelerated Stability Assessment Program Accelerated stress study (~ 1 month)
More informationUniversity of California, Berkeley
University of California, Berkeley U.C. Berkeley Division of Biostatistics Working Paper Series Year 2003 Paper 127 Rank Regression in Stability Analysis Ying Qing Chen Annpey Pong Biao Xing Division of
More informationProcess Simulation, Parameter Uncertainty, and Risk
Process Simulation, Parameter Uncertainty, and Risk in QbD and ICH Q8 Design Space John J. Peterson GlaxoSmithKline Pharmaceuticals john.peterson@gsk.com Poor quality Better quality LSL USL Target LSL
More informationQuality by Design and Analytical Methods
Quality by Design and Analytical Methods Isranalytica 2012 Tel Aviv, Israel 25 January 2012 Christine M. V. Moore, Ph.D. Acting Director ONDQA/CDER/FDA 1 Outline Introduction to Quality by Design (QbD)
More informationAutomated Stability Testing Delivering Rapid Product Development
Automated Stability Testing Delivering Rapid Product Development Dr Andy Brookes, AstraZeneca APS PharmSci 2-4 th September 2013 With thanks to: Pam Harrison, Helen Porter, Faye Turner, Arun Kamath The
More informationASAP concept and case studies
ASAP concept and case studies Sabine Thielges Stability Testing for Pharmaceuticals 20-21 March 2013 Background: Traditional Approaches to Stability Studies (ICH) 1. Long-term e.g. 25 C/60%RH, 30 C/75%RH
More informationwhy open access publication of stability date is essential Mark Santillo Regional QA Officer SW England
why open access publication of stability date is essential Mark Santillo Regional QA Officer SW England } Standards for stability testing and data interpretation } Considerations for small molecule antimicrobials
More informationPQRI Stability Shelf-Life Working Group Glossary Ver 1 0.doc
Term Accelerated testing Acceptance criteria: Batch Definition Studies designed to increase the rate of chemical degradation or physical change of a drug substance or drug product by using exaggerated
More informationContents. Part I: Fundamentals of Bayesian Inference 1
Contents Preface xiii Part I: Fundamentals of Bayesian Inference 1 1 Probability and inference 3 1.1 The three steps of Bayesian data analysis 3 1.2 General notation for statistical inference 4 1.3 Bayesian
More informationTo control the consistency and quality
Establishing Acceptance Criteria for Analytical Methods Knowing how method performance impacts out-of-specification rates may improve quality risk management and product knowledge. To control the consistency
More informationBayesian concept for combined Phase 2a/b trials
Bayesian concept for combined Phase 2a/b trials /////////// Stefan Klein 07/12/2018 Agenda Phase 2a: PoC studies Phase 2b: dose finding studies Simulation Results / Discussion 2 /// Bayer /// Bayesian
More informationASAP: Accelerated Stability Assessment Program. Improved Protocol and Data Analysis for Accelerated Shelf-Life Estimation of Solid.
ASAP: Accelerated Stability Assessment Program Improved Protocol and Data Analysis for Accelerated Shelf-Life Estimation of Solid Dosage Forms Kenneth C. Waterman Research Fellow, Pfizer Inc A more rapid
More informationserve the goal of analytical lmethod Its data reveals the quality, reliability and consistency of
Analytical method validation By Juree Charoenteeraboon, Ph.D Analytical method Goal consistent, reliable and accurate data Validation analytical method serve the goal of analytical lmethod Its data reveals
More informationMean Kinetic Temperature
Mean Kinetic Temperature History, Math, and Applications Paul Daniel Sr. Regulatory Compliance Expert paul.daniel@vaisala.com Practical Webinar Stuff Challenges with your audio connection? Try calling
More informationStability of pharmaceutical preparations
Stability of pharmaceutical preparations ICH Guidelines Q1A(R2) Stability Testing of New Drug Substances and Products Q1B Stability Testing : Photostability Testing of New Drug Substances and Products
More informationUnified Approach For Performing An Analytical Methods Comparability Study IVT s LAB WEEK Presenter: Peter M. Saama, Ph.D. Bayer HealthCare LLC
Unified Approach For Performing An Analytical Methods Comparability Study IVT s LAB WEEK 2015 Presenter: Peter M. Saama, Ph.D. Bayer HealthCare LLC Agenda/ Content Page 1 Analytical Methods for Performing
More informationPhysical Pharmacy. Chemical Kinetics and Stability
Physical Pharmacy Chemical Kinetics and Stability Chemical Kinetics and Stability The purpose of stability testing is to provide evidence on how the quality of a drug substance or drug product varies with
More informationINTERNATIONAL CONFERENCE ON HARMONISATION OF TECHNICAL REQUIREMENTS FOR REGISTRATION OF PHARMACEUTICALS FOR HUMAN USE
INTERNATIONAL CONFERENCE ON HARMONISATION OF TECHNICAL REQUIREMENTS FOR REGISTRATION OF PHARMACEUTICALS FOR HUMAN USE ICH HARMONISED TRIPARTITE GUIDELINE TEXT ON VALIDATION OF ANALYTICAL PROCEDURES Recommended
More informationGuidance for Industry
Guidance for Industry M4: Organization of the CTD U.S. Department of Health and Human Services Food and Drug Administration Center for Drug Evaluation and Research (CDER) Center for Biologics Evaluation
More informationVALIDATION OF ANALYTICAL METHODS FOR PHARMACEUTICAL ANALYSIS
EXCERPT FROM: VALIDATION OF ANALYTICAL METHODS FOR PHARMACEUTICAL ANALYSIS BY OONA MCPOLIN This excerpt is provided as a preview to enable the reader to sample the content and style of this book. You can
More informationTesting for bioequivalence of highly variable drugs from TR-RT crossover designs with heterogeneous residual variances
Testing for bioequivalence of highly variable drugs from T-T crossover designs with heterogeneous residual variances Christopher I. Vahl, PhD Department of Statistics Kansas State University Qing Kang,
More informationIndependence and Dependence in Calibration: A Discussion FDA and EMA Guidelines
ENGINEERING RESEARCH CENTER FOR STRUCTURED ORGANIC PARTICULATE SYSTEMS Independence and Dependence in Calibration: A Discussion FDA and EMA Guidelines Rodolfo J. Romañach, Ph.D. ERC-SOPS Puerto Rico Site
More informationBayesian Modeling of Accelerated Life Tests with Random Effects
Bayesian Modeling of Accelerated Life Tests with Random Effects Ramón V. León Avery J. Ashby Jayanth Thyagarajan Joint Statistical Meeting August, 00 Toronto, Canada Abstract We show how to use Bayesian
More informationBayesian Networks in Educational Assessment
Bayesian Networks in Educational Assessment Estimating Parameters with MCMC Bayesian Inference: Expanding Our Context Roy Levy Arizona State University Roy.Levy@asu.edu 2017 Roy Levy MCMC 1 MCMC 2 Posterior
More informationSTUDY OF THE APPLICABILTY OF CONTENT UNIFORMITY AND DISSOLUTION VARIATION TEST ON ROPINIROLE HYDROCHLORIDE TABLETS
& STUDY OF THE APPLICABILTY OF CONTENT UNIFORMITY AND DISSOLUTION VARIATION TEST ON ROPINIROLE HYDROCHLORIDE TABLETS Edina Vranić¹*, Alija Uzunović² ¹ Department of Pharmaceutical Technology, Faculty of
More informationS.O.P. no : ARO/126 Title: Preparation and validation of lyophilised Page : 1 of 6 urinary reference materials. Revision : 1 Date :
1. INTRODUCTION CONCEPT Title: Preparation and validation of lyophilised Page : 1 of 6 Reference materials are samples with a known concentration of analytes. It s the assignment of the Community Reference
More informationDEMONSTRATING CAPABILITY TO COMPLY WITH A TEST PROCEDURE: THE CONTENT UNIFORMITY AND DISSOLUTION ACCEPTANCE LIMITS (CUDAL) APPROACH
1 DEMONSTRATING CAPABILITY TO COMPLY WITH A TEST PROCEDURE: THE CONTENT UNIFORMITY AND DISSOLUTION ACCEPTANCE LIMITS (CUDAL) APPROACH Jim Bergum September 12, 2011 Key Responses For Batch Release 2 Potency
More informationStat 500 Midterm 2 8 November 2007 page 0 of 4
Stat 500 Midterm 2 8 November 2007 page 0 of 4 Please put your name on the back of your answer book. Do NOT put it on the front. Thanks. DO NOT START until I tell you to. You are welcome to read this front
More informationCase Study in the Use of Bayesian Hierarchical Modeling and Simulation for Design and Analysis of a Clinical Trial
Case Study in the Use of Bayesian Hierarchical Modeling and Simulation for Design and Analysis of a Clinical Trial William R. Gillespie Pharsight Corporation Cary, North Carolina, USA PAGE 2003 Verona,
More informationMechanistic Modelling of Product Degradation
Mechanistic Modelling of Product Degradation Henrique Reis Sardinha Abstract Physical and chemical stability of new pharmaceutical products represent two crucial factors for the pharmaceutical industry.
More informationDegradation data analysis for samples under unequal operating conditions: a case study on train wheels
Degradation data analysis for samples under unequal operating conditions: a case study on train wheels Marta AFONSO FREITAS UNIVERSIDADE FEDERAL DE MINAS GERAIS DEPARTMENT OF PRODUCTION ENGINEERING marta.afreitas@gmail.com
More informationQuantile POD for Hit-Miss Data
Quantile POD for Hit-Miss Data Yew-Meng Koh a and William Q. Meeker a a Center for Nondestructive Evaluation, Department of Statistics, Iowa State niversity, Ames, Iowa 50010 Abstract. Probability of detection
More informationEstimation of reliability parameters from Experimental data (Parte 2) Prof. Enrico Zio
Estimation of reliability parameters from Experimental data (Parte 2) This lecture Life test (t 1,t 2,...,t n ) Estimate θ of f T t θ For example: λ of f T (t)= λe - λt Classical approach (frequentist
More informationSwiss Medic Training Sampling
Swiss Medic Training Sampling Paul Sexton Sampling Preparation for Sampling Representative Sample Re-sampling Sampling Part I What to sample? Why sample? Where to sample? Who performs sampling? How to
More informationOnline Appendix to: Marijuana on Main Street? Estimating Demand in Markets with Limited Access
Online Appendix to: Marijuana on Main Street? Estating Demand in Markets with Lited Access By Liana Jacobi and Michelle Sovinsky This appendix provides details on the estation methodology for various speci
More informationThe Application of the Accelerated Stability Assessment Program (ASAP) to Quality by Design (QbD) for Drug Product Stability
AAPS PharmSciTech, Vol. 12, No. 3, September 2011 ( # 2011) DOI: 10.1208/s12249-011-9657-3 Review Article Theme: Quality by Design: Case Studies and Scientific Foundations Guest Editors: Robin Bogner,
More informationECMA st Edition / June Test Method for the Estimation of the Archival Lifetime of Optical Media
ECMA-379 1 st Edition / June 2007 Test Method for the Estimation of the Archival Lifetime of Optical Media COPYRIGHT PROTECTED DOCUMENT Ecma International 2007 Standard ECMA-379 1 st Edition / June 2007
More informationEstimation and Confidence Intervals for Parameters of a Cumulative Damage Model
United States Department of Agriculture Forest Service Forest Products Laboratory Research Paper FPL-RP-484 Estimation and Confidence Intervals for Parameters of a Cumulative Damage Model Carol L. Link
More informationSetting Attainable and Practical Particle Size Specifications
Setting Attainable and Practical Particle Size Specifications Mark Bumiller HORIBA Scientific mark.bumiller@horiba.com Why Set Specifications? Product release (quality) Communicate grade to buyers Internal
More informationMuncie, IN, May 25, 2010
: Impact on Drug Product Design and Stability Test Design William R. Porter 33rd ANNUAL MIDWEST BIPHARMACEUTICAL STATISTICS WRKSHP Muncie, IN, Stability by Design (SbD): Instability Assessment A Tool for
More informationVALIDATION OF ANALYTICAL METHODS. Presented by Dr. A. Suneetha Dept. of Pharm. Analysis Hindu College of Pharmacy
VALIDATION OF ANALYTICAL METHODS Presented by Dr. A. Suneetha Dept. of Pharm. Analysis Hindu College of Pharmacy According to FDA,validation is established documented evidence which provides a high degree
More informationUsing Historical Experimental Information in the Bayesian Analysis of Reproduction Toxicological Experimental Results
Using Historical Experimental Information in the Bayesian Analysis of Reproduction Toxicological Experimental Results Jing Zhang Miami University August 12, 2014 Jing Zhang (Miami University) Using Historical
More informationEnquiry. Demonstration of Uniformity of Dosage Units using Large Sample Sizes. Proposal for a new general chapter in the European Pharmacopoeia
Enquiry Demonstration of Uniformity of Dosage Units using Large Sample Sizes Proposal for a new general chapter in the European Pharmacopoeia In order to take advantage of increased batch control offered
More informationReview of the role of uncertainties in room acoustics
Review of the role of uncertainties in room acoustics Ralph T. Muehleisen, Ph.D. PE, FASA, INCE Board Certified Principal Building Scientist and BEDTR Technical Lead Division of Decision and Information
More informationChapter 4: Verification of compendial methods
Chapter 4: Verification of compendial methods Introduction In order to ensure accurate and reliable test results, the quality control laboratory (QCL) needs to use analytical methods (and accompanying
More informationThe STS Surgeon Composite Technical Appendix
The STS Surgeon Composite Technical Appendix Overview Surgeon-specific risk-adjusted operative operative mortality and major complication rates were estimated using a bivariate random-effects logistic
More informationA Flexible Class of Models for Data Arising from a thorough QT/QTc study
A Flexible Class of Models for Data Arising from a thorough QT/QTc study Suraj P. Anand and Sujit K. Ghosh Department of Statistics, North Carolina State University, Raleigh, North Carolina, USA Institute
More informationMeasurement error as missing data: the case of epidemiologic assays. Roderick J. Little
Measurement error as missing data: the case of epidemiologic assays Roderick J. Little Outline Discuss two related calibration topics where classical methods are deficient (A) Limit of quantification methods
More informationAnalytical Method Validation: An Updated Review
INTERNATIONAL JOURNAL OF ADVANCES IN PHARMACY, BIOLOGY AND CHEMISTRY Analytical Method Validation: An Updated Review G. Geetha, Karanam Naga Ganika Raju, B. Vignesh Kumar and M. Gnana Raja* Sankaralingam
More informationStat 5421 Lecture Notes Proper Conjugate Priors for Exponential Families Charles J. Geyer March 28, 2016
Stat 5421 Lecture Notes Proper Conjugate Priors for Exponential Families Charles J. Geyer March 28, 2016 1 Theory This section explains the theory of conjugate priors for exponential families of distributions,
More informationHow to develop an ATP.... and verify it has been met. Melissa Hanna-Brown Analytical R & D; Pfizer Global R+D, Sandwich UK.
How to develop an ATP... and verify it has been met Melissa Hanna-Brown Analytical R & D; Pfizer Global R+D, Sandwich UK. 06 May 2014 Overview Some context accuracy, precision and risk What is an ATP?
More informationEstimation of Operational Risk Capital Charge under Parameter Uncertainty
Estimation of Operational Risk Capital Charge under Parameter Uncertainty Pavel V. Shevchenko Principal Research Scientist, CSIRO Mathematical and Information Sciences, Sydney, Locked Bag 17, North Ryde,
More informationEVALUATION OF NON-ACCELERATED STABILITY DATA. Kathleen Karpenter Dietrich and Daniel L. Weiner Merrell Dow Pharmaceuticals Inc. 1. Y ij. 2. Y a. + bx.
VALUATIO OF O-ACCLRATD STABILITY DATA Kathleen Karpenter Dietrich and Daniel L. Weiner Merrell Dow Pharmaceuticals Inc. Abstract A package of SAS macros has been written which allows users to evaluate
More informationChapter 13 - Chemical Kinetics II. Integrated Rate Laws Reaction Rates and Temperature
Chapter 13 - Chemical Kinetics II Integrated Rate Laws Reaction Rates and Temperature Reaction Order - Graphical Picture A ->Products Integrated Rate Laws Zero Order Reactions Rate = k[a] 0 = k (constant
More informationSubject CS1 Actuarial Statistics 1 Core Principles
Institute of Actuaries of India Subject CS1 Actuarial Statistics 1 Core Principles For 2019 Examinations Aim The aim of the Actuarial Statistics 1 subject is to provide a grounding in mathematical and
More informationR-squared for Bayesian regression models
R-squared for Bayesian regression models Andrew Gelman Ben Goodrich Jonah Gabry Imad Ali 8 Nov 2017 Abstract The usual definition of R 2 (variance of the predicted values divided by the variance of the
More informationDefault Priors and Effcient Posterior Computation in Bayesian
Default Priors and Effcient Posterior Computation in Bayesian Factor Analysis January 16, 2010 Presented by Eric Wang, Duke University Background and Motivation A Brief Review of Parameter Expansion Literature
More information(5) Multi-parameter models - Gibbs sampling. ST440/540: Applied Bayesian Analysis
Summarizing a posterior Given the data and prior the posterior is determined Summarizing the posterior gives parameter estimates, intervals, and hypothesis tests Most of these computations are integrals
More informationOnline publication date: 22 March 2010
This article was downloaded by: [South Dakota State University] On: 25 March 2010 Access details: Access Details: [subscription number 919556249] Publisher Taylor & Francis Informa Ltd Registered in England
More informationNew Statistical Methods That Improve on MLE and GLM Including for Reserve Modeling GARY G VENTER
New Statistical Methods That Improve on MLE and GLM Including for Reserve Modeling GARY G VENTER MLE Going the Way of the Buggy Whip Used to be gold standard of statistical estimation Minimum variance
More informationContents. Preface to Second Edition Preface to First Edition Abbreviations PART I PRINCIPLES OF STATISTICAL THINKING AND ANALYSIS 1
Contents Preface to Second Edition Preface to First Edition Abbreviations xv xvii xix PART I PRINCIPLES OF STATISTICAL THINKING AND ANALYSIS 1 1 The Role of Statistical Methods in Modern Industry and Services
More informationSample size considerations in IM assays
Sample size considerations in IM assays Lanju Zhang Senior Principal Statistician, Nonclinical Biostatistics Midwest Biopharmaceutical Statistics Workshop May 6, 0 Outline Introduction Homogeneous population
More informationThe Bayesian Approach to Multi-equation Econometric Model Estimation
Journal of Statistical and Econometric Methods, vol.3, no.1, 2014, 85-96 ISSN: 2241-0384 (print), 2241-0376 (online) Scienpress Ltd, 2014 The Bayesian Approach to Multi-equation Econometric Model Estimation
More informationROBUSTNESS OF TWO-PHASE REGRESSION TESTS
REVSTAT Statistical Journal Volume 3, Number 1, June 2005, 1 18 ROBUSTNESS OF TWO-PHASE REGRESSION TESTS Authors: Carlos A.R. Diniz Departamento de Estatística, Universidade Federal de São Carlos, São
More informationAssessing the Effect of Prior Distribution Assumption on the Variance Parameters in Evaluating Bioequivalence Trials
Georgia State University ScholarWorks @ Georgia State University Mathematics Theses Department of Mathematics and Statistics 8--006 Assessing the Effect of Prior Distribution Assumption on the Variance
More informationBasic Sampling Methods
Basic Sampling Methods Sargur Srihari srihari@cedar.buffalo.edu 1 1. Motivation Topics Intractability in ML How sampling can help 2. Ancestral Sampling Using BNs 3. Transforming a Uniform Distribution
More informationTraceability, validation and measurement uncertainty 3 pillars for quality of measurement results. David MILDE
Traceability, validation and measurement uncertainty 3 pillars for quality of measurement results David MILDE Traceability Validation Measurement Uncerainty Outline Introduction to 3 pillars Recommended
More informationBayesian Analysis. Bayesian Analysis: Bayesian methods concern one s belief about θ. [Current Belief (Posterior)] (Prior Belief) x (Data) Outline
Bayesian Analysis DuBois Bowman, Ph.D. Gordana Derado, M. S. Shuo Chen, M. S. Department of Biostatistics and Bioinformatics Center for Biomedical Imaging Statistics Emory University Outline I. Introduction
More informationNonparametric Bayesian Methods - Lecture I
Nonparametric Bayesian Methods - Lecture I Harry van Zanten Korteweg-de Vries Institute for Mathematics CRiSM Masterclass, April 4-6, 2016 Overview of the lectures I Intro to nonparametric Bayesian statistics
More informationReports of the Institute of Biostatistics
Reports of the Institute of Biostatistics No 02 / 2008 Leibniz University of Hannover Natural Sciences Faculty Title: Properties of confidence intervals for the comparison of small binomial proportions
More informationBasic Statistics. 1. Gross error analyst makes a gross mistake (misread balance or entered wrong value into calculation).
Basic Statistics There are three types of error: 1. Gross error analyst makes a gross mistake (misread balance or entered wrong value into calculation). 2. Systematic error - always too high or too low
More informationUncertainty and Graphical Analysis
Uncertainty and Graphical Analysis Introduction Two measures of the quality of an experimental result are its accuracy and its precision. An accurate result is consistent with some ideal, true value, perhaps
More informationPredictive Microbiology (theory)
(theory) Kostas Koutsoumanis Head of Lab of Food Microbiology and Hygiene Dpt. Of Food Science and Technology Aristotle University of Thessaloniki Summer School In Silico Methods for Food Safety THE CONCEPT
More informationGuidance for Industry M4: The CTD General Questions and Answers
Guidance for Industry M4: The CTD General Questions and Answers U.S. Department of Health and Human Services Food and Drug Administration Center for Drug Evaluation and Research (CDER) Center for Biologics
More information(QALT): Data Analysis and Test Design
(QALT): Data Analysis and Test Design Huairui(Harry) Guo, Ph.D, CRE, CQE Thanasis Gerokostopoulos, CRE, CQE ASTR 2013, Oct. 9-11, San Diego, CA Qualitative vs. Quantitative ALT (Accelerated Life Tests)
More informationDesign of Multiregional Clinical Trials (MRCT) Theory and Practice
Design of Multiregional Clinical Trials (MRCT) Theory and Practice Gordon Lan Janssen R&D, Johnson & Johnson BASS 2015, Rockville, Maryland Collaborators: Fei Chen and Gang Li, Johnson & Johnson BASS 2015_LAN
More informationAnalysis of Regression and Bayesian Predictive Uncertainty Measures
Analysis of and Predictive Uncertainty Measures Dan Lu, Mary C. Hill, Ming Ye Florida State University, dl7f@fsu.edu, mye@fsu.edu, Tallahassee, FL, USA U.S. Geological Survey, mchill@usgs.gov, Boulder,
More informationNew Horizons in Accelerated Stability Modeling-- Tablet Dissolution, Tablet Disintegration and Product Appearance
New Horizons in Accelerated Stability Modeling-- Tablet Dissolution, Tablet Disintegration and Product Appearance Kenneth C. Waterman, Ph.D. FreeThink Technologies, Inc. Branford, CT ken.waterman@freethinktech.com
More informationApproximating Bayesian Posterior Means Using Multivariate Gaussian Quadrature
Approximating Bayesian Posterior Means Using Multivariate Gaussian Quadrature John A.L. Cranfield Paul V. Preckel Songquan Liu Presented at Western Agricultural Economics Association 1997 Annual Meeting
More informationAccelerated Stability of Peptides. Alisa K. Waterman, Ph.D. FreeThink Technologies, Inc.
Accelerated Stability of Peptides Alisa K. Waterman, Ph.D. FreeThink Technologies, Inc. 1 Small Molecules Subject to chemical modifications such as oxidation, deamination, hydrolysis Chemical modifications
More informationForward Problems and their Inverse Solutions
Forward Problems and their Inverse Solutions Sarah Zedler 1,2 1 King Abdullah University of Science and Technology 2 University of Texas at Austin February, 2013 Outline 1 Forward Problem Example Weather
More informationComparison of Three Calculation Methods for a Bayesian Inference of Two Poisson Parameters
Journal of Modern Applied Statistical Methods Volume 13 Issue 1 Article 26 5-1-2014 Comparison of Three Calculation Methods for a Bayesian Inference of Two Poisson Parameters Yohei Kawasaki Tokyo University
More informationComputer Vision Group Prof. Daniel Cremers. 10a. Markov Chain Monte Carlo
Group Prof. Daniel Cremers 10a. Markov Chain Monte Carlo Markov Chain Monte Carlo In high-dimensional spaces, rejection sampling and importance sampling are very inefficient An alternative is Markov Chain
More informationProvläsningsexemplar / Preview INTERNATIONAL STANDARD ISO Second edition
INTERNATIONAL STANDARD ISO 6974-1 Second edition 2012-05-15 Natural gas Determination of composition and associated uncertainty by gas chromatography Part 1: General guidelines and calculation of composition
More informationBUGS Bayesian inference Using Gibbs Sampling
BUGS Bayesian inference Using Gibbs Sampling Glen DePalma Department of Statistics May 30, 2013 www.stat.purdue.edu/~gdepalma 1 / 20 Bayesian Philosophy I [Pearl] turned Bayesian in 1971, as soon as I
More informationRule of Thumb Think beyond simple ANOVA when a factor is time or dose think ANCOVA.
May 003: Think beyond simple ANOVA when a factor is time or dose think ANCOVA. Case B: Factorial ANOVA (New Rule, 6.3). A few corrections have been inserted in blue. [At times I encounter information that
More informationZero Order First Order Second Order. ln[a]t = k t + ln[a]0. [A]t = [A]0 e -kt 1 /[A]t = k t + ( 1 /[A]0 )
Kinetics (ch 14) Cañada College Chem 220 Name: Student ID: Zero Order First Order Second Order Dif Law Rate = k Rate = k [A] Rate = k [A] 2 Int Law [A]t = k t + [A]0 ln[a]t = k t + ln[a]0 [A]t = [A]0 e
More informationBayesian Methods in Multilevel Regression
Bayesian Methods in Multilevel Regression Joop Hox MuLOG, 15 september 2000 mcmc What is Statistics?! Statistics is about uncertainty To err is human, to forgive divine, but to include errors in your design
More informationValidation of analytical methods. Adrian Covaci Toxicological Center, University of Antwerp
Validation of analytical methods Adrian Covaci Toxicological Center, University of Antwerp Why validate? The objective of validating an analytical method is to demonstrate its suitability for its intended
More informationBAYESIAN ANALYSIS OF DOSE-RESPONSE CALIBRATION CURVES
Libraries Annual Conference on Applied Statistics in Agriculture 2005-17th Annual Conference Proceedings BAYESIAN ANALYSIS OF DOSE-RESPONSE CALIBRATION CURVES William J. Price Bahman Shafii Follow this
More informationCOMP90051 Statistical Machine Learning
COMP90051 Statistical Machine Learning Semester 2, 2017 Lecturer: Trevor Cohn 17. Bayesian inference; Bayesian regression Training == optimisation (?) Stages of learning & inference: Formulate model Regression
More informationThe Bayesian Choice. Christian P. Robert. From Decision-Theoretic Foundations to Computational Implementation. Second Edition.
Christian P. Robert The Bayesian Choice From Decision-Theoretic Foundations to Computational Implementation Second Edition With 23 Illustrations ^Springer" Contents Preface to the Second Edition Preface
More informationForecasting. Copyright 2015 Pearson Education, Inc.
5 Forecasting To accompany Quantitative Analysis for Management, Twelfth Edition, by Render, Stair, Hanna and Hale Power Point slides created by Jeff Heyl Copyright 2015 Pearson Education, Inc. LEARNING
More information