A comparison of estimation methods in nonlinear mixed effects models using a blind analysis
|
|
- Imogene Fields
- 6 years ago
- Views:
Transcription
1 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, INSERM U738, University Paris VII, Paris France Project objectives Compare various recently developed and classical algorithms for non-linear mixed effect models using ML parametric estimate: Bias Precision Standard Error estimates PAGE 05 Pampeluna 2
2 The simulation model: one compartment, 1 st order absorption reparametrized to avoid flip-flop With CP Dose * Ka V *( Ka-Ke) () t = ( exp( - Ke t) - exp( - Ka t) )*exp( ε ) V = θ 1 EXP(η 1 ) KE = θ 2 EXP(η 2 ) KA = KE+ θ 3 EXP(η 3 ) θ 1, θ 2, θ 3 fixed effect parameters, greater than 0, for V, KA and KA-KE; η 1, η 2, η 3 are random effects representing between individual variability normally distributed, with a block VARIANCE matrix Ω : (η 1, η 2, η 3 ) ~ N ( 0,0,0), 2 ϖ 1 ϖ 21 ϖ 31 ϖ 31 ϖ 32 2 ϖ 3 ε random normally distributed effect representing within residual variability with ε ~N(0, σ 2 ). ϖ ϖ ϖ PAGE 05 Pampeluna 3 Simulation design simulated datasets individuals (ID) / dataset 4 samples / ID at times 1, 4, 12, 24 hrs + random shift: ε 1 ~ N(0, SD=0.25) at 1 and 4 hrs ε 2 ~ N(0, SD=2) at 12 and 24 hrs Missing points at random P(missing)=0.25 PAGE 05 Pampeluna 4
3 Parameters used for simulation V Ke Ka-Ke Value CV IIV 47% 81% 15% Corr IIV V Ke Ka-Ke V Ke Ka-Ke CV Resid = 25% PAGE 05 Pampeluna 5 Example of 4 simulated datasets datasim.txt Run = 1 datasim.txt Run = 2 Cp Cp Time datasim.txt Run = Time datasim.txt Run = 4 Cp Cp Time Time PAGE 05 Pampeluna 6
4 Blinded analysis Exact model was given to each participant who ignored the true simulation parameters Suggested initial parameter estimates provided only for fixed effects: V = 30 (true=27.2) Ke = 0.3 (true = 0.232) Ka-Ke =0.5 (true=0.304) Data sets blinded and send on 7th March 2005; Results received between March 25 th and June 10 th 05 PAGE 05 Pampeluna 7 Algorithms involved in the comparison Approximate Likelihood SAS proc nlmixed FO approximation NONMEM V and VI FOCE and (Splus 6.2) Lindstrom MJ, Bates DM. Nonlinear mixed effects models for repeated measures data. Biometrics 1990; 46: ITBS (MW\Pharm & MulTifit) Mentré F, Gomeni R. A two-step iterative algorithm for estimation in nonlinear mixed-effect models with an evaluation in population pharmacokinetics. J Biopharm Stat 1995; 5: Exact Likelihood 3 Stochastic methods PEM MCPEM SAEM One method based on approximated integral of exact likelihood SAS proc nlmixed Adaptive Gaussian Quadrature: 3 points per dimension (27 points in all) used to integrate the multidimensional exact likelihood PAGE 05 Pampeluna 8
5 Compilation of results: Relative error, bias and Root Mean Square Error For Ω, computation of 3 1 Det( Ω with true=1 est Ω true ) Relative Estimation error (RER) for each simulation presented as boxplots: RER= (Est True) / True Wilcoxon test for testing H0: RER=0 vs. H1: RER 0 Bias % = mean(est True) / True Precision: relative Root Mean Square Error (RMSE) : 2 ( Est True) ) 2 / % RMSE = mean True Results send out to participants on June 6 th PAGE 05 Pampeluna 9 RESULTS
6 All algorithms converged and provided estimates for all simulated data sets except FOCE (NONMEM and ) 1 st try 1 2 nd try NONMEM Version V 2 NONMEM Version VI 2 (Beta) Only Est. step successful % 88% % 49% Est and COV successful % 70% 47% 49% 1. Biased estimates with 1 st set of results 2 nd set of results were 2. repeated search with perturbed initial estimates Almost perfect matching on failed run# 7, 9, 28, 45, 47, 53, 68, 69, 76, 80, 90 PAGE 05 Pampeluna 11 V true= Pval Wilcox= N= NM V NM VI MWPharm MultiFit SAS FO SAS GQ MCPEM PEM SAEM Ke true= Pval Wilcox= N= NM V NM VI MWPharm MultiFit SAS FO SAS GQ MCPEM PEM SAEM Ka-Ke true= NM V NM VI MWPharm MultiFit SAS FO SAS GQ MCPEM PEM SAEM Pval Wilcox= N= PAGE 05 Pampeluna 12
7 om11 V true= Pval Wilcox= N= NM V NM VI MWPharm MultiFit SAS FO SAS GQ MCPEM PEM SAEM om22 Ke true= Pval Wilcox= N= NM V NM VI MWPharm MultiFit SAS FO SAS GQ MCPEM PEM SAEM om33 Ka-Ke true= NM V NM VI MWPharm MultiFit SAS FO SAS GQ MCPEM PEM SAEM Pval Wilcox= N= PAGE 05 Pampeluna 13 Cor12 V.Ke true= Pval Wilcox= N= NM V NM VI MWPharm MultiFit SAS FO SAS GQ MCPEM PEM SAEM Cor13 V.KaKe true= Pval Wilcox= N= NM V NM VI MWPharm MultiFit SAS FO SAS GQ MCPEM PEM SAEM Cor23 Ke.KaKe true= NM V NM VI MWPharm MultiFit SAS FO SAS GQ MCPEM PEM SAEM Pval Wilcox= N= PAGE 05 Pampeluna 14
8 Det(OM.est*inv(OM.true))**(1/3) true= Pval Wilcox= N= NM V NM VI MWPharm MultiFit SAS FO SAS GQ MCPEM PEM SAEM S**2 true= Pval Wilcox= 0.01 N= NM V NM VI MWPharm MultiFit SAS FO SAS GQ MCPEM PEM SAEM PAGE 05 Pampeluna 15 Bias % for fixed effects +<5% +5-10% % +>20% TH1 (V=27.2) TH2 (Ke=0.23) TH3 (Ka-Ke=0.30) NM V NM VI * ** MWPharm MultiFit SAS FO SAS GQ MCPEM PEM SAEM * NM VI based on 49 successful runs out of ** 2 nd chance and 88 successful runs out of PAGE 05 Pampeluna 16
9 Bias % for variance omega**2 +<20% % +50-% +>% om V (tr=0.22) om Ke (tr=0.65) om33 (tr=0.02) NM V NM VI * ** MWPharm MultiFit SAS FO SAS GQ MCPEM PEM SAEM * NM VI based on 49 successful runs out of ** 2 nd chance and 88 successful runs out of PAGE 05 Pampeluna 17 Bias % for correlation coefficients +<20% +20-% +-500% +>500% V-Ke (tr=0.91) V-KaKe (tr=0.33) Ke-KaKe (tr=0.01) NM V NM VI * ** MWPharm MultiFit SAS FO SAS GQ MCPEM PEM SAEM * NM VI based on 49 successful runs out of ** 2 nd chance and 88 successful runs out of PAGE 05 Pampeluna 18
10 Bias % for det( Ω.inv(Ω tr ) ) and S 2 +<10% % % +>50% det( OM.inv(OMtr) ) (true=1) S**2 (true=0.06) NM V NM VI * ** MWPharm MultiFit SAS FO SAS GQ MCPEM PEM SAEM * NM VI based on 49 successful runs out of ** 2 nd chance and 88 successful runs out of PAGE 05 Pampeluna 19 RMSE Results
11 RMSE % for fixed effects <10% 10-20% 20-50% >50% TH1 (V=27.2) TH2 (Ke=0.23) TH3 (Ka-Ke=0.30) NM V NM VI * ** MWPharm MultiFit SAS FO SAS GQ MCPEM PEM SAEM * NM VI based on 49 successful runs out of ** 2 nd chance and 88 successful runs out of PAGE 05 Pampeluna 21 RMSE % for variance omega**2 <20% 20-50% 50-% >% om V (tr=0.22) om Ke (tr=0.65) om33 (tr=0.02) NM V NM VI * ** MWPharm MultiFit SAS FO SAS GQ MCPEM PEM SAEM * NM VI based on 49 successful runs out of ** 2 nd chance and 88 successful runs out of PAGE 05 Pampeluna 22
12 RMSE % for correlation coefficients <20% 20-% -500% >500% V-Ke (tr=0.91) V-KaKe (tr=0.33) Ke-KaKe (tr=0.01) NM V NM VI * ** MWPharm MultiFit SAS FO SAS GQ MCPEM PEM SAEM * NM VI based on 49 successful runs out of ** 2 nd chance and 88 successful runs out of PAGE 05 Pampeluna 23 RMSE % for det( Ω.inv(Ω tr ) ) and S 2 <20% 20-50% 50-% >% det( OM.inv(OMtr) ) (true=1) S**2 (true=0.06) NM V NM VI * MWPharm MultiFit SAS FO SAS GQ MCPEM PEM SAEM * NM VI based on 49 successful runs out of ** 2 nd chance and 88 successful runs out of PAGE 05 Pampeluna 24
13 Comparison between 1 st set and 2 nd one for : «Because we were surprised with the blinded results regarding NLME, we went back to the code we used to produce the fits and realized that we had specified the model in a unstable way that created convergence problems in NLME. We re-defined the model formulation in S+ and re-ran the analyses. The fits were considerably more stable than before.» Hsu & Pinheiro V Ke Ka-Ke st set 2nd set 1st set 2nd set 1st set 2nd set om11 V om22 Ke om33 Ka-Ke st set 2nd set st set 2nd set st set 2nd set PAGE 05 Pampeluna 25 Use of non-parametric ranks to classify all methods Rank TH1 (V=27.2) NM V NM VI * ** MWPharm MultiFit SAS FO SAS GQ MCPEM PEM 0.9 SAEM Quote = rank k= 1 10 Quote k * NM VI is beta version: only 49% convergence ** : only 88% convergence with 2 nd try PAGE 05 Pampeluna 26
14 Classification of methods based on their sum of ranks for bias and RMSE SAEM PEM BIAS SAS GQ PEM 84 SAS GQ SAEM PEM MCPEM NM V NM VI* SAEM RMSE ** MWPharm MultiFit MCPEM SAS FO SAEM 89 MCPEM 78 PEM 73 SAS GQ 70 NM V 64 NM VI 59 MWPharm 58 MultiFit SAS FO 17 PAGE 05 Pampeluna 27 CPU Times Adj_time=time_one_run GHz SAEM 1 PEM MCPEM SAS GQ SAS FO MultiFit MWPhar m Nlme NM VI NM V 3 GHz CPU time for 1 run (sec) 6 Adjusted CPU times SAEM, 32 PEM, 612 MCPEM, 360 SAS GQ, 19 SAS FO, 2 MultiFit, 13 MWPharm, 40, 11 NM VI, 6 NM V, Time (sec) 1. A compiled version would probably be considerably faster than actual implementation in Matlab for PEM and SAEM 2. For a successful nls + run 3. For a successful NONMEM run (EST & COV) PAGE 05 Pampeluna 28
15 Conclusions Methods without approximation of the likelihood (either stochastic methods or GQ) show better performances than others with an approximation. Reduction of BIAS and RMSE with stochastic methods is obtained only with slight CPU time increase for SAEM FO in SAS shows, as expected, large bias FOCE methods, either implemented in or in NONMEM have either convergence problems or computational ones for covariance matrix with these data sets Future NONMEM VI does not solve NM V issues but rather seems to be even more stringent than version V Datasets will be made available soon on a web site there are a number of buts and ifs when you make comparisons of this sort (Niclas) PAGE 05 Pampeluna 29 Acknowledgements for accepting to participate into this blind comparison (in alphabetical order) Serge Guzy (MCPEM, Xoma corp., USA) Niclas Jonsson (NONMEM FOCE V & VI. Univ. Uppsala, Sweden) Marc Lavielle (SAEM, Univ. Paris Sud Orsay, FR) Bob Leary (PEM; Pharsight, USA) José Pinheiro and Chyi-Hung Hsu (Splus. Novartis, USA) Hans Proost (MW\Pharm and MultiFit. Univ. Groningen, NL) Russ Wolfinger (SAS PROC NLMIXED, USA) PAGE 05 Pampeluna 30
Stochastic 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 informationPharmacokinetics Introduction to
Wikimedia Commons Korinna Creative Commons Attribution-ShareAlike 3. Unported Introduction to Population PK Pharmacokinetics Introduction to Introduction to Population PK Helmut Schütz BEBAC Workshop Bucarest,
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 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 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 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 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 informationEvaluation 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 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 informationPopulation Approach. Pharmacokinetics. Analysis Approaches. Describing The Signal and the Noise
1 Population Approach Describing The Signal and the Noise 2 Pharmacometrics Science Pharmacology Models Pharmacokinetics Pharmacodynamics Experiment Dose Concentration Effect 3 Analysis Approaches ANOVA
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 informationApproaches to Modeling Menstrual Cycle Function
Approaches to Modeling Menstrual Cycle Function Paul S. Albert (albertp@mail.nih.gov) Biostatistics & Bioinformatics Branch Division of Epidemiology, Statistics, and Prevention Research NICHD SPER Student
More informationON THE CONSEQUENCES OF MISSPECIFING ASSUMPTIONS CONCERNING RESIDUALS DISTRIBUTION IN A REPEATED MEASURES AND NONLINEAR MIXED MODELLING CONTEXT
ON THE CONSEQUENCES OF MISSPECIFING ASSUMPTIONS CONCERNING RESIDUALS DISTRIBUTION IN A REPEATED MEASURES AND NONLINEAR MIXED MODELLING CONTEXT Rachid el Halimi and Jordi Ocaña Departament d Estadística
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 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 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 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 informationMETHODS FOR POPULATION PHARMACOKINETICS AND PHARMACODYNAMICS
METHODS FOR POPULATION PHARMACOKINETICS AND PHARMACODYNAMICS Emily Anne Colby A dissertation submitted to the faculty of the University of North Carolina at Chapel Hill in partial fulfillment of the requirements
More informationRobust design in model-based analysis of longitudinal clinical data
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 Context Optimal
More informationParameter Estimation in Nonlinear Mixed Effect Models Using saemix, an R Implementation of the SAEM Algorithm
Parameter Estimation in Nonlinear Mixed Effect Models Using saemix, an R Implementation of the SAEM Algorithm Emmanuelle Comets, Audrey Lavenu, Marc Lavielle To cite this version: Emmanuelle Comets, Audrey
More informationData Abstraction Form for population PK, PD publications
Data Abstracton Form for populaton PK/PD publcatons Brendel K. 1*, Dartos C. 2*, Comets E. 1, Lemenuel-Dot A. 3, Laffont C.M. 3, Lavelle C. 4, Grard P. 2, Mentré F. 1 1 INSERM U738, Pars, France 2 EA3738,
More informationMS&E 226: Small Data. Lecture 11: Maximum likelihood (v2) Ramesh Johari
MS&E 226: Small Data Lecture 11: Maximum likelihood (v2) Ramesh Johari ramesh.johari@stanford.edu 1 / 18 The likelihood function 2 / 18 Estimating the parameter This lecture develops the methodology behind
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 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 informationAlgebra of Random Variables: Optimal Average and Optimal Scaling Minimising
Review: Optimal Average/Scaling is equivalent to Minimise χ Two 1-parameter models: Estimating < > : Scaling a pattern: Two equivalent methods: Algebra of Random Variables: Optimal Average and Optimal
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 informationCO2 Handout. t(cbind(co2$type,co2$treatment,model.matrix(~type*treatment,data=co2)))
CO2 Handout CO2.R: library(nlme) CO2[1:5,] plot(co2,outer=~treatment*type,layout=c(4,1)) m1co2.lis
More informationFitting PK Models with SAS NLMIXED Procedure Halimu Haridona, PPD Inc., Beijing
PharmaSUG China 1 st Conference, 2012 Fitting PK Models with SAS NLMIXED Procedure Halimu Haridona, PPD Inc., Beijing ABSTRACT Pharmacokinetic (PK) models are important for new drug development. Statistical
More informationKriging models with Gaussian processes - covariance function estimation and impact of spatial sampling
Kriging models with Gaussian processes - covariance function estimation and impact of spatial sampling François Bachoc former PhD advisor: Josselin Garnier former CEA advisor: Jean-Marc Martinez Department
More informationStat 5100 Handout #26: Variations on OLS Linear Regression (Ch. 11, 13)
Stat 5100 Handout #26: Variations on OLS Linear Regression (Ch. 11, 13) 1. Weighted Least Squares (textbook 11.1) Recall regression model Y = β 0 + β 1 X 1 +... + β p 1 X p 1 + ε in matrix form: (Ch. 5,
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 informationSAS Code: Joint Models for Continuous and Discrete Longitudinal Data
CHAPTER 14 SAS Code: Joint Models for Continuous and Discrete Longitudinal Data We show how models of a mixed type can be analyzed using standard statistical software. We mainly focus on the SAS procedures
More informationProduct Held at Accelerated Stability Conditions. José G. Ramírez, PhD Amgen Global Quality Engineering 6/6/2013
Modeling Sub-Visible Particle Data Product Held at Accelerated Stability Conditions José G. Ramírez, PhD Amgen Global Quality Engineering 6/6/2013 Outline Sub-Visible Particle (SbVP) Poisson Negative Binomial
More informationMaximum Likelihood Estimation. only training data is available to design a classifier
Introduction to Pattern Recognition [ Part 5 ] Mahdi Vasighi Introduction Bayesian Decision Theory shows that we could design an optimal classifier if we knew: P( i ) : priors p(x i ) : class-conditional
More informationPerformance Comparison of K-Means and Expectation Maximization with Gaussian Mixture Models for Clustering EE6540 Final Project
Performance Comparison of K-Means and Expectation Maximization with Gaussian Mixture Models for Clustering EE6540 Final Project Devin Cornell & Sushruth Sastry May 2015 1 Abstract In this article, we explore
More informationLecture Stat Information Criterion
Lecture Stat 461-561 Information Criterion Arnaud Doucet February 2008 Arnaud Doucet () February 2008 1 / 34 Review of Maximum Likelihood Approach We have data X i i.i.d. g (x). We model the distribution
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 informationChapter 3: Maximum-Likelihood & Bayesian Parameter Estimation (part 1)
HW 1 due today Parameter Estimation Biometrics CSE 190 Lecture 7 Today s lecture was on the blackboard. These slides are an alternative presentation of the material. CSE190, Winter10 CSE190, Winter10 Chapter
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 informationComputational statistics
Computational statistics Lecture 3: Neural networks Thierry Denœux 5 March, 2016 Neural networks A class of learning methods that was developed separately in different fields statistics and artificial
More informationDS-GA 1002 Lecture notes 11 Fall Bayesian statistics
DS-GA 100 Lecture notes 11 Fall 016 Bayesian statistics In the frequentist paradigm we model the data as realizations from a distribution that depends on deterministic parameters. In contrast, in Bayesian
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 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 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 informationCMU-Q Lecture 24:
CMU-Q 15-381 Lecture 24: Supervised Learning 2 Teacher: Gianni A. Di Caro SUPERVISED LEARNING Hypotheses space Hypothesis function Labeled Given Errors Performance criteria Given a collection of input
More informationW-BASED VS LATENT VARIABLES SPATIAL AUTOREGRESSIVE MODELS: EVIDENCE FROM MONTE CARLO SIMULATIONS
1 W-BASED VS LATENT VARIABLES SPATIAL AUTOREGRESSIVE MODELS: EVIDENCE FROM MONTE CARLO SIMULATIONS An Liu University of Groningen Henk Folmer University of Groningen Wageningen University Han Oud Radboud
More informationModel structure. Lecture Note #3 (Chap.6) Identification of time series model. ARMAX Models and Difference Equations
System Modeling and Identification Lecture ote #3 (Chap.6) CHBE 70 Korea University Prof. Dae Ryoo Yang Model structure ime series Multivariable time series x [ ] x x xm Multidimensional time series (temporal+spatial)
More informationMustafa H. Tongarlak Bruce E. Ankenman Barry L. Nelson
Proceedings of the 0 Winter Simulation Conference S. Jain, R. R. Creasey, J. Himmelspach, K. P. White, and M. Fu, eds. RELATIVE ERROR STOCHASTIC KRIGING Mustafa H. Tongarlak Bruce E. Ankenman Barry L.
More informationModel Assumptions; Predicting Heterogeneity of Variance
Model Assumptions; Predicting Heterogeneity of Variance Today s topics: Model assumptions Normality Constant variance Predicting heterogeneity of variance CLP 945: Lecture 6 1 Checking for Violations of
More informationdans les modèles à vraisemblance non explicite par des algorithmes gradient-proximaux perturbés
Inférence pénalisée dans les modèles à vraisemblance non explicite par des algorithmes gradient-proximaux perturbés Gersende Fort Institut de Mathématiques de Toulouse, CNRS and Univ. Paul Sabatier Toulouse,
More informationA Novel Screening Method Using Score Test for Efficient Covariate Selection in Population Pharmacokinetic Analysis
A Novel Screening Method Using Score Test for Efficient Covariate Selection in Population Pharmacokinetic Analysis Yixuan Zou 1, Chee M. Ng 1 1 College of Pharmacy, University of Kentucky, Lexington, KY
More informationABC random forest for parameter estimation. Jean-Michel Marin
ABC random forest for parameter estimation Jean-Michel Marin Université de Montpellier Institut Montpelliérain Alexander Grothendieck (IMAG) Institut de Biologie Computationnelle (IBC) Labex Numev! joint
More informationAGEC 661 Note Fourteen
AGEC 661 Note Fourteen Ximing Wu 1 Selection bias 1.1 Heckman s two-step model Consider the model in Heckman (1979) Y i = X iβ + ε i, D i = I {Z iγ + η i > 0}. For a random sample from the population,
More informationTemporal context calibrates interval timing
Temporal context calibrates interval timing, Mehrdad Jazayeri & Michael N. Shadlen Helen Hay Whitney Foundation HHMI, NPRC, Department of Physiology and Biophysics, University of Washington, Seattle, Washington
More informationSAS PROC NLMIXED Mike Patefield The University of Reading 12 May
SAS PROC NLMIXED Mike Patefield The University of Reading 1 May 004 E-mail: w.m.patefield@reading.ac.uk non-linear mixed models maximum likelihood repeated measurements on each subject (i) response vector
More informationEfficient Estimation for the Partially Linear Models with Random Effects
A^VÇÚO 1 33 ò 1 5 Ï 2017 c 10 Chinese Journal of Applied Probability and Statistics Oct., 2017, Vol. 33, No. 5, pp. 529-537 doi: 10.3969/j.issn.1001-4268.2017.05.009 Efficient Estimation for the Partially
More informationParametric Models. Dr. Shuang LIANG. School of Software Engineering TongJi University Fall, 2012
Parametric Models Dr. Shuang LIANG School of Software Engineering TongJi University Fall, 2012 Today s Topics Maximum Likelihood Estimation Bayesian Density Estimation Today s Topics Maximum Likelihood
More informationA New Approach to Modeling Covariate Effects and Individualization in Population Pharmacokinetics-Pharmacodynamics
Journal of Pharmacokinetics and Pharmacodynamics, Vol. 33, No. 1, February 2006 ( 2006) DOI: 10.1007/s10928-005-9000-2 A New Approach to Modeling Covariate Effects and Individualization in Population Pharmacokinetics-Pharmacodynamics
More informationPrediction of double gene knockout measurements
Prediction of double gene knockout measurements Sofia Kyriazopoulou-Panagiotopoulou sofiakp@stanford.edu December 12, 2008 Abstract One way to get an insight into the potential interaction between a pair
More informationAssessment of uncertainty in computer experiments: from Universal Kriging to Bayesian Kriging. Céline Helbert, Delphine Dupuy and Laurent Carraro
Assessment of uncertainty in computer experiments: from Universal Kriging to Bayesian Kriging., Delphine Dupuy and Laurent Carraro Historical context First introduced in the field of geostatistics (Matheron,
More informationNon-linear mixed-effects pharmacokinetic/pharmacodynamic modelling in NLME using differential equations
Computer Methods and Programs in Biomedicine () 76, 31 Non-linear mixed-effects pharmacokinetic/pharmacodynamic modelling in NLME using differential equations Christoffer W.Tornøe a, *, Henrik Agersø a,
More informationPattern Recognition. Parameter Estimation of Probability Density Functions
Pattern Recognition Parameter Estimation of Probability Density Functions Classification Problem (Review) The classification problem is to assign an arbitrary feature vector x F to one of c classes. The
More informationA discrete semi-markov model for the effect of need-based treatments on the disease states Ekkehard Glimm Novartis Pharma AG
A discrete semi-markov model for the effect of need-based treatments on the disease states Ekkehard Glimm Novartis Pharma AG Basel, Switzerland Lillian Yau Sandoz Biopharmaceuticals Munich, Germany Agenda
More informationRestricted Maximum Likelihood in Linear Regression and Linear Mixed-Effects Model
Restricted Maximum Likelihood in Linear Regression and Linear Mixed-Effects Model Xiuming Zhang zhangxiuming@u.nus.edu A*STAR-NUS Clinical Imaging Research Center October, 015 Summary This report derives
More informationOutline Lecture 2 2(32)
Outline Lecture (3), Lecture Linear Regression and Classification it is our firm belief that an understanding of linear models is essential for understanding nonlinear ones Thomas Schön Division of Automatic
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 informationPh.D. Preliminary Examination Statistics June 2, 2014
Ph.D. Preliminary Examination Statistics June, 04 NOTES:. The exam is worth 00 points.. Partial credit may be given for partial answers if possible.. There are 5 pages in this exam paper. I have neither
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 informationOn Identification of Cascade Systems 1
On Identification of Cascade Systems 1 Bo Wahlberg Håkan Hjalmarsson Jonas Mårtensson Automatic Control and ACCESS, School of Electrical Engineering, KTH, SE-100 44 Stockholm, Sweden. (bo.wahlberg@ee.kth.se
More informationWhen is a nonlinear mixed-effects model identifiable?
When is a nonlinear mixed-effects model identifiable? O. G. Nunez and D. Concordet Universidad Carlos III de Madrid and Ecole Nationale Vétérinaire de Toulouse Abstract We consider the identifiability
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 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 informationMaster s Written Examination - Solution
Master s Written Examination - Solution Spring 204 Problem Stat 40 Suppose X and X 2 have the joint pdf f X,X 2 (x, x 2 ) = 2e (x +x 2 ), 0 < x < x 2
More information36-720: Linear Mixed Models
36-720: Linear Mixed Models Brian Junker October 8, 2007 Review: Linear Mixed Models (LMM s) Bayesian Analogues Facilities in R Computational Notes Predictors and Residuals Examples [Related to Christensen
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 informationA. Motivation To motivate the analysis of variance framework, we consider the following example.
9.07 ntroduction to Statistics for Brain and Cognitive Sciences Emery N. Brown Lecture 14: Analysis of Variance. Objectives Understand analysis of variance as a special case of the linear model. Understand
More informationResearch Article. Jacob Leander, 1,2 Joachim Almquist, 1,3 Christine Ahlström, 4 Johan Gabrielsson, 5 and Mats Jirstrand 1,6
The AAPS Journal ( # 2015) DOI: 10.1208/s12248-015-9718-8 Research Article Mixed Effects Modeling Using Stochastic Differential Equations: Illustrated by Pharmacokinetic Data of Nicotinic Acid in Obese
More informationGeneral Regression Model
Scott S. Emerson, M.D., Ph.D. Department of Biostatistics, University of Washington, Seattle, WA 98195, USA January 5, 2015 Abstract Regression analysis can be viewed as an extension of two sample statistical
More informationMaximum Likelihood Conjoint Measurement in R
Maximum Likelihood Conjoint Measurement in R Kenneth Knoblauch¹, Blaise Tandeau de Marsac¹, Laurence T. Maloney² 1. Inserm, U846 Stem Cell and Brain Research Institute Dept. Integrative Neurosciences Bron,
More informationNemours Biomedical Research Biostatistics Core Statistics Course Session 4. Li Xie March 4, 2015
Nemours Biomedical Research Biostatistics Core Statistics Course Session 4 Li Xie March 4, 2015 Outline Recap: Pairwise analysis with example of twosample unpaired t-test Today: More on t-tests; Introduction
More informationECE521 lecture 4: 19 January Optimization, MLE, regularization
ECE521 lecture 4: 19 January 2017 Optimization, MLE, regularization First four lectures Lectures 1 and 2: Intro to ML Probability review Types of loss functions and algorithms Lecture 3: KNN Convexity
More informationComparing MLE, MUE and Firth Estimates for Logistic Regression
Comparing MLE, MUE and Firth Estimates for Logistic Regression Nitin R Patel, Chairman & Co-founder, Cytel Inc. Research Affiliate, MIT nitin@cytel.com Acknowledgements This presentation is based on joint
More informationStatistics 262: Intermediate Biostatistics Non-parametric Survival Analysis
Statistics 262: Intermediate Biostatistics Non-parametric Survival Analysis Jonathan Taylor & Kristin Cobb Statistics 262: Intermediate Biostatistics p.1/?? Overview of today s class Kaplan-Meier Curve
More informationResults of a simulation of modeling and nonparametric methodology for count data (from patient diary) in drug studies
Results of a simulation of modeling and nonparametric methodology for count data (from patient diary) in drug studies Erhard Quebe-Fehling Workshop Stat. Methoden für korrelierte Daten Bochum 24-Nov-06
More informationSAS Example 3: Deliberately create numerical problems
SAS Example 3: Deliberately create numerical problems Four experiments 1. Try to fit this model, failing the parameter count rule. 2. Set φ 12 =0 to pass the parameter count rule, but still not identifiable.
More informationNonrecursive Models Highlights Richard Williams, University of Notre Dame, https://www3.nd.edu/~rwilliam/ Last revised April 6, 2015
Nonrecursive Models Highlights Richard Williams, University of Notre Dame, https://www3.nd.edu/~rwilliam/ Last revised April 6, 2015 This lecture borrows heavily from Duncan s Introduction to Structural
More informationEstimation in Generalized Linear Models with Heterogeneous Random Effects. Woncheol Jang Johan Lim. May 19, 2004
Estimation in Generalized Linear Models with Heterogeneous Random Effects Woncheol Jang Johan Lim May 19, 2004 Abstract The penalized quasi-likelihood (PQL) approach is the most common estimation procedure
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 informationPublished online: 10 Apr 2012.
This article was downloaded by: Columbia University] On: 23 March 215, At: 12:7 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 172954 Registered office: Mortimer
More informationOn Fitting Generalized Linear Mixed Effects Models for Longitudinal Binary Data Using Different Correlation
On Fitting Generalized Linear Mixed Effects Models for Longitudinal Binary Data Using Different Correlation Structures Authors: M. Salomé Cabral CEAUL and Departamento de Estatística e Investigação Operacional,
More informationThe First Thing You Ever Do When Receive a Set of Data Is
The First Thing You Ever Do When Receive a Set of Data Is Understand the goal of the study What are the objectives of the study? What would the person like to see from the data? Understand the methodology
More informationCovariance function estimation in Gaussian process regression
Covariance function estimation in Gaussian process regression François Bachoc Department of Statistics and Operations Research, University of Vienna WU Research Seminar - May 2015 François Bachoc Gaussian
More informationA Nonlinear Random Coefficients Model. for Degradation Testing
A Nonlinear Random Coefficients Model for Degradation Testing Suk Joo Bae Paul H. Kvam Georgia Institute of Technology February 8, 2004 Abstract As an alternative to traditional life testing, degradation
More informationCovariate Model Building in Nonlinear Mixed Effects Models
Digital Comprehensive Summaries of Uppsala Dissertations from the Faculty of Pharmacy 59 Covariate Model uilding in Nonlinear Mixed Effects Models JKO RIING CT UNIVERSITTIS UPSLIENSIS UPPSL 2007 ISSN 1651-6192
More informationAlgebra of Random Variables: Optimal Average and Optimal Scaling Minimising
Review: Optimal Average/Scaling is equivalent to Minimise χ Two 1-parameter models: Estimating < > : Scaling a pattern: Two equivalent methods: Algebra of Random Variables: Optimal Average and Optimal
More informationRandom Features for Large Scale Kernel Machines
Random Features for Large Scale Kernel Machines Andrea Vedaldi From: A. Rahimi and B. Recht. Random features for large-scale kernel machines. In Proc. NIPS, 2007. Explicit feature maps Fast algorithm for
More informationMeasuring relationships among multiple responses
Measuring relationships among multiple responses Linear association (correlation, relatedness, shared information) between pair-wise responses is an important property used in almost all multivariate analyses.
More informationDIAGNOSTICS FOR STRATIFIED CLINICAL TRIALS IN PROPORTIONAL ODDS MODELS
DIAGNOSTICS FOR STRATIFIED CLINICAL TRIALS IN PROPORTIONAL ODDS MODELS Ivy Liu and Dong Q. Wang School of Mathematics, Statistics and Computer Science Victoria University of Wellington New Zealand Corresponding
More information