Robust design in model-based analysis of longitudinal clinical data
|
|
- Gervase Lucas
- 5 years ago
- Views:
Transcription
1 Robust design in model-based analysis of longitudinal clinical data Giulia Lestini, Sebastian Ueckert, France Mentré IAME UMR 1137, INSERM, University Paris Diderot, France PODE, June 0 016
2 Context Optimal design in Nonlinear Mixed Effects Models (NLMEM) Choosing a good design for a planned study is essential Number of subjects Number of sampling times for each subject Sampling times (allocation in time) Optimal design depends on prior information (model and parameters) Adaptive design Robust design (robustness on parameters) Atkinson, Optimum Experimental Designs Dodds, Hooker and Vicini, J Pharmacokinet Pharmacodyn. 005 Foo L-K, Duffull S. J Biopharm Stat. 010
3 Objectives Objectives To compare various robust design criteria in NLMEM for two examples: Pharmacokinetic/Pharmacodynamic (PKPD) model with continuous data Longitudinal binary model using a new method for the evaluation of the Fisher information matrix (FIM) for NLMEM with discrete data 3
4 Method Basic mixed effect model Individual model (one continuous response) y i =f(φ i, ξ i ) + ε i vector of n i observations ξ i : individual sampling times t ij j=1, n i φ i : individual parameters (size p) f: nonlinear function defining the structural model ε i : gaussian zero mean random error var (ε i ) = diag (σ inter + σ slope f(φ i,ξ i )) combined error model Random-effects model φ i = μ exp(b i ) or μ + b i b i ~N 0, Ω here Ω diagonal: ω k = Var(b ik ) Population parameters: Ψ (size P) μ (fixed effects) unknowns in Ω (variance of random effects) σ inter and/or σ slope (error model parameters) 4
5 Method Basic population design Assumption N individuals i same elementary design ξ in all N subjects (ξ i = ξ) with t 1,, t n sampling times n tot = N n Population design Ξ = ξ, N 5
6 Method Fisher Information Matrix (FIM) in NLMEM Elementary FIM: M F Ψ, ξ = E L(y;Ψ) Ψ Ψ T No analytical expression for FIM Continuous data FO approximation o M F is implemented in the R function «PFIM» 1 Discrete data New method: Adaptive Gaussian Quadrature (AGQ) and Quasi Random Monte Carlo,3 (QRMC) o M F is implemented in an R program. Population FIM for one group design M F Ψ, Ξ = N M F Ψ, ξ Standard D- criterion: M F Ψ, Ξ 1/P where Ψ = Ψ (fixed values) 1 Ueckert S, Mentré F. Computational and Methodological Statistics (CMStatistics) Ueckert S, Mentré F. Population Optimum Design of Experiments (PODE)
7 Method Criteria for optimal robust designs For robust design, a distribution for Ψ, p Ψ, is assumed Optimal robust designs ξ DE ξ ED ξ EID ξ ELD ξ MM = argmax ξ Criteria E Ψ M F Ψ, ξ E Ψ M F Ψ, ξ E Ψ M F Ψ, ξ E Ψ log M F Ψ, ξ min Ψ M F Ψ, ξ 1 1 7
8 Method Criteria for robust optimal designs Criteria are computed by Monte Carlo simulations (MC) K= Total number of MC samples Ψ 1... Ψ k... Ψ K M F Ψ 1, ξ... M F Ψ k, ξ... M F Ψ K, ξ Compute Robust Criteria 8
9 Outline Part I Part Comparison of robust design criteria in NLMEM with continuous data Designs were evaluated using (i) D-criterion and predicted Relative Standard Errors (RSE) (ii) Relative Root Mean Squared Errors (RRMSE) derived from Clinical Trial Simulations (CTS) Comparison of robust design criteria in NLMEM with discrete data (NEW way to compute FIM) Designs were evaluated using (i) D-criterion and predicted Relative Standard Errors (RSE) 9
10 Objectives Objectives To compare various robust design criteria in NLMEM for Pharmacokinetic/Pharmacodynamic (PKPD) model with continuous data 10
11 Method Part I - Example: PKPD model with continuous data responses model, for a biomarker in oncology (TGF β) PK: concentration f PK (φ, t) = DOSE V k = CL V k a k a k e kt e k at, ka Concentration CL/V Parameters: k a, V, CL PD: relative inhibition of TGF-β df PD (φ, t) dt = k out I max f PK (φ, t) f PK (φ, t)+ic 50 k out f PD (φ, t), ksyn Effect kout I max = 1 Parameters: k out, IC 50 The model is implemented in the DDMoRe model repository Gueorguieva et al., Comput Methods Programs Biomed. 007 Bueno et al., Eur J Cancer. 008 Gueorguieva et al., Br J Clin Pharmacol. 014 Lestini et al., Pharm Res
12 Method Part I - Example: PKPD model with continuous data PK Parameters Ψ* p(ψ) μ ka h 1 μ V L μ CL Lh ω V ω CL σ slope,pk PD Parameters Ψ* p(ψ) μ kout h 1 0. logn log(0.), 0.8 μ IC50 mg/l 0.3 logn log(0.3), 0.8 ω kout ω IC50 σ inter,pd
13 Method For each optimal design, and for a fixed equispaced design ξ ES compute D-criterion 13 and predicted RSE(%) for each set simulated parameters Ψ k, k = 1,.., K Design optimization Constraints N = 50 patients n = 3 observations per patient For PK, times fixed to 0.1, 4, 1 h For PD, 3 sampling times among possible times: 1,, 3, 4, 6, 9,10,15,, 3, 4 h 11 3 = 165 elementary designs FIM Computation FIM was computed in PFIM with FO approximation Optimization Find the D-optimal design ξ D using D-criterion for Ψ* Find optimal robust designs ξ DE, ξ ED, ξ EID, ξ ELD, ξ MM using MC with K=1000 Evaluation
14 Results D-criterion for optimal designs Robust designs Non-Robust designs Best design with highest median: ξ ELD Whiskers of boxplot: 10 th and 90 th percentiles 14
15 Results Predicted Relative Standard Errors (RSE) μ kout μ IC 50 Robust designs Non-Robust designs Whiskers of boxplot: 10 th and 90 th percentiles 15
16 Results Predicted Relative Standard Errors (RSE) ω k out ω IC 50 Robust designs Non-Robust designs Whiskers of boxplot: 10 th and 90 th percentiles 16
17 Method Simulation Study For each design ξ Ψ 1 Y 1 Ψ 1... Ψ k... K=1000 Simulations N=50 subjects... Y k Ψ k... Compute Relative Root Mean Squared Error (RRMSE) Ψ K Y K Ψ K ξ = ξ DE, ξ ED, ξ EID, ξ ELD, ξ MM, ξ D, ξ ES Parameter estimation: SAEM algorithm in MONOLIX chains, initial estimates: Ψ* 17
18 Results Relative Root Mean Squared Errors (RRMSE) = Means across parameters of RRMSEs standardized to ξ ELD RRMSEs 18
19 Conclusion Conclusion (Part 1) All criteria led to various designs rather different ξ ELD performed globally the best in terms of median of D-criteria across the 1000 MC simulations RRMSE obtained from the CTS study confirm the results, showing ξ ELD being globally the more robust design of the 1000 simulations 19
20 Objectives Objectives To compare various robust design criteria in NLMEM for Longitudinal binary model using a new method for the evaluation of the Fisher information matrix (FIM) for NLMEM with discrete data 0
21 Method Part - Example: NLMEM with binary data Logistic model for repeated binary response with treatment increasing the slope of the logit of the response with time logit π = β 1 + β (1 + μ 3 δ)t where β = g μ, b = μ + b; π is the probability of success treatment groups (δ = 0 & δ = 1) Parameters Ψ* p(ψ) μ μ (months) 0.09 N 0.09, 0. μ 3 5 N 5, ω (months) ω 1 Ogungbenro K, Aarons L.. J Pharmacokinet Pharmacodyn. 011 Ueckert S, Mentré F. Computational and Methodological Statistics (CMStatistics)
22 Method Design optimization Constraints N = N T (treatment) + N C (control) = 100 patients n = 4 observations per patient including 0 and 1 months Possible intermediate times: 1,, 3, 4, 5, 6, 7, 8, 9, 10, 11 months 11 = 55 elementary designs (first and last times are fixed) FIM Computation o FIM was computed with the new method by Ueckert and Mentré, 015, based on AGQ and QRMC, with 3 nodes and 500 integrations samples Optimization Find the D-optimal design ξ D using D-criterion for Ψ* Find optimal robust designs ξ DE, ξ ED, ξ EID, ξ ELD, ξ MM using MC with K=1000 Evaluation For each optimal design, and for a fixed equispaced design ξ ES compute D- criterion and predicted RSE(%) for each set simulated parameters Ψ k, k = 1,.., K
23 Results D-criterion distribution for selected designs Robust designs Non-Robust designs Best designs with highest median: ξ EID and ξ ELD Whiskers of boxplot: 10 th and 90 th percentiles
24 Results Predicted Relative Standard Errors (RSE) Robust designs Non-Robust designs μ 1 μ μ 3 Whiskers of boxplot: 10 th and 90 th percentiles 4
25 Results Predicted Relative Standard Errors (RSE) ω 1 ω Robust designs Non-Robust designs Whiskers of boxplot: 10 th and 90 th percentiles 5
26 Conclusion Conclusion (Part ) The evaluation of FIM with the new approach is rather fast, allowing for the first time robust design optimization for discrete longitudinal models From median of 1000 simulated D-criteria, ξ ELD and ξ EID performed globally the best ξ ES has efficiency of 0.8 compare to ξ D. When prior uncertainty is assumed, the efficiency is 0.66 compared to ξ ELD 6
27 Conclusion General conclusions All these robust criteria were never systematically compared in NLMEM Different criteria led to various optimal designs, with different impact on predicted RSE From median of 1000 simulated D-criteria For PKPD model: ELD led to the best designs For longitudinal binary model: ELD and EID led to the best designs (which are very close) Perspectives Perform CTS for binary data example (Part ) Perform robust and adaptive design Perform model averaging 7
28 Thank You! Thank you for your attention! The research leading to these results has received support from the Innovative Medicines Initiative Joint Undertaking under grant agreement n , resources of which are composed of financial contributions from the European Union's Seventh Framework Programme (FP7/ ) and EFPIA companies in kind contribution. The DDMoRe project is also financially supported by contributions from Academic and SME partners 8
29 Back-up 9
30 Method Criteria for robust optimal designs For robust design, a distribution for Ψ, p Ψ, is assumed Optimal designs Criteria Compute Criteria ξ DE E Ψ M F Ψ, ξ 1 K K k=1 M F Ψ k, ξ 1 K ξ ED E Ψ M F Ψ, ξ K M F Ψ k, ξ k=1 = argmax ξ EID ξ E Ψ M F Ψ, ξ K K M F Ψ k, ξ 1 k=1 ξ ELD E Ψ log M F Ψ, ξ 1 K K log M F Ψ k, ξ k=1 ξ MM min Ψ M F Ψ, ξ min k(k=1, K) M F Ψ k, ξ 1 30
Evaluation of the Fisher information matrix in nonlinear mixed eect models without linearization
Evaluation of the Fisher information matrix in nonlinear mixed eect models without linearization Sebastian Ueckert, Marie-Karelle Riviere, and France Mentré INSERM, IAME, UMR 1137, F-75018 Paris, France;
More informationF. Combes (1,2,3) S. Retout (2), N. Frey (2) and F. Mentré (1) PODE 2012
Prediction of shrinkage of individual parameters using the Bayesian information matrix in nonlinear mixed-effect models with application in pharmacokinetics F. Combes (1,2,3) S. Retout (2), N. Frey (2)
More informationDESIGN EVALUATION AND OPTIMISATION IN CROSSOVER PHARMACOKINETIC STUDIES ANALYSED BY NONLINEAR MIXED EFFECTS MODELS
DESIGN EVALUATION AND OPTIMISATION IN CROSSOVER PHARMACOKINETIC STUDIES ANALYSED BY NONLINEAR MIXED EFFECTS MODELS Thu Thuy Nguyen, Caroline Bazzoli, France Mentré UMR 738 INSERM - University Paris Diderot,
More informationSebastian Ueckert, France Mentré. To cite this version: HAL Id: inserm
A new method for evaluation of the Fisher information matrix for discrete mixed effect models using Monte Carlo sampling and adaptive Gaussian quadrature Sebastian Ueckert, France Mentré To cite this version:
More informationEvaluation of the Fisher information matrix in nonlinear mixed effect models using adaptive Gaussian quadrature
Evaluation of the Fisher information matrix in nonlinear mixed effect models using adaptive Gaussian quadrature Thu Thuy Nguyen, France Mentré To cite this version: Thu Thuy Nguyen, France Mentré. Evaluation
More informationStochastic approximation EM algorithm in nonlinear mixed effects model for viral load decrease during anti-hiv treatment
Stochastic approximation EM algorithm in nonlinear mixed effects model for viral load decrease during anti-hiv treatment Adeline Samson 1, Marc Lavielle and France Mentré 1 1 INSERM E0357, Department of
More informationAn Approach for Identifiability of Population Pharmacokinetic-Pharmacodynamic Models
An Approach for Identifiability of Population Pharmacokinetic-Pharmacodynamic Models Vittal Shivva 1 * Julia Korell 12 Ian Tucker 1 Stephen Duffull 1 1 School of Pharmacy University of Otago Dunedin New
More informationPopulation Design in Nonlinear Mixed Effects Multiple Response Models: extension of PFIM and evaluation by simulation with NONMEM and MONOLIX
Populaton Desgn n Nonlnear Mxed Effects Multple Response Models: extenson of PFIM and evaluaton by smulaton wth NONMEM and MONOLIX May 4th 007 Carolne Bazzol, Sylve Retout, France Mentré Inserm U738 Unversty
More informationOPTIMAL DESIGN FOR POPULATION PK/PD MODELS. 1. Introduction
ØÑ Å Ø Ñ Ø Ð ÈÙ Ð Ø ÓÒ DOI: 10.2478/v10127-012-0012-1 Tatra Mt. Math. Publ. 51 (2012), 115 130 OPTIMAL DESIGN FOR POPULATION PK/PD MODELS Sergei Leonov Alexander Aliev ABSTRACT. We provide some details
More informationAn effective approach for obtaining optimal sampling windows for population pharmacokinetic experiments
An effective approach for obtaining optimal sampling windows for population pharmacokinetic experiments Kayode Ogungbenro and Leon Aarons Centre for Applied Pharmacokinetic Research School of Pharmacy
More informationEstimation and Model Selection in Mixed Effects Models Part I. Adeline Samson 1
Estimation and Model Selection in Mixed Effects Models Part I Adeline Samson 1 1 University Paris Descartes Summer school 2009 - Lipari, Italy These slides are based on Marc Lavielle s slides Outline 1
More informationExtension of the SAEM algorithm for nonlinear mixed models with 2 levels of random effects Panhard Xavière
Extension of the SAEM algorithm for nonlinear mixed models with 2 levels of random effects Panhard Xavière * Modèles et mé thodes de l'évaluation thérapeutique des maladies chroniques INSERM : U738, Universit
More informationGeneralizing the MCPMod methodology beyond normal, independent data
Generalizing the MCPMod methodology beyond normal, independent data José Pinheiro Joint work with Frank Bretz and Björn Bornkamp Novartis AG Trends and Innovations in Clinical Trial Statistics Conference
More informationDescription of UseCase models in MDL
Description of UseCase models in MDL A set of UseCases (UCs) has been prepared to illustrate how different modelling features can be implemented in the Model Description Language or MDL. These UseCases
More informationHandling parameter constraints in complex/mechanistic population pharmacokinetic models
Handling parameter constraints in complex/mechanistic population pharmacokinetic models An application of the multivariate logistic-normal distribution Nikolaos Tsamandouras Centre for Applied Pharmacokinetic
More informationCorrection of the likelihood function as an alternative for imputing missing covariates. Wojciech Krzyzanski and An Vermeulen PAGE 2017 Budapest
Correction of the likelihood function as an alternative for imputing missing covariates Wojciech Krzyzanski and An Vermeulen PAGE 2017 Budapest 1 Covariates in Population PKPD Analysis Covariates are defined
More informationModelling a complex input process in a population pharmacokinetic analysis: example of mavoglurant oral absorption in healthy subjects
Modelling a complex input process in a population pharmacokinetic analysis: example of mavoglurant oral absorption in healthy subjects Thierry Wendling Manchester Pharmacy School Novartis Institutes for
More informationDesign of preclinical combination studies
University of Manchester, UK Talk Outline 1 Preliminaries What is a drug combination? Data structures. Ray design 2 Combination index Sources of variability 3 Two drugs Analysis of data 4 Types of studies
More informationHeterogeneous shedding of influenza by human subjects and. its implications for epidemiology and control
1 2 3 4 5 6 7 8 9 10 11 Heterogeneous shedding of influenza by human subjects and its implications for epidemiology and control Laetitia Canini 1*, Mark EJ Woolhouse 1, Taronna R. Maines 2, Fabrice Carrat
More informationSensitivity study of dose-finding methods
to of dose-finding methods Sarah Zohar 1 John O Quigley 2 1. Inserm, UMRS 717,Biostatistic Department, Hôpital Saint-Louis, Paris, France 2. Inserm, Université Paris VI, Paris, France. NY 2009 1 / 21 to
More informationGeneralizing the MCPMod methodology beyond normal, independent data
Generalizing the MCPMod methodology beyond normal, independent data José Pinheiro Joint work with Frank Bretz and Björn Bornkamp Novartis AG ASA NJ Chapter 35 th Annual Spring Symposium June 06, 2014 Outline
More informationExtension of the SAEM algorithm for nonlinear mixed. models with two levels of random effects
Author manuscript, published in "Biostatistics 29;1(1):121-35" DOI : 1.193/biostatistics/kxn2 Extension of the SAEM algorithm for nonlinear mixed models with two levels of random effects Xavière Panhard
More informationApproximate analysis of covariance in trials in rare diseases, in particular rare cancers
Approximate analysis of covariance in trials in rare diseases, in particular rare cancers Stephen Senn (c) Stephen Senn 1 Acknowledgements This work is partly supported by the European Union s 7th Framework
More informationA Clinical Trial Simulation System, Its Applications, and Future Challenges. Peter Westfall, Texas Tech University Kuenhi Tsai, Merck Research Lab
A Clinical Trial Simulation System, Its Applications, and Future Challenges Peter Westfall, Texas Tech University Kuenhi Tsai, Merck Research Lab Acknowledgement Major collaborators Stephan Ogenstand,
More informationComparison of multiple imputation methods for systematically and sporadically missing multilevel data
Comparison of multiple imputation methods for systematically and sporadically missing multilevel data V. Audigier, I. White, S. Jolani, T. Debray, M. Quartagno, J. Carpenter, S. van Buuren, M. Resche-Rigon
More informationExtension of the SAEM algorithm to left-censored data in nonlinear mixed-effects model: application to HIV dynamics model
Extension of the SAEM algorithm to left-censored data in nonlinear mixed-effects model: application to HIV dynamics model Adeline Samson 1, Marc Lavielle 2, France Mentré 1 1 INSERM U738, Paris, France;
More informationDesign of HIV Dynamic Experiments: A Case Study
Design of HIV Dynamic Experiments: A Case Study Cong Han Department of Biostatistics University of Washington Kathryn Chaloner Department of Biostatistics University of Iowa Nonlinear Mixed-Effects Models
More informationA comparison of estimation methods in nonlinear mixed effects models using a blind analysis
A comparison of estimation methods in nonlinear mixed effects models using a blind analysis Pascal Girard, PhD EA 3738 CTO, INSERM, University Claude Bernard Lyon I, Lyon France Mentré, PhD, MD Dept Biostatistics,
More informationImproved conditional approximations of the population Fisher Information Matrix
Improved condtonal approxmatons of the populaton Fsher Informaton Matrx Joakm Nyberg, Sebastan Ueckert and Andrew C. Hooker Pharmacometrcs research group Department of Pharmaceutcal Boscences Uppsala Unversty
More informationRandomized dose-escalation design for drug combination cancer trials with immunotherapy
Randomized dose-escalation design for drug combination cancer trials with immunotherapy Pavel Mozgunov, Thomas Jaki, Xavier Paoletti Medical and Pharmaceutical Statistics Research Unit, Department of Mathematics
More informationApproximation of the Fisher information and design in nonlinear mixed effects models
of the Fisher information and design in nonlinear mixed effects models Tobias Mielke Optimum for Mixed Effects Non-Linear and Generalized Linear PODE - 2011 Outline 1 2 3 Mixed Effects Similar functions
More informationEvaluation of bootstrap methods for estimating uncertainty of parameters in nonlinear mixedeffects
Evaluation of bootstrap methods for estimating uncertainty of parameters in nonlinear mixed-effects models: a simulation study in population pharmacokinetics Hoai-Thu Thai, France Mentré, Nick Holford,
More informationPharmacometrics : Nonlinear mixed effect models in Statistics. Department of Statistics Ewha Womans University Eun-Kyung Lee
Pharmacometrics : Nonlinear mixed effect models in Statistics Department of Statistics Ewha Womans University Eun-Kyung Lee 1 Clinical Trials Preclinical trial: animal study Phase I trial: First in human,
More information1Non Linear mixed effects ordinary differential equations models. M. Prague - SISTM - NLME-ODE September 27,
GDR MaMoVi 2017 Parameter estimation in Models with Random effects based on Ordinary Differential Equations: A bayesian maximum a posteriori approach. Mélanie PRAGUE, Daniel COMMENGES & Rodolphe THIÉBAUT
More informationImproving a safety of the Continual Reassessment Method via a modified allocation rule
Improving a safety of the Continual Reassessment Method via a modified allocation rule Pavel Mozgunov, Thomas Jaki Medical and Pharmaceutical Statistics Research Unit, Department of Mathematics and Statistics,
More informationEfficient algorithms for calculating optimal designs in pharmacokinetics and dose finding studies
Efficient algorithms for calculating optimal designs in pharmacokinetics and dose finding studies Tim Holland-Letz Ruhr-Universität Bochum Medizinische Fakultät 44780 Bochum, Germany email: tim.holland-letz@rub.de
More informationIntegration of SAS and NONMEM for Automation of Population Pharmacokinetic/Pharmacodynamic Modeling on UNIX systems
Integration of SAS and NONMEM for Automation of Population Pharmacokinetic/Pharmacodynamic Modeling on UNIX systems Alan J Xiao, Cognigen Corporation, Buffalo NY Jill B Fiedler-Kelly, Cognigen Corporation,
More informationJoint modeling of longitudinal and repeated time-to-event data using nonlinear mixed-effects models and the SAEM algorithm
Joint modeling of longitudinal and repeated time-to-event data using nonlinear mixed-effects models and the SAEM algorithm Cyprien Mbogning, Kevin Bleakley, Marc Lavielle To cite this version: Cyprien
More informationGroup sequential designs with negative binomial data
Group sequential designs with negative binomial data Ekkehard Glimm 1 Tobias Mütze 2,3 1 Statistical Methodology, Novartis, Basel, Switzerland 2 Department of Medical Statistics, University Medical Center
More informationHierarchical expectation propagation for Bayesian aggregation of average data
Hierarchical expectation propagation for Bayesian aggregation of average data Andrew Gelman, Columbia University Sebastian Weber, Novartis also Bob Carpenter, Daniel Lee, Frédéric Bois, Aki Vehtari, and
More informationAP-Optimum Designs for Minimizing the Average Variance and Probability-Based Optimality
AP-Optimum Designs for Minimizing the Average Variance and Probability-Based Optimality Authors: N. M. Kilany Faculty of Science, Menoufia University Menoufia, Egypt. (neveenkilany@hotmail.com) W. A. Hassanein
More informationOptimal Adaptive Designs for Dose Finding in Early Phase Clinical Trials
Optimal Adaptive Designs for Dose Finding in Early Phase Clinical Trials Alam, Muhammad Iftakhar The copyright of this thesis rests with the author and no quotation from it or information derived from
More informationCompound Optimal Designs for Percentile Estimation in Dose-Response Models with Restricted Design Intervals
Compound Optimal Designs for Percentile Estimation in Dose-Response Models with Restricted Design Intervals Stefanie Biedermann 1, Holger Dette 1, Wei Zhu 2 Abstract In dose-response studies, the dose
More informationPersonalized Treatment Selection Based on Randomized Clinical Trials. Tianxi Cai Department of Biostatistics Harvard School of Public Health
Personalized Treatment Selection Based on Randomized Clinical Trials Tianxi Cai Department of Biostatistics Harvard School of Public Health Outline Motivation A systematic approach to separating subpopulations
More informationWeb-based Supplementary Material for A Two-Part Joint. Model for the Analysis of Survival and Longitudinal Binary. Data with excess Zeros
Web-based Supplementary Material for A Two-Part Joint Model for the Analysis of Survival and Longitudinal Binary Data with excess Zeros Dimitris Rizopoulos, 1 Geert Verbeke, 1 Emmanuel Lesaffre 1 and Yves
More informationPractical issues related to the use of biomarkers in a seamless Phase II/III design
Practical issues related to the use of biomarkers in a seamless Phase II/III design Tomasz Burzykowski International Drug Development Institute Louvain-la-Neuve, Belgium tomasz.burzykowski@iddi.be www.iddi.com
More informationNonlinear Mixed Effects Models
Nonlinear Mixed Effects Modeling Department of Mathematics Center for Research in Scientific Computation Center for Quantitative Sciences in Biomedicine North Carolina State University July 31, 216 Introduction
More informationUncertainty quantification and visualization for functional random variables
Uncertainty quantification and visualization for functional random variables MascotNum Workshop 2014 S. Nanty 1,3 C. Helbert 2 A. Marrel 1 N. Pérot 1 C. Prieur 3 1 CEA, DEN/DER/SESI/LSMR, F-13108, Saint-Paul-lez-Durance,
More informationResearch Article. Multiple Imputation of Missing Covariates in NONMEM and Evaluation of the Method s Sensitivity to η-shrinkage
The AAPS Journal, Vol. 15, No. 4, October 2013 ( # 2013) DOI: 10.1208/s12248-013-9508-0 Research Article Multiple Imputation of Missing Covariates in NONMEM and Evaluation of the Method s Sensitivity to
More informationConfidence and Prediction Intervals for Pharmacometric Models
Citation: CPT Pharmacometrics Syst. Pharmacol. (218) 7, 36 373; VC 218 ASCPT All rights reserved doi:1.12/psp4.12286 TUTORIAL Confidence and Prediction Intervals for Pharmacometric Models Anne K ummel
More informationInterim Monitoring of Clinical Trials: Decision Theory, Dynamic Programming. and Optimal Stopping
Interim Monitoring of Clinical Trials: Decision Theory, Dynamic Programming and Optimal Stopping Christopher Jennison Department of Mathematical Sciences, University of Bath, UK http://people.bath.ac.uk/mascj
More informationRANDOMIZATIONN METHODS THAT
RANDOMIZATIONN METHODS THAT DEPEND ON THE COVARIATES WORK BY ALESSANDRO BALDI ANTOGNINI MAROUSSA ZAGORAIOU ALESSANDRA G GIOVAGNOLI (*) 1 DEPARTMENT OF STATISTICAL SCIENCES UNIVERSITY OF BOLOGNA, ITALY
More information3 Results. Part I. 3.1 Base/primary model
3 Results Part I 3.1 Base/primary model For the development of the base/primary population model the development dataset (for data details see Table 5 and sections 2.1 and 2.2), which included 1256 serum
More informationOptimal input design for nonlinear dynamical systems: a graph-theory approach
Optimal input design for nonlinear dynamical systems: a graph-theory approach Patricio E. Valenzuela Department of Automatic Control and ACCESS Linnaeus Centre KTH Royal Institute of Technology, Stockholm,
More informationOptimization Tools in an Uncertain Environment
Optimization Tools in an Uncertain Environment Michael C. Ferris University of Wisconsin, Madison Uncertainty Workshop, Chicago: July 21, 2008 Michael Ferris (University of Wisconsin) Stochastic optimization
More informationA weighted differential entropy based approach for dose-escalation trials
A weighted differential entropy based approach for dose-escalation trials Pavel Mozgunov, Thomas Jaki Medical and Pharmaceutical Statistics Research Unit, Department of Mathematics and Statistics, Lancaster
More informationOptimising Group Sequential Designs. Decision Theory, Dynamic Programming. and Optimal Stopping
: Decision Theory, Dynamic Programming and Optimal Stopping Christopher Jennison Department of Mathematical Sciences, University of Bath, UK http://people.bath.ac.uk/mascj InSPiRe Conference on Methodology
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 informationTHE IMPORTANCE OF THE SIMULATION EXPECTATION AS A GOODNESS OF FIT DIAGNOSTIC FOR CATEGORICAL POPULATION PHARMACODYNAMIC MODELS
THE IMPORTANCE OF THE SIMULATION EXPECTATION AS A GOODNESS OF FIT DIAGNOSTIC FOR CATEGORICAL POPULATION PHARMACODYNAMIC MODELS Marc R. Gastonguay, Jeffrey T. Hane Metrum Research Group LLC, Avon, CT 06001
More informationOn Assessing Bioequivalence and Interchangeability between Generics Based on Indirect Comparisons
On Assessing Bioequivalence and Interchangeability between Generics Based on Indirect Comparisons Jiayin Zheng 1, Shein-Chung Chow 1 and Mengdie Yuan 2 1 Department of Biostatistics & Bioinformatics, Duke
More informationType I error rate control in adaptive designs for confirmatory clinical trials with treatment selection at interim
Type I error rate control in adaptive designs for confirmatory clinical trials with treatment selection at interim Frank Bretz Statistical Methodology, Novartis Joint work with Martin Posch (Medical University
More informationA Note on Cochran Test for Homogeneity in Two Ways ANOVA and Meta-Analysis
Open Journal of Statistics, 05, 5, 787-796 Published Online December 05 in SciRes. http://www.scirp.org/journal/ojs http://dx.doi.org/0.436/ojs.05.57078 A Note on Cochran Test for Homogeneity in Two Ways
More informationSTA 4273H: Statistical Machine Learning
STA 4273H: Statistical Machine Learning Russ Salakhutdinov Department of Statistics! rsalakhu@utstat.toronto.edu! http://www.utstat.utoronto.ca/~rsalakhu/ Sidney Smith Hall, Room 6002 Lecture 3 Linear
More informationBACKGROUNDS. Two Models of Deformable Body. Distinct Element Method (DEM)
BACKGROUNDS Two Models of Deformable Body continuum rigid-body spring deformation expressed in terms of field variables assembly of rigid-bodies connected by spring Distinct Element Method (DEM) simple
More informationBayesian variable selection and classification with control of predictive values
Bayesian variable selection and classification with control of predictive values Eleni Vradi 1, Thomas Jaki 2, Richardus Vonk 1, Werner Brannath 3 1 Bayer AG, Germany, 2 Lancaster University, UK, 3 University
More informationAccurate Maximum Likelihood Estimation for Parametric Population Analysis. Bob Leary UCSD/SDSC and LAPK, USC School of Medicine
Accurate Maximum Likelihood Estimation for Parametric Population Analysis Bob Leary UCSD/SDSC and LAPK, USC School of Medicine Why use parametric maximum likelihood estimators? Consistency: θˆ θ as N ML
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 informationChecking model assumptions with regression diagnostics
@graemeleehickey www.glhickey.com graeme.hickey@liverpool.ac.uk Checking model assumptions with regression diagnostics Graeme L. Hickey University of Liverpool Conflicts of interest None Assistant Editor
More informationRandomized Quasi-Monte Carlo for MCMC
Randomized Quasi-Monte Carlo for MCMC Radu Craiu 1 Christiane Lemieux 2 1 Department of Statistics, Toronto 2 Department of Statistics, Waterloo Third Workshop on Monte Carlo Methods Harvard, May 2007
More informationPrerequisite: STATS 7 or STATS 8 or AP90 or (STATS 120A and STATS 120B and STATS 120C). AP90 with a minimum score of 3
University of California, Irvine 2017-2018 1 Statistics (STATS) Courses STATS 5. Seminar in Data Science. 1 Unit. An introduction to the field of Data Science; intended for entering freshman and transfers.
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 informationFitting Complex PK/PD Models with Stan
Fitting Complex PK/PD Models with Stan Michael Betancourt @betanalpha, @mcmc_stan University of Warwick Bayes Pharma Novartis Basel Campus May 20, 2015 Because of the clinical applications, precise PK/PD
More informationBIOSTATISTICAL METHODS
BIOSTATISTICAL METHODS FOR TRANSLATIONAL & CLINICAL RESEARCH Cross-over Designs #: DESIGNING CLINICAL RESEARCH The subtraction of measurements from the same subject will mostly cancel or minimize effects
More informationASME 2013 IDETC/CIE 2013 Paper number: DETC A DESIGN ORIENTED RELIABILITY METHODOLOGY FOR FATIGUE LIFE UNDER STOCHASTIC LOADINGS
ASME 2013 IDETC/CIE 2013 Paper number: DETC2013-12033 A DESIGN ORIENTED RELIABILITY METHODOLOGY FOR FATIGUE LIFE UNDER STOCHASTIC LOADINGS Zhen Hu, Xiaoping Du Department of Mechanical & Aerospace Engineering
More informationUncertainty Propagation
Setting: Uncertainty Propagation We assume that we have determined distributions for parameters e.g., Bayesian inference, prior experiments, expert opinion Ṫ 1 = 1 - d 1 T 1 - (1 - ")k 1 VT 1 Ṫ 2 = 2 -
More informationOptimal Designs for Some Selected Nonlinear Models
Optimal Designs for Some Selected Nonlinear Models A.S. Hedayat,, Ying Zhou Min Yang Department of Mathematics, Statistics, and Computer Science University of Illinois at Chicago Science and Engineering
More informationSAEMIX, an R version of the SAEM algorithm for parameter estimation in nonlinear mixed effect models
SAEMIX, an R version of the SAEM algorithm for parameter estimation in nonlinear mixed effect models Emmanuelle Comets 1,2 & Audrey Lavenu 2 & Marc Lavielle 3 1 CIC 0203, U. Rennes-I, CHU Pontchaillou,
More informationPractice of SAS Logistic Regression on Binary Pharmacodynamic Data Problems and Solutions. Alan J Xiao, Cognigen Corporation, Buffalo NY
Practice of SAS Logistic Regression on Binary Pharmacodynamic Data Problems and Solutions Alan J Xiao, Cognigen Corporation, Buffalo NY ABSTRACT Logistic regression has been widely applied to population
More informationA comparison of different Bayesian design criteria to compute efficient conjoint choice experiments
Faculty of Business and Economics A comparison of different Bayesian design criteria to compute efficient conjoint choice experiments J. Yu, P. Goos and M. Vandebroek DEPARTMENT OF DECISION SCIENCES AND
More informationSample Size Calculations for Group Randomized Trials with Unequal Sample Sizes through Monte Carlo Simulations
Sample Size Calculations for Group Randomized Trials with Unequal Sample Sizes through Monte Carlo Simulations Ben Brewer Duke University March 10, 2017 Introduction Group randomized trials (GRTs) are
More informationADAPTIVE EXPERIMENTAL DESIGNS. Maciej Patan and Barbara Bogacka. University of Zielona Góra, Poland and Queen Mary, University of London
ADAPTIVE EXPERIMENTAL DESIGNS FOR SIMULTANEOUS PK AND DOSE-SELECTION STUDIES IN PHASE I CLINICAL TRIALS Maciej Patan and Barbara Bogacka University of Zielona Góra, Poland and Queen Mary, University of
More informationChapter 4 Multi-factor Treatment Designs with Multiple Error Terms 93
Contents Preface ix Chapter 1 Introduction 1 1.1 Types of Models That Produce Data 1 1.2 Statistical Models 2 1.3 Fixed and Random Effects 4 1.4 Mixed Models 6 1.5 Typical Studies and the Modeling Issues
More informationImplementing Response-Adaptive Randomization in Multi-Armed Survival Trials
Implementing Response-Adaptive Randomization in Multi-Armed Survival Trials BASS Conference 2009 Alex Sverdlov, Bristol-Myers Squibb A.Sverdlov (B-MS) Response-Adaptive Randomization 1 / 35 Joint work
More informationMixture modelling of recurrent event times with long-term survivors: Analysis of Hutterite birth intervals. John W. Mac McDonald & Alessandro Rosina
Mixture modelling of recurrent event times with long-term survivors: Analysis of Hutterite birth intervals John W. Mac McDonald & Alessandro Rosina Quantitative Methods in the Social Sciences Seminar -
More informationOn the efficiency of two-stage adaptive designs
On the efficiency of two-stage adaptive designs Björn Bornkamp (Novartis Pharma AG) Based on: Dette, H., Bornkamp, B. and Bretz F. (2010): On the efficiency of adaptive designs www.statistik.tu-dortmund.de/sfb823-dp2010.html
More informationPredicting the Onset of Nonlinear Pharmacokinetics
Citation: CPT Pharmacometrics Syst Pharmacol (208) 7, 670 677; doi:0002/psp4236 ARTICLE Predicting the Onset of Nonlinear Pharmacokinetics Andrew M Stein and Lambertus A Peletier 2 When analyzing the pharmacokinetics
More informationBayesian Estimation of Input Output Tables for Russia
Bayesian Estimation of Input Output Tables for Russia Oleg Lugovoy (EDF, RANE) Andrey Polbin (RANE) Vladimir Potashnikov (RANE) WIOD Conference April 24, 2012 Groningen Outline Motivation Objectives Bayesian
More informationIntroduction to Crossover Trials
Introduction to Crossover Trials Stat 6500 Tutorial Project Isaac Blackhurst A crossover trial is a type of randomized control trial. It has advantages over other designed experiments because, under certain
More information1. Describing patient variability, 2. Minimizing patient variability, 3. Describing and Minimizing Intraindividual
1. Describing patient variability, 2. Minimizing patient variability, 3. Describing and Minimizing Intraindividual Variability. 4. Improving the reporting of lab assay results. Roger Jelliffe MD, Alan
More informationTest Strategies for Experiments with a Binary Response and Single Stress Factor Best Practice
Test Strategies for Experiments with a Binary Response and Single Stress Factor Best Practice Authored by: Sarah Burke, PhD Lenny Truett, PhD 15 June 2017 The goal of the STAT COE is to assist in developing
More informationUse of frequentist and Bayesian approaches for extrapolating from adult efficacy data to design and interpret confirmatory trials in children
Use of frequentist and Bayesian approaches for extrapolating from adult efficacy data to design and interpret confirmatory trials in children Lisa Hampson, Franz Koenig and Martin Posch Department of Mathematics
More informationLatent Variable Model for Weight Gain Prevention Data with Informative Intermittent Missingness
Journal of Modern Applied Statistical Methods Volume 15 Issue 2 Article 36 11-1-2016 Latent Variable Model for Weight Gain Prevention Data with Informative Intermittent Missingness Li Qin Yale University,
More informationA NOVEL OPTIMAL PROBABILITY DENSITY FUNCTION TRACKING FILTER DESIGN 1
A NOVEL OPTIMAL PROBABILITY DENSITY FUNCTION TRACKING FILTER DESIGN 1 Jinglin Zhou Hong Wang, Donghua Zhou Department of Automation, Tsinghua University, Beijing 100084, P. R. China Control Systems Centre,
More informationEstimating functional uncertainty using polynomial chaos and adjoint equations
0. Estimating functional uncertainty using polynomial chaos and adjoint equations February 24, 2011 1 Florida State University, Tallahassee, Florida, Usa 2 Moscow Institute of Physics and Technology, Moscow,
More informationInfluence of covariate distribution on the predictive performance of pharmacokinetic models in paediatric research
British Journal of Clinical Pharmacology DOI:1.1111/bcp.13 Influence of covariate distribution on the predictive performance of pharmacokinetic models in paediatric research Chiara Piana, 1 Meindert Danhof
More informationScore test for random changepoint in a mixed model
Score test for random changepoint in a mixed model Corentin Segalas and Hélène Jacqmin-Gadda INSERM U1219, Biostatistics team, Bordeaux GDR Statistiques et Santé October 6, 2017 Biostatistics 1 / 27 Introduction
More informationEstimation, Model Selection and Optimal Design in Mixed Eects Models Applications to pharmacometrics. Marc Lavielle 1. Cemracs CIRM logo
Estimation, Model Selection and Optimal Design in Mixed Eects Models Applications to pharmacometrics Marc Lavielle 1 1 INRIA Saclay Cemracs 2009 - CIRM Outline 1 Introduction 2 Some pharmacokinetics-pharmacodynamics
More informationML estimation: Random-intercepts logistic model. and z
ML estimation: Random-intercepts logistic model log p ij 1 p = x ijβ + υ i with υ i N(0, συ) 2 ij Standardizing the random effect, θ i = υ i /σ υ, yields log p ij 1 p = x ij β + σ υθ i with θ i N(0, 1)
More informationSolving the steady state diffusion equation with uncertainty Final Presentation
Solving the steady state diffusion equation with uncertainty Final Presentation Virginia Forstall vhfors@gmail.com Advisor: Howard Elman elman@cs.umd.edu Department of Computer Science May 6, 2012 Problem
More informationHow Critical is the Duration of the Sampling Scheme for the Determination of Half-Life, Characterization of Exposure and Assessment of Bioequivalence?
How Critical is the Duration of the Sampling Scheme for the Determination of Half-Life, Characterization of Exposure and Assessment of Bioequivalence? Philippe Colucci 1,2 ; Jacques Turgeon 1,3 ; and Murray
More information