Semiparametric Mixed Effects Models with Flexible Random Effects Distribution
|
|
- Rudolf Rodgers
- 5 years ago
- Views:
Transcription
1 Semiparametric Mixed Effects Models with Flexible Random Effects Distribution Marie Davidian North Carolina State University davidian Joint work with A. Tsiatis, J. Chen, X. Song, and D. Zhang
2 Outline 1. Introduction and motivation 2. Semiparametric mixed models 3. The class H and the SNP representation 4. Implementation 5. Simulations 6. Examples, revisited 7. Discussion
3 1. Introduction and motivation Longitudinal studies: Clinical trials, epidemiological investigations Repeated measures on some variable on each subject intermittently over time Survival endpoint (e.g., death, time to disease progression) Objectives: Inference on Within-subject patterns of change of repeated measurement variable, association with covariates Relationship between repeated measurement variable and survival, association with covariates
4 Example 1: Framingham study Cholesterol measurements every 2 years for 2634 participants over 10-year period Objectives: Change in cholesterol over time and association with age at baseline, gender Consider a subset of 200 subjects for illustrative purposes Individual profiles: Approximate straight line
5 Time in years Cholesterol level Cholesterol levels over time
6 Standard model: Linear mixed effects model For subject i at time t ij, j = 1,..., n i, i = 1,..., m Y ij = β 1 age i + β 2 gender i + β 3 age i t ij + β 4 gender i t ij + b 0i + b 1i t ij + e ij Usual assumptions: e i = (e i1,..., e ini ) T N(0, σ 2 I) b i = (b 0i, b 1i ) T N(µ b, D) Relevance of assumptions: Individual regression fits Pooled residuals, subject-specific slopes appear normal However, subject-specific intercepts do not...
7 Percentage Estimates of subject specific intercepts
8 Example 2: Six cities study Binary indicator of respiratory infection recored annually for 537 Ohio children at ages 7 10 Objectives: Within-subject change in respiratory status, association with baseline maternal smoking Standard model: Generalized linear mixed effects model For subject i at age a ij, j = 1,..., n i, i = 1,..., m E(Y ij b i ) = exp(β 1smoke + β 2 a ij + b i ) 1 + exp(β 1 smoke + β 2 a ij + b i ), var(y ij b i ) = E(Y ij b i ){1 E(Y ij b i )} Usual assumption: b i N(µ b, σ 2 b )
9 Example 3: ACTG 175 Clinical trial with 2467 subjects to compare 4 antiretroviral regimens Main objective: Compare on basis of time to AIDS or death Also, CD4 counts approximately every 12 weeks Subsequent objective: Characterize within-subject patterns of CD4 change complicated by informative censoring Subsequent objective: Characterize relationship between features of CD4 profiles and survival
10 log CD week Death/progression, censoring throughout
11 Standard model: Joint mixed effects-proportional hazards model Longitudinal data model: For subject i at weeks t i = (t i1,..., t ini ) T W ij = X i (t ij ) + e ij, X i (u) = b 0i + b 1i u X i (u) = inherent trajectory e i = (e i1,..., e ini ) T N(0, σ 2 I) Survival data model: For subject i, observe V i = min(t i, C i ), i = I(T i C i ), covariates S i λ i (u) = lim du 0 du 1 P {u T i < u + du T i u, α i, S i, C i, e i, t i } = λ 0 (u) exp{γx i (u) + η T S i } Usual assumption: b i = (b i0, b i1 ) T N(µ b, D)
12 intercept slope
13 2. Semiparametric mixed models Theme: The foregoing examples suggest that A simple parametric model may be adequate to describe subject-specific profiles in terms of random effects b i However, the relevance of the usual normality assumption on random effects is questionable Concern: Sensitivity of inferences to departures from normality
14 Needed: Relax the normality assumption on b i Semiparametric model Completely nonparametric (e.g., Mallet, 1986; Butler and Louis, 1992) includes unusual, discrete distributions Alternatively: Impose some realistic yet not overly restrictive conditions Restrict to a smooth class (e.g., Davidian and Gallant, 1993; Madger and Zeger, 1996; Verbeke and Lesaffre, 1996; Tao et al., 1999) Here: Assume b i have distribution with density in a smooth class H
15 3. The class H and the SNP representation Assume: b i = g(µ, S i ) + RZ i, Z i has density h H, R lower triangular with distinct elements θ E.g., for ACTG 175, g(µ, S i ) = µ 0 (1 S i ) + µ 1 S i S i = I(Trt=ZDV) H is a class of smooth densities studied by Gallant and Nychka (1987) Densities in H are sufficiently differentiable to rule out kinks, jumps, violent oscillation But can be skewed, multi-modal, fat- or thin-tailed relative to the normal (and the normal is H)
16 Formally: h H for Z (q 1) may be written as h(z) = P (z)ϕ 2 q (z) + small lower bound for tail behavior P (z) is an infinite-dimensional polynomial ϕ q (z) is q-variate standard normal density Practically speaking: Suggests approximating h H by truncation h K (z) = PK(z)ϕ 2 q (z) P K (z) is Kth order polynomial; e.g., for K = 2 P K (z) = a 00 + a 10 z 1 + a 01 z 2 + a 20 z a 02 z a 11 z 1 z 2 Vector of coefficients a must satisfy h K (z) dz = 1 K = 0 is standard normal b i N{g(µ, S i ), RR T }
17 Imposing h K (z) dz = 1: a (d 1), d depends on K h K (z) dz = 1 E{P 2 K(U)} = 1, U N(0, I) E{P 2 K (U)} = at Aa (A p.d.) = a T BBa = c T c = 1, c = Ba Polar coordinate transformation c 1 = sin(φ 1 ), c 2 = cos(φ 1 ) sin(φ 2 ),. c d 1 = cos(φ 1 ) cos(φ 2 ) cos(φ d 2 ) sin(φ d 1 ), c d = cos(φ 1 ) cos(φ 2 ) cos(φ d 2 ) cos(φ d 1 ), π/2 < φ r π/2, r = 1,..., d 1. Parameterize h K (z) in terms of φ = (φ 1,..., φ d 1 ) T.
18 Result: For fixed K, may represent density of b i in terms of (µ T, θ T, φ T ) T Likelihood for Ω = (µ T, θ T, φ T ) T plus any other model parameters (e.g., β, γ, η) is usual, finite-dimensional problem In principle, can use standard optimization methods to estimate Ω (coming up... ) Lingo: Seminonparametric Choosing tuning parameter K: K controls degree of flexibility and departure from normality (like a bandwidth )
19 Adaptive choice of K based on information criteria: If l K ( Ω) is maximized loglikelihood for fixed K, N = total number of observations, Ω (p 1), minimize { l K ( Ω) + pc(n)}/n AIC, C(N) = 1; BIC, C(N) = log N/2; Hannan-Quinn (HQ), C(N) = log log N AIC prefers larger models, BIC smaller, HQ intermediate Confidence intervals fixing K at choice achieve nominal coverage (Eastwood and Gallant, 1991)
20 4. Implementation Linear mixed effects model: For normal e i, Ω = (µ T, θ T, φ T, β T, σ) T, can write loglikelihood l K (Ω; Y ) in a closed form Maximize in Ω using standard optimization routines, e.g., SAS nlpqn Starting values chosen by grid search or penalized loglikelihood SEs, confidence intervals usual inverse of observed information for chosen K Zhang and Davidian (2001, Biometrics)
21 Generalized linear mixed effects model: Other models (e.g., binomial, Poisson), Ω = (µ T, θ T, φ T, β T ) T, cannot write l K (Ω; Y ) in a closed form Gallant and Tauchen (1992) provide efficient rejection sampling algorithm from estimated h K (z) (acceptance rate > 50%) Facilitates use of MCEM algorithm (e.g., McCulloch, 1997; Booth and Hobert, 1999) SEs, confidence intervals MC approx observed information for chosen K Chen, Zhang, and Davidian (2002, Biostatistics)
22 Joint longitudinal-survival model: Ω = (µ T, θ T, φ T, γ, η T, λ 0 ) T Under assumptions, l K (Ω; V,, W, t, Z) = log L(Ω; V,, W, t, Z) n m i L(Ω; V,, W, t, S) = p(v i, i b i, S i, γ, η, λ 0 ) p(w ij b i, σ 2, t ij ) i=1 p(b i Z i, µ, θ, φ)db i j=1 p(w ij b i, σ 2, t ij ) = { 1 exp (W ij b i0 b i1 t ij ) 2 } 2πσ 2 2σ 2, p(v i, i b i, S i, γ, η, λ 0 ) = [λ 0 (V i ) exp{γ(b i0 + b i1 V i ) + ηz i }] i [ exp Vi 0 λ 0 (u) exp{γ(α i0 + α i1 u) + ηs i }du, p(b i S i, µ, θ, φ) = h K [R 1 {b i g(µ, S i )}] R 1 ]
23 EM algorithm: L(Ω; V,, W, t, Z) is maximized when λ 0 (u) is non-zero only at death times, and Ω maximizing L(Ω; V,, W, t, Z) exists E-step: Intractable integration carried out via Gauss-Hermite quadrature M-step: Maximization in (µ T, θ T, φ T ) T and (γ, η T, σ, λ 0 ) T separates One-step Newton-Raphson update for (γ, η T ) T SEs and confidence intervals: Profile likelihood Song, Davidian, and Tsiatis (2002, Biometrics)
24 5. Simulations Linear mixed model: 100 data sets, each fit with K = 0, 1, 2 Y ij = b i + β 1 t ij + β 2 w i + e ij, i = 1,..., 100, j = 1,..., 5 t ij = j 3, β 1 = 2, w i = I(i 50), β 2 = 1, e ij N(0, ) Case 1: b i 0.7N( 3, 1) + 0.3N(2, 1) (mixture of normals) AIC preferred K = 1, 2 35%, 65% of time (BIC: 76%, 24%; HQ: 56%, 44%) Case 2: b i N( 1.5, 6.25) AIC preferred K = 0, 1, 2 84%, 7%, 9% (BIC: 97%, 3%, 0%; HQ: 89%, 5%, 6%)
25 K = 0 (Normality) Preferred by HQ MC Ave. MC SD Ave. SE MC Ave. MC SD Ave. SE RE Case 1: Mixture Scenario β 1 (2) β 2 (1) E(b) ( 1.5) var(b) (6.25) σ (0.5) Case 2: Normal Scenario β 1 (2) β 2 (1) E(b) ( 1.5) var(b) (6.25) σ (0.5)
26 (a) (b) Densities Densities x x (a) Average of 100 estimates (K = 0 and AIC, BIC, HQ) and truth (b) 100 estimates chosen by HQ
27 Joint model: 200 data sets, each fit with K = 0, 1, 2, 2-point quadrature in E-step W ij = b 0i + b 01 t ij + e ij, i = 1,..., 200, e ij N(0, 0.60) t ij = 0, 2, 4, 8, 16, 24, 32, 40, 48, 56, 64, 72, 80, 10% missingness {var(b i0 ), cov(b i0, b i1 ), var(b i1 )} = (4.96, , 0.012) λ i (u) λ 0 (u) exp{γx i (u)}, λ 0 (u) 1 for u 16, = 0 ow; γ = 1.0 Exponential (110) censoring (53% censored data) Case 1: b i bivariate mixture of normals Case 2: b i bivariate normal
28 Estimation of b i distributions: Similar to linear mixed model Estimation of hazard parameter γ: K = 0 AIC HQ BIC Case 1: Mixture Scenario γ SD SE CP %K = (0, 1, 2) (0.0,70.5,29.5) (0.0,93.5,6.5) (0.0,100.0,0.0) Case 2: Normal Scenario γ SD SE CP %K = (0, 1, 2) (91.5,5.0,3.5) (99.0,1.0,0.5) (100.0,0.0,0.0) Amazing robustness to distribution of b i
29 density density b b b b
30 density b0 density b1
31 density b0 density b1
32 6. Examples, revisited Example 1: Framingham cholesterol data Y ij = β 1 age i + β 2 gender i + β 3 age i t ij + β 4 gender i t ij + b 0i + b 1i t ij + e ij e i N(0, σ 2 I), b i = µ + RZ i, h H Subset of 200 subjects, each every 2 years Fit using K = 0, 1, 2 All criteria select K = 1
33 0 2 4 Density Slope Intercept Estimated joint density for K = 1
34 Density Intercept Estimated marginals for intercepts, K = 0 and 1
35 Example 3: ACTG 175 joint longitudinal-survival model, W ij = X i (t ij ) + e ij, X i (u) = b 0i + b 1i u, e i N(0, σ 2 ) b i = µ 0 (1 S i ) + µ 1 S i + RZ i, S i = I(Trt=ZDV) h H λ i (u) = lim du 0 du 1 P {u T i < u + du T i u, α i, S i, C i, e i, t i } = λ 0 (u) exp{γx i (u) + η T S i } 2467 subjects Fit for K = 0, 1, 2, 3, 4; K = 3 or 4 chosen
36 800 0 density b1 Estimated joint density for K = b0
37 density b1 Estimated marginals for slope for K = 0, 2, 3, 4
38 7. Discussion Potential to gain efficiency in parameters associated with subject-level covariates in longitudinal models; other parameters robust Remarkable robustness of inference on hazard parameters to misspecification of random effects distribution in joint model Insight on population heterogeneity, possible omitted covariates Implementation only mildly more difficult than assuming normal random effects
SEMIPARAMETRIC APPROACHES TO INFERENCE IN JOINT MODELS FOR LONGITUDINAL AND TIME-TO-EVENT DATA
SEMIPARAMETRIC APPROACHES TO INFERENCE IN JOINT MODELS FOR LONGITUDINAL AND TIME-TO-EVENT DATA Marie Davidian and Anastasios A. Tsiatis http://www.stat.ncsu.edu/ davidian/ http://www.stat.ncsu.edu/ tsiatis/
More informationSome New Methods for Latent Variable Models and Survival Analysis. Latent-Model Robustness in Structural Measurement Error Models.
Some New Methods for Latent Variable Models and Survival Analysis Marie Davidian Department of Statistics North Carolina State University 1. Introduction Outline 3. Empirically checking latent-model robustness
More informationMixed model analysis of censored longitudinal data with flexible random-effects density
Biostatistics (2012), 13, 1, pp. 61 73 doi:10.1093/biostatistics/kxr026 Advance Access publication on September 13, 2011 Mixed model analysis of censored longitudinal data with flexible random-effects
More informationA Sampling of IMPACT Research:
A Sampling of IMPACT Research: Methods for Analysis with Dropout and Identifying Optimal Treatment Regimes Marie Davidian Department of Statistics North Carolina State University http://www.stat.ncsu.edu/
More informationImproving Efficiency of Inferences in Randomized Clinical Trials Using Auxiliary Covariates
Improving Efficiency of Inferences in Randomized Clinical Trials Using Auxiliary Covariates Anastasios (Butch) Tsiatis Department of Statistics North Carolina State University http://www.stat.ncsu.edu/
More informationLatent-model Robustness in Joint Models for a Primary Endpoint and a Longitudinal Process
Latent-model Robustness in Joint Models for a Primary Endpoint and a Longitudinal Process Xianzheng Huang, 1, * Leonard A. Stefanski, 2 and Marie Davidian 2 1 Department of Statistics, University of South
More informationUsing Estimating Equations for Spatially Correlated A
Using Estimating Equations for Spatially Correlated Areal Data December 8, 2009 Introduction GEEs Spatial Estimating Equations Implementation Simulation Conclusion Typical Problem Assess the relationship
More informationDynamic Prediction of Disease Progression Using Longitudinal Biomarker Data
Dynamic Prediction of Disease Progression Using Longitudinal Biomarker Data Xuelin Huang Department of Biostatistics M. D. Anderson Cancer Center The University of Texas Joint Work with Jing Ning, Sangbum
More informationMultivariate Survival Analysis
Multivariate Survival Analysis Previously we have assumed that either (X i, δ i ) or (X i, δ i, Z i ), i = 1,..., n, are i.i.d.. This may not always be the case. Multivariate survival data can arise in
More informationDiscussion of Missing Data Methods in Longitudinal Studies: A Review by Ibrahim and Molenberghs
Discussion of Missing Data Methods in Longitudinal Studies: A Review by Ibrahim and Molenberghs Michael J. Daniels and Chenguang Wang Jan. 18, 2009 First, we would like to thank Joe and Geert for a carefully
More informationLongitudinal + Reliability = Joint Modeling
Longitudinal + Reliability = Joint Modeling Carles Serrat Institute of Statistics and Mathematics Applied to Building CYTED-HAROSA International Workshop November 21-22, 2013 Barcelona Mainly from Rizopoulos,
More informationApproximation of Survival Function by Taylor Series for General Partly Interval Censored Data
Malaysian Journal of Mathematical Sciences 11(3): 33 315 (217) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES Journal homepage: http://einspem.upm.edu.my/journal Approximation of Survival Function by Taylor
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 informationREGRESSION ANALYSIS FOR TIME-TO-EVENT DATA THE PROPORTIONAL HAZARDS (COX) MODEL ST520
REGRESSION ANALYSIS FOR TIME-TO-EVENT DATA THE PROPORTIONAL HAZARDS (COX) MODEL ST520 Department of Statistics North Carolina State University Presented by: Butch Tsiatis, Department of Statistics, NCSU
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 informationThe consequences of misspecifying the random effects distribution when fitting generalized linear mixed models
The consequences of misspecifying the random effects distribution when fitting generalized linear mixed models John M. Neuhaus Charles E. McCulloch Division of Biostatistics University of California, San
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 informationBiost 518 Applied Biostatistics II. Purpose of Statistics. First Stage of Scientific Investigation. Further Stages of Scientific Investigation
Biost 58 Applied Biostatistics II Scott S. Emerson, M.D., Ph.D. Professor of Biostatistics University of Washington Lecture 5: Review Purpose of Statistics Statistics is about science (Science in the broadest
More informationModels for Multivariate Panel Count Data
Semiparametric Models for Multivariate Panel Count Data KyungMann Kim University of Wisconsin-Madison kmkim@biostat.wisc.edu 2 April 2015 Outline 1 Introduction 2 3 4 Panel Count Data Motivation Previous
More informationPower and Sample Size Calculations with the Additive Hazards Model
Journal of Data Science 10(2012), 143-155 Power and Sample Size Calculations with the Additive Hazards Model Ling Chen, Chengjie Xiong, J. Philip Miller and Feng Gao Washington University School of Medicine
More informationSmooth Semiparametric Regression Analysis for Arbitrarily Censored Time-to-Event Data
Smooth Semiparametric Regression Analysis for Arbitrarily Censored Time-to-Event Data Min Zhang and Marie Davidian Department of Statistics, North Carolina State University, Raleigh, North Carolina 27695-8203,
More informationIntroduction to Statistical Analysis
Introduction to Statistical Analysis Changyu Shen Richard A. and Susan F. Smith Center for Outcomes Research in Cardiology Beth Israel Deaconess Medical Center Harvard Medical School Objectives Descriptive
More informationFaculty of Health Sciences. Regression models. Counts, Poisson regression, Lene Theil Skovgaard. Dept. of Biostatistics
Faculty of Health Sciences Regression models Counts, Poisson regression, 27-5-2013 Lene Theil Skovgaard Dept. of Biostatistics 1 / 36 Count outcome PKA & LTS, Sect. 7.2 Poisson regression The Binomial
More informationA TWO-STAGE LINEAR MIXED-EFFECTS/COX MODEL FOR LONGITUDINAL DATA WITH MEASUREMENT ERROR AND SURVIVAL
A TWO-STAGE LINEAR MIXED-EFFECTS/COX MODEL FOR LONGITUDINAL DATA WITH MEASUREMENT ERROR AND SURVIVAL Christopher H. Morrell, Loyola College in Maryland, and Larry J. Brant, NIA Christopher H. Morrell,
More informationA Bayesian Nonparametric Approach to Monotone Missing Data in Longitudinal Studies with Informative Missingness
A Bayesian Nonparametric Approach to Monotone Missing Data in Longitudinal Studies with Informative Missingness A. Linero and M. Daniels UF, UT-Austin SRC 2014, Galveston, TX 1 Background 2 Working model
More informationChapter 1 Statistical Inference
Chapter 1 Statistical Inference causal inference To infer causality, you need a randomized experiment (or a huge observational study and lots of outside information). inference to populations Generalizations
More informationWeb-based Supplementary Materials for A Robust Method for Estimating. Optimal Treatment Regimes
Biometrics 000, 000 000 DOI: 000 000 0000 Web-based Supplementary Materials for A Robust Method for Estimating Optimal Treatment Regimes Baqun Zhang, Anastasios A. Tsiatis, Eric B. Laber, and Marie Davidian
More informationLikelihood and Conditional Likelihood Inference for Generalized Additive Mixed Models for Clustered Data
Likelihood and Conditional Likelihood Inference for Generalized Additive Mied Models for Clustered Data Daowen Zhang and Marie Davidian Department of Statistics, North Carolina State University Bo 8203,
More informationRegularization in Cox Frailty Models
Regularization in Cox Frailty Models Andreas Groll 1, Trevor Hastie 2, Gerhard Tutz 3 1 Ludwig-Maximilians-Universität Munich, Department of Mathematics, Theresienstraße 39, 80333 Munich, Germany 2 University
More informationHierarchical Hurdle Models for Zero-In(De)flated Count Data of Complex Designs
for Zero-In(De)flated Count Data of Complex Designs Marek Molas 1, Emmanuel Lesaffre 1,2 1 Erasmus MC 2 L-Biostat Erasmus Universiteit - Rotterdam Katholieke Universiteit Leuven The Netherlands Belgium
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 informationAnalysing geoadditive regression data: a mixed model approach
Analysing geoadditive regression data: a mixed model approach Institut für Statistik, Ludwig-Maximilians-Universität München Joint work with Ludwig Fahrmeir & Stefan Lang 25.11.2005 Spatio-temporal regression
More informationBiostatistics Workshop Longitudinal Data Analysis. Session 4 GARRETT FITZMAURICE
Biostatistics Workshop 2008 Longitudinal Data Analysis Session 4 GARRETT FITZMAURICE Harvard University 1 LINEAR MIXED EFFECTS MODELS Motivating Example: Influence of Menarche on Changes in Body Fat Prospective
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 informationGoodness-of-fit tests in mixed models DEPARTMENT OF DECISION SCIENCES AND INFORMATION MANAGEMENT (KBI) G. Claeskens and J.D. Hart
Faculty of Business and Economics Goodness-of-fit tests in mixed models G. Claeskens and J.D. Hart DEPARTMENT OF DECISION SCIENCES AND INFORMATION MANAGEMENT (KBI) KBI 0902 Goodness-of-fit tests in mixed
More informationPENALIZED LIKELIHOOD PARAMETER ESTIMATION FOR ADDITIVE HAZARD MODELS WITH INTERVAL CENSORED DATA
PENALIZED LIKELIHOOD PARAMETER ESTIMATION FOR ADDITIVE HAZARD MODELS WITH INTERVAL CENSORED DATA Kasun Rathnayake ; A/Prof Jun Ma Department of Statistics Faculty of Science and Engineering Macquarie University
More informationConditional Estimation for Generalized Linear Models When Covariates Are Subject-specific Parameters in a Mixed Model for Longitudinal Measurements
Conditional Estimation for Generalized Linear Models When Covariates Are Subject-specific Parameters in a Mixed Model for Longitudinal Measurements Erning Li, Daowen Zhang, and Marie Davidian Department
More informationFall 2017 STAT 532 Homework Peter Hoff. 1. Let P be a probability measure on a collection of sets A.
1. Let P be a probability measure on a collection of sets A. (a) For each n N, let H n be a set in A such that H n H n+1. Show that P (H n ) monotonically converges to P ( k=1 H k) as n. (b) For each n
More informationStat 542: Item Response Theory Modeling Using The Extended Rank Likelihood
Stat 542: Item Response Theory Modeling Using The Extended Rank Likelihood Jonathan Gruhl March 18, 2010 1 Introduction Researchers commonly apply item response theory (IRT) models to binary and ordinal
More informationMixed models in R using the lme4 package Part 7: Generalized linear mixed models
Mixed models in R using the lme4 package Part 7: Generalized linear mixed models Douglas Bates University of Wisconsin - Madison and R Development Core Team University of
More informationComputationally Efficient Estimation of Multilevel High-Dimensional Latent Variable Models
Computationally Efficient Estimation of Multilevel High-Dimensional Latent Variable Models Tihomir Asparouhov 1, Bengt Muthen 2 Muthen & Muthen 1 UCLA 2 Abstract Multilevel analysis often leads to modeling
More informationMultilevel Statistical Models: 3 rd edition, 2003 Contents
Multilevel Statistical Models: 3 rd edition, 2003 Contents Preface Acknowledgements Notation Two and three level models. A general classification notation and diagram Glossary Chapter 1 An introduction
More information6 Pattern Mixture Models
6 Pattern Mixture Models A common theme underlying the methods we have discussed so far is that interest focuses on making inference on parameters in a parametric or semiparametric model for the full data
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 informationPart 8: GLMs and Hierarchical LMs and GLMs
Part 8: GLMs and Hierarchical LMs and GLMs 1 Example: Song sparrow reproductive success Arcese et al., (1992) provide data on a sample from a population of 52 female song sparrows studied over the course
More informationSTAT 6350 Analysis of Lifetime Data. Probability Plotting
STAT 6350 Analysis of Lifetime Data Probability Plotting Purpose of Probability Plots Probability plots are an important tool for analyzing data and have been particular popular in the analysis of life
More informationOutline. Mixed models in R using the lme4 package Part 5: Generalized linear mixed models. Parts of LMMs carried over to GLMMs
Outline Mixed models in R using the lme4 package Part 5: Generalized linear mixed models Douglas Bates University of Wisconsin - Madison and R Development Core Team UseR!2009,
More informationPairwise rank based likelihood for estimating the relationship between two homogeneous populations and their mixture proportion
Pairwise rank based likelihood for estimating the relationship between two homogeneous populations and their mixture proportion Glenn Heller and Jing Qin Department of Epidemiology and Biostatistics Memorial
More informationUNIVERSITY OF CALIFORNIA, SAN DIEGO
UNIVERSITY OF CALIFORNIA, SAN DIEGO Estimation of the primary hazard ratio in the presence of a secondary covariate with non-proportional hazards An undergraduate honors thesis submitted to the Department
More informationLatent Variable Models for Binary Data. Suppose that for a given vector of explanatory variables x, the latent
Latent Variable Models for Binary Data Suppose that for a given vector of explanatory variables x, the latent variable, U, has a continuous cumulative distribution function F (u; x) and that the binary
More informationPackage Rsurrogate. October 20, 2016
Type Package Package Rsurrogate October 20, 2016 Title Robust Estimation of the Proportion of Treatment Effect Explained by Surrogate Marker Information Version 2.0 Date 2016-10-19 Author Layla Parast
More informationSurvival Analysis Math 434 Fall 2011
Survival Analysis Math 434 Fall 2011 Part IV: Chap. 8,9.2,9.3,11: Semiparametric Proportional Hazards Regression Jimin Ding Math Dept. www.math.wustl.edu/ jmding/math434/fall09/index.html Basic Model Setup
More informationFULL LIKELIHOOD INFERENCES IN THE COX MODEL
October 20, 2007 FULL LIKELIHOOD INFERENCES IN THE COX MODEL BY JIAN-JIAN REN 1 AND MAI ZHOU 2 University of Central Florida and University of Kentucky Abstract We use the empirical likelihood approach
More informationReview Article Analysis of Longitudinal and Survival Data: Joint Modeling, Inference Methods, and Issues
Journal of Probability and Statistics Volume 2012, Article ID 640153, 17 pages doi:10.1155/2012/640153 Review Article Analysis of Longitudinal and Survival Data: Joint Modeling, Inference Methods, and
More informationEstimation of Optimal Treatment Regimes Via Machine Learning. Marie Davidian
Estimation of Optimal Treatment Regimes Via Machine Learning Marie Davidian Department of Statistics North Carolina State University Triangle Machine Learning Day April 3, 2018 1/28 Optimal DTRs Via ML
More informationGeneralized Linear. Mixed Models. Methods and Applications. Modern Concepts, Walter W. Stroup. Texts in Statistical Science.
Texts in Statistical Science Generalized Linear Mixed Models Modern Concepts, Methods and Applications Walter W. Stroup CRC Press Taylor & Francis Croup Boca Raton London New York CRC Press is an imprint
More informationJoint Modeling of Survival and Longitudinal Data: Likelihood Approach Revisited
Biometrics 62, 1037 1043 December 2006 DOI: 10.1111/j.1541-0420.2006.00570.x Joint Modeling of Survival and Longitudinal Data: Likelihood Approach Revisited Fushing Hsieh, 1 Yi-Kuan Tseng, 2 and Jane-Ling
More informationConstrained Maximum Likelihood Estimation for Model Calibration Using Summary-level Information from External Big Data Sources
Constrained Maximum Likelihood Estimation for Model Calibration Using Summary-level Information from External Big Data Sources Yi-Hau Chen Institute of Statistical Science, Academia Sinica Joint with Nilanjan
More informationIntroduction to Empirical Processes and Semiparametric Inference Lecture 01: Introduction and Overview
Introduction to Empirical Processes and Semiparametric Inference Lecture 01: Introduction and Overview Michael R. Kosorok, Ph.D. Professor and Chair of Biostatistics Professor of Statistics and Operations
More informationConstrained estimation for binary and survival data
Constrained estimation for binary and survival data Jeremy M. G. Taylor Yong Seok Park John D. Kalbfleisch Biostatistics, University of Michigan May, 2010 () Constrained estimation May, 2010 1 / 43 Outline
More informationMixed models in R using the lme4 package Part 5: Generalized linear mixed models
Mixed models in R using the lme4 package Part 5: Generalized linear mixed models Douglas Bates Madison January 11, 2011 Contents 1 Definition 1 2 Links 2 3 Example 7 4 Model building 9 5 Conclusions 14
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 informationCausal inference in epidemiological practice
Causal inference in epidemiological practice Willem van der Wal Biostatistics, Julius Center UMC Utrecht June 5, 2 Overview Introduction to causal inference Marginal causal effects Estimating marginal
More informationRonald Christensen. University of New Mexico. Albuquerque, New Mexico. Wesley Johnson. University of California, Irvine. Irvine, California
Texts in Statistical Science Bayesian Ideas and Data Analysis An Introduction for Scientists and Statisticians Ronald Christensen University of New Mexico Albuquerque, New Mexico Wesley Johnson University
More informationFrailty Modeling for clustered survival data: a simulation study
Frailty Modeling for clustered survival data: a simulation study IAA Oslo 2015 Souad ROMDHANE LaREMFiQ - IHEC University of Sousse (Tunisia) souad_romdhane@yahoo.fr Lotfi BELKACEM LaREMFiQ - IHEC University
More informationBayesian Multivariate Logistic Regression
Bayesian Multivariate Logistic Regression Sean M. O Brien and David B. Dunson Biostatistics Branch National Institute of Environmental Health Sciences Research Triangle Park, NC 1 Goals Brief review of
More informationSemiparametric Generalized Linear Models
Semiparametric Generalized Linear Models North American Stata Users Group Meeting Chicago, Illinois Paul Rathouz Department of Health Studies University of Chicago prathouz@uchicago.edu Liping Gao MS Student
More informationOn Estimating the Relationship between Longitudinal Measurements and Time-to-Event Data Using a Simple Two-Stage Procedure
Biometrics DOI: 10.1111/j.1541-0420.2009.01324.x On Estimating the Relationship between Longitudinal Measurements and Time-to-Event Data Using a Simple Two-Stage Procedure Paul S. Albert 1, and Joanna
More informationPQL Estimation Biases in Generalized Linear Mixed Models
PQL Estimation Biases in Generalized Linear Mixed Models Woncheol Jang Johan Lim March 18, 2006 Abstract The penalized quasi-likelihood (PQL) approach is the most common estimation procedure for the generalized
More informationSurvival Analysis. Lu Tian and Richard Olshen Stanford University
1 Survival Analysis Lu Tian and Richard Olshen Stanford University 2 Survival Time/ Failure Time/Event Time We will introduce various statistical methods for analyzing survival outcomes What is the survival
More informationGeneralized Linear Mixed-Effects Models. Copyright c 2015 Dan Nettleton (Iowa State University) Statistics / 58
Generalized Linear Mixed-Effects Models Copyright c 2015 Dan Nettleton (Iowa State University) Statistics 510 1 / 58 Reconsideration of the Plant Fungus Example Consider again the experiment designed to
More informationBIOS 2083: Linear Models
BIOS 2083: Linear Models Abdus S Wahed September 2, 2009 Chapter 0 2 Chapter 1 Introduction to linear models 1.1 Linear Models: Definition and Examples Example 1.1.1. Estimating the mean of a N(μ, σ 2
More informationPackage JointModel. R topics documented: June 9, Title Semiparametric Joint Models for Longitudinal and Counting Processes Version 1.
Package JointModel June 9, 2016 Title Semiparametric Joint Models for Longitudinal and Counting Processes Version 1.0 Date 2016-06-01 Author Sehee Kim Maintainer Sehee Kim
More informationType I and type II error under random-effects misspecification in generalized linear mixed models Link Peer-reviewed author version
Type I and type II error under random-effects misspecification in generalized linear mixed models Link Peer-reviewed author version Made available by Hasselt University Library in Document Server@UHasselt
More informationLecture 5 Models and methods for recurrent event data
Lecture 5 Models and methods for recurrent event data Recurrent and multiple events are commonly encountered in longitudinal studies. In this chapter we consider ordered recurrent and multiple events.
More informationAccelerated Failure Time Models
Accelerated Failure Time Models Patrick Breheny October 12 Patrick Breheny University of Iowa Survival Data Analysis (BIOS 7210) 1 / 29 The AFT model framework Last time, we introduced the Weibull distribution
More informationGOODNESS-OF-FIT TEST ISSUES IN GENERALIZED LINEAR MIXED MODELS. A Dissertation NAI-WEI CHEN
GOODNESS-OF-FIT TEST ISSUES IN GENERALIZED LINEAR MIXED MODELS A Dissertation by NAI-WEI CHEN Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the requirements
More informationare censored. Because artificial censoring reduces the information available and leads to nonsmooth estimating functions, prior research has noted
ABSTRACT VOCK, DAVID MICHAEL. Advanced Statistical Methods for Complex Longitudinal Data. (Under the direction of Marie Davidian and Anastasios A. Tsiatis.) Longitudinally collected data play an important
More informationMixed models in R using the lme4 package Part 5: Generalized linear mixed models
Mixed models in R using the lme4 package Part 5: Generalized linear mixed models Douglas Bates 2011-03-16 Contents 1 Generalized Linear Mixed Models Generalized Linear Mixed Models When using linear mixed
More informationAnalysis of Longitudinal Data. Patrick J. Heagerty PhD Department of Biostatistics University of Washington
Analsis of Longitudinal Data Patrick J. Heagert PhD Department of Biostatistics Universit of Washington 1 Auckland 2008 Session Three Outline Role of correlation Impact proper standard errors Used to weight
More informationStatistical Methods for Alzheimer s Disease Studies
Statistical Methods for Alzheimer s Disease Studies Rebecca A. Betensky, Ph.D. Department of Biostatistics, Harvard T.H. Chan School of Public Health July 19, 2016 1/37 OUTLINE 1 Statistical collaborations
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 informationLikelihood Construction, Inference for Parametric Survival Distributions
Week 1 Likelihood Construction, Inference for Parametric Survival Distributions In this section we obtain the likelihood function for noninformatively rightcensored survival data and indicate how to make
More informationLinear Regression Models P8111
Linear Regression Models P8111 Lecture 25 Jeff Goldsmith April 26, 2016 1 of 37 Today s Lecture Logistic regression / GLMs Model framework Interpretation Estimation 2 of 37 Linear regression Course started
More informationT E C H N I C A L R E P O R T A SANDWICH-ESTIMATOR TEST FOR MISSPECIFICATION IN MIXED-EFFECTS MODELS. LITIERE S., ALONSO A., and G.
T E C H N I C A L R E P O R T 0658 A SANDWICH-ESTIMATOR TEST FOR MISSPECIFICATION IN MIXED-EFFECTS MODELS LITIERE S., ALONSO A., and G. MOLENBERGHS * I A P S T A T I S T I C S N E T W O R K INTERUNIVERSITY
More informationSemiparametric Regression
Semiparametric Regression Patrick Breheny October 22 Patrick Breheny Survival Data Analysis (BIOS 7210) 1/23 Introduction Over the past few weeks, we ve introduced a variety of regression models under
More informationCTDL-Positive Stable Frailty Model
CTDL-Positive Stable Frailty Model M. Blagojevic 1, G. MacKenzie 2 1 Department of Mathematics, Keele University, Staffordshire ST5 5BG,UK and 2 Centre of Biostatistics, University of Limerick, Ireland
More information[Part 2] Model Development for the Prediction of Survival Times using Longitudinal Measurements
[Part 2] Model Development for the Prediction of Survival Times using Longitudinal Measurements Aasthaa Bansal PhD Pharmaceutical Outcomes Research & Policy Program University of Washington 69 Biomarkers
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 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 informationAnalysis of competing risks data and simulation of data following predened subdistribution hazards
Analysis of competing risks data and simulation of data following predened subdistribution hazards Bernhard Haller Institut für Medizinische Statistik und Epidemiologie Technische Universität München 27.05.2013
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 informationModelling Survival Events with Longitudinal Data Measured with Error
Modelling Survival Events with Longitudinal Data Measured with Error Hongsheng Dai, Jianxin Pan & Yanchun Bao First version: 14 December 29 Research Report No. 16, 29, Probability and Statistics Group
More informationNovember 2002 STA Random Effects Selection in Linear Mixed Models
November 2002 STA216 1 Random Effects Selection in Linear Mixed Models November 2002 STA216 2 Introduction It is common practice in many applications to collect multiple measurements on a subject. Linear
More informationReview. Timothy Hanson. Department of Statistics, University of South Carolina. Stat 770: Categorical Data Analysis
Review Timothy Hanson Department of Statistics, University of South Carolina Stat 770: Categorical Data Analysis 1 / 22 Chapter 1: background Nominal, ordinal, interval data. Distributions: Poisson, binomial,
More informationSTAT 705 Generalized linear mixed models
STAT 705 Generalized linear mixed models Timothy Hanson Department of Statistics, University of South Carolina Stat 705: Data Analysis II 1 / 24 Generalized Linear Mixed Models We have considered random
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 information1 The problem of survival analysis
1 The problem of survival analysis Survival analysis concerns analyzing the time to the occurrence of an event. For instance, we have a dataset in which the times are 1, 5, 9, 20, and 22. Perhaps those
More informationLecture 14: Introduction to Poisson Regression
Lecture 14: Introduction to Poisson Regression Ani Manichaikul amanicha@jhsph.edu 8 May 2007 1 / 52 Overview Modelling counts Contingency tables Poisson regression models 2 / 52 Modelling counts I Why
More informationModelling counts. Lecture 14: Introduction to Poisson Regression. Overview
Modelling counts I Lecture 14: Introduction to Poisson Regression Ani Manichaikul amanicha@jhsph.edu Why count data? Number of traffic accidents per day Mortality counts in a given neighborhood, per week
More information