Heterogeneous shedding of influenza by human subjects and. its implications for epidemiology and control
|
|
- Randell Ellis
- 5 years ago
- Views:
Transcription
1 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 3,4,5 1 Centre for Immunity, Infection & Evolution, University of Edinburgh, UK 2 Influenza Division, National Center for Immunization and Respiratory Disease, Centers for Disease Control and Prevention, Atlanta, Georgia, USA 3 INSERM, UMR_S 1136, Institut Pierre Louis d Epidémiologie et de Santé Publique, Paris, France 4 Sorbonne Universités, UPMC Univ Paris 06, UMR_S 1136, Institut Pierre Louis d Epidémiologie et de Santé Publique, Paris, France 5 Assistance Publique Hôpitaux de Paris, Hôpital Saint Antoine, Paris, France Supplementary material Text S1: Supplementary methods Description of the nonlinear mixed effect models. Nonlinear mixed effect models (or the population approach) allows a description of population characteristics (fixed effect) as well as the inter-individual variability (IIV, random effect) 1. In this method, a function f describing the variables being modeled, e.g., the viral load, depends nonlinearly on θ i, a vector of the p parameters of subject i. A vector ξ i representing the times at which samples are collected from subject i, ξ i = (t i1 ; t i2 ; ; t in ), is also considered. The statistical model for subject i is then given by: 23 y i f ; a bf ; i i i i i 1
2 where y i is a vector with n i observations of subject i, with i varying from 1 to N, i is the vector of the residual errors which is the part of the observations unexplained by the model f. It is assumed that the errors i are independent from one observation to another and that their distribution is Gaussian i ~ N(0; I ni ), where I ni is an identity matrix of dimension n i. a and b are two parameters characterizing the error model variance. To fit the data, we used 6 an additive error model: y f a, where y is the observed response, f is the model function, a is the additive error term, and the error,, is normally distributed following N(0,1). In nonlinear mixed effect models, the model f is common to all the subjects, but the vector of parameters θ i for subject i may vary from one subject to another. The interindividual variability is modeled with the vector of random effect parameters i. The vector of parameters θ i for the subject i can then be expressed as a second-level model which links with the function g, the vector of fixed effect parameters common for all subject and the vector of random effects b i specific for subject i: θ i = g( ; i ). The vector of random effect is assumed to follow a Gaussian distribution i ~ N(0; Ω), i and i are assumed to be independent for subject i and i i is assumed independent from one subject to another. Ω is the matrix of random effect variance and covariances. The function g is here an exponential model. The vector of parameters is hence written as θ = βe η i. Parameter estimation and statistical methods. Population parameter estimates and inter-individual variability (IIV) estimates were obtained using a maximum-likelihood method implemented in MONOLIX version ( which uses 22 the stochastic approximation expectation-approximation (SAEM) algorithm 2 to estimate population parameters. The VK model was fitted to log 10 viral load. HCV genotype, cohort and patient type (treatment naïve vs. non-responder) were included as covariates in the model 2
3 to study their effect on the PK/VK parameters. Individual parameters were estimated using the empirical Bayes method 3. The model was simultaneously fitted to VK, systemic and respiratory symptom data since simultaneous fitting has been shown to provide more precise estimates 4,5. Parameters estimated were b, p,,,, s, and. For each parameter, we report the population estimates and their relative standard errors, as well as the inter-individual variability in percent with its relative standard error (Table 1). Individual parameters were estimated using the empirical Bayes method 3 and are presented in Tables S1. In the following we present only the significant results Dimension analysis The ODEs compartments are the following: T: unit = cell I: unit = cell V: unit=tcid50/ml F L and F S : unit= arbitrarily set to F (this represent the cytokines levels) N: cytotoxic activity: unit= Cell -1.d -1 Cell.d -1 =beta.cell.tcid50/ml Cell.d -1 =beta.cell.tcid50/ml delta.cell eta.cell. Cell -1.d -1 TCID50/mL.d -1 =p.cell/(1+psi.f) s.tcid50/ml - sigma.beta.cell.tcid50/ml F.d -1 =gamma.cell alpha L.F F.d -1 = rho.f - alpha S.F Cell -1.d -2 = k.f We can then deduce that: beta = (TCID50/mL) -1 d -1 delta= d -1 eta = Cell psi = F -1 c = d -1 sigma = Cell -1.TCID 50 /ml p = TCID50/mL.d -1.Cell -1 gamma = F.Cell -1.d -1 alpha L = d -1 rho=d -1 alpha S =d -1 3
4 k= F -1.Cell -1.d -2 Practical identifiability results To choose the values for these parameters, we re-estimated the parameters for every combination of the following V 0 =0.1, 0.5 and 1.0 TCID 50 /ml, =0.5, 1.0 and 2.0 d -1 and c=1.0, 3.0 and 10.0 d -1. These values are in the same range as the parameter estimates obtained in previous human vivo studies 6-8. We chose the set that provided the smallest Bayesian information criteria (BIC), where V 0 =0.5 TCID 50 /ml, =0.5 d -1 and c=1.0 d -1 (Fig. S1). We found that four parameters, p,, S, S did not vary significantly between the different fits and that infectivity rate,, tended to be lower whereas production of local cytokines ( L ) and diffusion of cytokines to circulation ( ) tended to be higher with lower virus clearance, c. The effect of cytokines ( ) tended to be higher with lower infected cell mortality rate ( ) (Fig. S1)
5 Supplementary Information references 1 Karlsson, M. & Sheiner, L. The importance of modeling interoccasion variability in population pharmacokinetic analyses. J Pharmacokinet Pharmacodyn 21, (1993). 2 Kuhn, E. & Lavielle, M. Maximum likelihood estimation in nonlinear mixed effects models. Comput Stat Data Anal 49, (2005). 3 Pinheiro, J. & Bates, D. (New York: Springer Verlag, 2000). 4 Zhang, L., Beal, S. L. & Sheiner, L. B. Simultaneous vs. sequential analysis for population PK/PD data I: Best-case performance. J Pharmacokinet Pharmacodyn 30, , doi: /b:jopa f (2003). 5 Savic, R. & Karlsson, M. Importance of shrinkage in empirical Bayes estimates for diagnostics: problems and solutions. AAPS J 11, (2009). 6 Baccam, P., Beauchemin, C., Macken, C. A., Hayden, F. G. & Perelson, A. S. Kinetics of influenza A virus infection in humans. Journal of virology 80, , doi: /jvi (2006). 7 Canini, L. & Carrat, F. Population modeling of influenza A/H1N1 virus kinetics and symptom dynamics. Journal of virology 85, (2011). 8 Hancioglu, B., Swigon, D. & Clermont, G. A dynamical model of human immune response to influenza A virus infection. Journal of theoretical biology 246, 70-86, doi: /j.jtbi (2007). 5
6 ID Infectivity rate (TCID50/mL) -1 d -1 x10-8 Virus production rate P (TCID50/mL) -1 d - 1 Cell -1 Table S1: Parameter estimates (mode) for individual patients. Individual parameters are Empirical Bayes Estimates Local cytokines effect F -1 Cytotoxic activity [-] Local cytokines clearance rate L Fd -1 ] Systemic cytokines production rate [F -1 ] Systemic cytokines production rate s [d -1 ] Local cytokines production rate [F.Cell -1.d -1 ]
7
8 Table S2: Estimated correlation matrix of the random effects. The relative standard error (rse) is the ratio of the standard error divided by the average value and is computed for each population estimate and IIV. p p L S (80) (32) (31) (21) (349) (18) L (33) (33) (5) (32) (208) (28) (128) (216) (94) S (121) (37) (31) (117) (22) (35) (9) (414) (20) (29) (14) (105) (37) 8
9 Figure S1. Boxplots of BIC and parameter estimates obtained after fitting the model for different values of (orange), V 0 (blue) and c (green). We re-estimated the model parameters for different fixed values of, V 0 and c. We tested our model for =1.0, 3.0 and 10.0 d -1, for V 0 =0.1, 0.5 and 1 log 10 TCID 50 /ml and for c=0.1, 0.5 and 1.0 d -1. We represent here the BIC as well as the population parameters distribution depending on the value of the fixed parameters. 9
10 A B 10
11 C Figure S2. Individual fits (blue line) to observation (black dots) for viral titre (panel A), systemic symptoms (panel B) and respiratory symptoms (panel C) 11
12 Figure S3. Goodness of fit plots for viral titre (left), systemic symptoms (centre) and respiratory symptoms (right). The upper panels represent the observations vs. predictions, the black line is the line of identity, the blue line the spline and show that the model fits well the data. The central panels represent the Normalised prediction distribution errors (NPDE) vs. the predictions and show no trends in the distribution of the residual depending on the predictions and the lower panels represent the NPDE vs. time and show no trends in the distribution of the residual depending on the sample time. Overall these goodness of fit plots suggest that the model fits well the data. 12
13 Figure S4. Effect of parameters b and s on the infectiousness rank. Infectiousness is computed as the area under the curve of J(t) = V(t)(b + sr(t)). Rank #1 is for the most infectious subject. Each line shows the evolution of a subject infectiousness rank for the different scenarios tested. 13
14 Monolix code: DESCRIPTION: Influenza viral kinetic, systemic and respiratory symptoms dynamics INPUT: parameter = {V0, beta, delta, p, psi, eta, c, sigma, alpha, gamma_r, alphasys, gamma_s} ; Parameters to estimate EQUATION: odetype = stiff ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ; compute the initial values t0 = 0 Tt_0 = 4e8 ; Fixing initial number of target cells to 4x10^8 I_0 = 0 V_0 = V0 ; Estimating initial viral titre F_0 = 0 F_s = 0 N_0 = 0 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ; System of ODE ddt_tt = -beta*tt*v ; Target cells ddt_i = beta*tt*v - delta*i -eta*i*n ; Infected cells ddt_v = p/(1+psi*f)*i - c*v sigma*beta*tt*v ; Free virus ddt_f = gamma_s*i - alpha*f ; Local pro-inflammatory cytokines ddt_fs = gamma_r*f-alphasys*fs ; Systemic pro-inflammatory cytokines ddt_n = Fs ; Cytotoxicity LVL =max(log10(v),0) ; Log-transformed viral titre logists =12*(Fs/(1+Fs)) ; Transformation of systemic pro-inflammatory cytokines levels into systemic symptom score logistr = 12*(F/(1+F)) ; Transformation of local pro-inflammatory cytokines levels into respiratory symptom score OUTPUT: output = {LVL, logists, logistr} ; Fitting simultaneously viral titre and symptoms scores 14
F. 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 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 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 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 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 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 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 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 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 informationBayesian Inference on Joint Mixture Models for Survival-Longitudinal Data with Multiple Features. Yangxin Huang
Bayesian Inference on Joint Mixture Models for Survival-Longitudinal Data with Multiple Features Yangxin Huang Department of Epidemiology and Biostatistics, COPH, USF, Tampa, FL yhuang@health.usf.edu January
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 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 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 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 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 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 informationInferring biomarkers for Mycobacterium avium subsp. paratuberculosis infection and disease progression in cattle using experimental data
Inferring biomarkers for Mycobacterium avium subsp. paratuberculosis infection and disease progression in cattle using experimental data 1,2,4, Gesham Magombedze, Tinevimbo Shiri, Shigetoshi Eda 4, 5,
More informationNonlinear mixed-effects models using Stata
Nonlinear mixed-effects models using Stata Yulia Marchenko Executive Director of Statistics StataCorp LP 2017 German Stata Users Group meeting Yulia Marchenko (StataCorp) 1 / 48 Outline What is NLMEM?
More informationEarly Detection and Assessment of Epidemics by Particle Filtering
Early Detection and Assessment of Epidemics by Particle Filtering C. Jégat, F. Carrat, C. Lajaunie and H. Wackernagel Abstract Influenza infects from 5% to 20% of the population during an epidemic episode,
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 informationMA/ST 810 Mathematical-Statistical Modeling and Analysis of Complex Systems
MA/ST 810 Mathematical-Statistical Modeling and Analysis of Complex Systems Integrating Mathematical and Statistical Models Recap of mathematical models Models and data Statistical models and sources of
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 informationBetween-Subject and Within-Subject Model Mixtures for Classifying HIV Treatment Response
Progress in Applied Mathematics Vol. 4, No. 2, 2012, pp. [148 166] DOI: 10.3968/j.pam.1925252820120402.S0801 ISSN 1925-251X [Print] ISSN 1925-2528 [Online] www.cscanada.net www.cscanada.org Between-Subject
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 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 informationMixed effect model for the spatiotemporal analysis of longitudinal manifold value data
Mixed effect model for the spatiotemporal analysis of longitudinal manifold value data Stéphanie Allassonnière with J.B. Schiratti, O. Colliot and S. Durrleman Université Paris Descartes & Ecole Polytechnique
More informationA Regression Model For Recurrent Events With Distribution Free Correlation Structure
A Regression Model For Recurrent Events With Distribution Free Correlation Structure J. Pénichoux(1), A. Latouche(2), T. Moreau(1) (1) INSERM U780 (2) Université de Versailles, EA2506 ISCB - 2009 - Prague
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 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 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 informationBayesian Regression (1/31/13)
STA613/CBB540: Statistical methods in computational biology Bayesian Regression (1/31/13) Lecturer: Barbara Engelhardt Scribe: Amanda Lea 1 Bayesian Paradigm Bayesian methods ask: given that I have observed
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 informationTransmission of Hand, Foot and Mouth Disease and Its Potential Driving
Transmission of Hand, Foot and Mouth Disease and Its Potential Driving Factors in Hong Kong, 2010-2014 Bingyi Yang 1, Eric H. Y. Lau 1*, Peng Wu 1, Benjamin J. Cowling 1 1 WHO Collaborating Centre for
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 informationPreface Introduction to Statistics and Data Analysis Overview: Statistical Inference, Samples, Populations, and Experimental Design The Role of
Preface Introduction to Statistics and Data Analysis Overview: Statistical Inference, Samples, Populations, and Experimental Design The Role of Probability Sampling Procedures Collection of Data Measures
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 informationApproximate Bayesian Computation
Approximate Bayesian Computation Michael Gutmann https://sites.google.com/site/michaelgutmann University of Helsinki and Aalto University 1st December 2015 Content Two parts: 1. The basics of approximate
More informationMedian Cross-Validation
Median Cross-Validation Chi-Wai Yu 1, and Bertrand Clarke 2 1 Department of Mathematics Hong Kong University of Science and Technology 2 Department of Medicine University of Miami IISA 2011 Outline Motivational
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 information1 Disease Spread Model
Technical Appendix for The Impact of Mass Gatherings and Holiday Traveling on the Course of an Influenza Pandemic: A Computational Model Pengyi Shi, Pinar Keskinocak, Julie L Swann, Bruce Y Lee December
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 informationMixed-Effects Biological Models: Estimation and Inference
Mixed-Effects Biological Models: Estimation and Inference Hulin Wu, Ph.D. Dr. D.R. Seth Family Professor & Associate Chair Department of Biostatistics School of Public Health University of Texas Health
More informationAsymptotic equivalence of paired Hotelling test and conditional logistic regression
Asymptotic equivalence of paired Hotelling test and conditional logistic regression Félix Balazard 1,2 arxiv:1610.06774v1 [math.st] 21 Oct 2016 Abstract 1 Sorbonne Universités, UPMC Univ Paris 06, CNRS
More informationStatistical Methods for Bridging Experimental Data and Dynamic Models with Biomedical Applications
Statistical Methods for Bridging Experimental Data and Dynamic Models with Biomedical Applications Hulin Wu, Ph.D. Dr. D.R. Seth Family Professor & Associate Chair Department of Biostatistics, School of
More informationIntroduction to lnmle: An R Package for Marginally Specified Logistic-Normal Models for Longitudinal Binary Data
Introduction to lnmle: An R Package for Marginally Specified Logistic-Normal Models for Longitudinal Binary Data Bryan A. Comstock and Patrick J. Heagerty Department of Biostatistics University of Washington
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 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 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 informationOn the Spread of Epidemics in a Closed Heterogeneous Population
On the Spread of Epidemics in a Closed Heterogeneous Population Artem Novozhilov Applied Mathematics 1 Moscow State University of Railway Engineering (MIIT) the 3d Workshop on Mathematical Models and Numerical
More informationStochastic Analogues to Deterministic Optimizers
Stochastic Analogues to Deterministic Optimizers ISMP 2018 Bordeaux, France Vivak Patel Presented by: Mihai Anitescu July 6, 2018 1 Apology I apologize for not being here to give this talk myself. I injured
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 informationInflammation and organ failure severely affect midazolam clearance in critically ill children
Online Data Supplement Inflammation and organ failure severely affect midazolam clearance in critically ill children Nienke J. Vet, MD 1,2 * Janneke M. Brussee, MSc 3* Matthijs de Hoog, MD, PhD 1,2 Miriam
More informationEfficient Likelihood-Free Inference
Efficient Likelihood-Free Inference Michael Gutmann http://homepages.inf.ed.ac.uk/mgutmann Institute for Adaptive and Neural Computation School of Informatics, University of Edinburgh 8th November 2017
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 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 informationA Hierarchical Bayesian Approach for Parameter Estimation in HIV Models
A Hierarchical Bayesian Approach for Parameter Estimation in HIV Models H. T. Banks, Sarah Grove, Shuhua Hu, and Yanyuan Ma Center for Research in Scientific Computation North Carolina State University
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 informationA State Space Model Approach for HIV Infection Dynamics
A State Space Model Approach for HIV Infection Dynamics Jiabin Wang a, Hua Liang b, and Rong Chen a Mathematical models have been proposed and developed to model HIV dynamics over the past few decades,
More informationCSE 473: Artificial Intelligence Autumn Topics
CSE 473: Artificial Intelligence Autumn 2014 Bayesian Networks Learning II Dan Weld Slides adapted from Jack Breese, Dan Klein, Daphne Koller, Stuart Russell, Andrew Moore & Luke Zettlemoyer 1 473 Topics
More informationSpecifying Latent Curve and Other Growth Models Using Mplus. (Revised )
Ronald H. Heck 1 University of Hawai i at Mānoa Handout #20 Specifying Latent Curve and Other Growth Models Using Mplus (Revised 12-1-2014) The SEM approach offers a contrasting framework for use in analyzing
More informationDSGE Methods. Estimation of DSGE models: GMM and Indirect Inference. Willi Mutschler, M.Sc.
DSGE Methods Estimation of DSGE models: GMM and Indirect Inference Willi Mutschler, M.Sc. Institute of Econometrics and Economic Statistics University of Münster willi.mutschler@wiwi.uni-muenster.de Summer
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 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 informationBasic Concepts in Population Modeling, Simulation, and Model-Based Drug Development Part 2: Introduction to Pharmacokinetic Modeling Methods
Tutorial Citation: CPT: Pharmacometrics & Systems Pharmacology (3), e38; doi:.38/psp.3.4 3 ASCPT All rights reserved 63-836/ www.nature.com/psp Basic Concepts in Population Modeling, Simulation, and Model-Based
More informationPMR Learning as Inference
Outline PMR Learning as Inference Probabilistic Modelling and Reasoning Amos Storkey Modelling 2 The Exponential Family 3 Bayesian Sets School of Informatics, University of Edinburgh Amos Storkey PMR Learning
More informationInference for partially observed stochastic dynamic systems: A new algorithm, its theory and applications
Inference for partially observed stochastic dynamic systems: A new algorithm, its theory and applications Edward Ionides Department of Statistics, University of Michigan ionides@umich.edu Statistics Department
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 informationStat 5101 Lecture Notes
Stat 5101 Lecture Notes Charles J. Geyer Copyright 1998, 1999, 2000, 2001 by Charles J. Geyer May 7, 2001 ii Stat 5101 (Geyer) Course Notes Contents 1 Random Variables and Change of Variables 1 1.1 Random
More informationAnalysis of Longitudinal Data. Patrick J. Heagerty PhD Department of Biostatistics University of Washington
Analysis of Longitudinal Data Patrick J Heagerty PhD Department of Biostatistics University of Washington Auckland 8 Session One Outline Examples of longitudinal data Scientific motivation Opportunities
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 informationMulti-level Models: Idea
Review of 140.656 Review Introduction to multi-level models The two-stage normal-normal model Two-stage linear models with random effects Three-stage linear models Two-stage logistic regression with random
More informationEarly Detection of a Change in Poisson Rate After Accounting For Population Size Effects
Early Detection of a Change in Poisson Rate After Accounting For Population Size Effects School of Industrial and Systems Engineering, Georgia Institute of Technology, 765 Ferst Drive NW, Atlanta, GA 30332-0205,
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 informationMaster's thesis Longitudinal modeling of antibody dynamics during herpes zoster infection
2014 2015 FACULTY OF SCIENCES Master of Statistics Master's thesis Longitudinal modeling of antibody dynamics during herpes zoster infection Supervisor : Prof. dr. Niel HENS Supervisor : Dr. BENSON OGUNJIMI
More informationComparison of Stochastic Soybean Yield Response Functions to Phosphorus Fertilizer
"Science Stays True Here" Journal of Mathematics and Statistical Science, Volume 016, 111 Science Signpost Publishing Comparison of Stochastic Soybean Yield Response Functions to Phosphorus Fertilizer
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 informationNon-Parametric Bayes
Non-Parametric Bayes Mark Schmidt UBC Machine Learning Reading Group January 2016 Current Hot Topics in Machine Learning Bayesian learning includes: Gaussian processes. Approximate inference. Bayesian
More informationSYDE 372 Introduction to Pattern Recognition. Probability Measures for Classification: Part I
SYDE 372 Introduction to Pattern Recognition Probability Measures for Classification: Part I Alexander Wong Department of Systems Design Engineering University of Waterloo Outline 1 2 3 4 Why use probability
More informationMODELING THE INTERACTION OF CYTOTOXIC T LYMPHOCYTES AND INFLUENZA VIRUS INFECTED EPITHELIAL CELLS. Abdessamad Tridane. Yang Kuang
MATHEMATICAL BIOSCIENCES doi:.3934/mbe.2.7.75 AND ENGINEERING Volume 7, Number, January 2 pp. 75 89 MODELING THE INTERACTION OF CYTOTOXIC T LYMPHOCYTES AND INFLUENZA VIRUS INFECTED EPITHELIAL CELLS Abdessamad
More informationStatistical Inference for Stochastic Epidemic Models
Statistical Inference for Stochastic Epidemic Models George Streftaris 1 and Gavin J. Gibson 1 1 Department of Actuarial Mathematics & Statistics, Heriot-Watt University, Riccarton, Edinburgh EH14 4AS,
More information. Also, in this case, p i = N1 ) T, (2) where. I γ C N(N 2 2 F + N1 2 Q)
Supplementary information S7 Testing for association at imputed SPs puted SPs Score tests A Score Test needs calculations of the observed data score and information matrix only under the null hypothesis,
More informationModelling geoadditive survival data
Modelling geoadditive survival data Thomas Kneib & Ludwig Fahrmeir Department of Statistics, Ludwig-Maximilians-University Munich 1. Leukemia survival data 2. Structured hazard regression 3. Mixed model
More informationEMPIRICAL DATA SETS in human fatigue and performance. Nonlinear Mixed-Effects Modeling: Individualization and Prediction
onlinear Mixed-Effects Modeling: Individualization and Prediction Erik Olofsen, David F. Dinges, and Hans P.A. Van Dongen OLOFSE E, DIGES DF, VA DOGE HPA. onlinear mixedeffects modeling: individualization
More informationMachine Learning Linear Classification. Prof. Matteo Matteucci
Machine Learning Linear Classification Prof. Matteo Matteucci Recall from the first lecture 2 X R p Regression Y R Continuous Output X R p Y {Ω 0, Ω 1,, Ω K } Classification Discrete Output X R p Y (X)
More informationJoint Modeling of Longitudinal Item Response Data and Survival
Joint Modeling of Longitudinal Item Response Data and Survival Jean-Paul Fox University of Twente Department of Research Methodology, Measurement and Data Analysis Faculty of Behavioural Sciences Enschede,
More informationStatistics 203: Introduction to Regression and Analysis of Variance Course review
Statistics 203: Introduction to Regression and Analysis of Variance Course review Jonathan Taylor - p. 1/?? Today Review / overview of what we learned. - p. 2/?? General themes in regression models Specifying
More informationPlausible Values for Latent Variables Using Mplus
Plausible Values for Latent Variables Using Mplus Tihomir Asparouhov and Bengt Muthén August 21, 2010 1 1 Introduction Plausible values are imputed values for latent variables. All latent variables can
More informationPackage LBLGXE. R topics documented: July 20, Type Package
Type Package Package LBLGXE July 20, 2015 Title Bayesian Lasso for detecting Rare (or Common) Haplotype Association and their interactions with Environmental Covariates Version 1.2 Date 2015-07-09 Author
More informationRate Maps and Smoothing
Rate Maps and Smoothing Luc Anselin Spatial Analysis Laboratory Dept. Agricultural and Consumer Economics University of Illinois, Urbana-Champaign http://sal.agecon.uiuc.edu Outline Mapping Rates Risk
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 informationCost Efficiency, Asymmetry and Dependence in US electricity industry.
Cost Efficiency, Asymmetry and Dependence in US electricity industry. Graziella Bonanno bonanno@diag.uniroma1.it Department of Computer, Control, and Management Engineering Antonio Ruberti - Sapienza University
More informationChallenges in modelling the pharmacokinetics of isoniazid in South African tuberculosis patients
PAGE 2005 Pamplona, Spain Challenges in modelling the pharmacokinetics of isoniazid in South African tuberculosis patients Justin J Wilkins 1, Grant Langdon 1, Helen McIlleron 1, Goonaseelan Pillai 2,
More informationAnalysing repeated measurements whilst accounting for derivative tracking, varying within-subject variance and autocorrelation: the xtiou command
Analysing repeated measurements whilst accounting for derivative tracking, varying within-subject variance and autocorrelation: the xtiou command R.A. Hughes* 1, M.G. Kenward 2, J.A.C. Sterne 1, K. Tilling
More informationBIAS OF MAXIMUM-LIKELIHOOD ESTIMATES IN LOGISTIC AND COX REGRESSION MODELS: A COMPARATIVE SIMULATION STUDY
BIAS OF MAXIMUM-LIKELIHOOD ESTIMATES IN LOGISTIC AND COX REGRESSION MODELS: A COMPARATIVE SIMULATION STUDY Ingo Langner 1, Ralf Bender 2, Rebecca Lenz-Tönjes 1, Helmut Küchenhoff 2, Maria Blettner 2 1
More informationPENALIZING YOUR MODELS
PENALIZING YOUR MODELS AN OVERVIEW OF THE GENERALIZED REGRESSION PLATFORM Michael Crotty & Clay Barker Research Statisticians JMP Division, SAS Institute Copyr i g ht 2012, SAS Ins titut e Inc. All rights
More informationBayesian Models in Machine Learning
Bayesian Models in Machine Learning Lukáš Burget Escuela de Ciencias Informáticas 2017 Buenos Aires, July 24-29 2017 Frequentist vs. Bayesian Frequentist point of view: Probability is the frequency of
More informationLearning Bayesian belief networks
Lecture 4 Learning Bayesian belief networks Milos Hauskrecht milos@cs.pitt.edu 5329 Sennott Square Administration Midterm: Monday, March 7, 2003 In class Closed book Material covered by Wednesday, March
More informationLinear models for the joint analysis of multiple. array-cgh profiles
Linear models for the joint analysis of multiple array-cgh profiles F. Picard, E. Lebarbier, B. Thiam, S. Robin. UMR 5558 CNRS Univ. Lyon 1, Lyon UMR 518 AgroParisTech/INRA, F-75231, Paris Statistics for
More informationA simple kinetic equation of swarm formation: blow up and global existence
A simple kinetic equation of swarm formation: blow up and global existence Miroslaw Lachowicz, Henryk Leszczyński, Martin Parisot To cite this version: Miroslaw Lachowicz, Henryk Leszczyński, Martin Parisot.
More information