Approaches to Modeling Menstrual Cycle Function
|
|
- Ralph Simpson
- 6 years ago
- Views:
Transcription
1 Approaches to Modeling Menstrual Cycle Function Paul S. Albert Biostatistics & Bioinformatics Branch Division of Epidemiology, Statistics, and Prevention Research NICHD SPER Student Workshop June 22,
2 Menstrual Cycle Function 2
3 Scientific Interest To characterize the typical menstrual cycle pattern in a population of women. What is the effect of a subject-specific covariate on typical menstrual cycle? To characterize the variation in menstrual cycle function across women and across consecutive cycles on the same woman. What is the inter-relationship between multiple hormones across the menstrual cycle? To characterize the relationship between menstrual cycle function and menstrual cycle length. 3
4 Longitudinal Menstrual Cycle Pattern Menstrual Early Follicular Ovulation Early Luteal Mid Luteal Late Luteal Urine measurements taken daily -or- Serum measurements taken at regularly scheduled time points. 4
5 Time Scale for Menstrual Cycle Function 1. Assign visit days to phases in the cycle 2. Actual visit days 3. Actual visit days standardized by cycle length 4. Registered cycles: Standardize by day of ovulation for time period before ovulation. Standardize by time between ovulation and menstruation for time period between ovulation and menstruation. 5
6 Linear Mixed Model Pinheiro and Bates-Springer, 2000 For follow-up data at fixed scheduled time points: i b Z X Y ε β + + = Fixed Effects Random Effects = = = 8 2 ) ( 1 0 d i d d j i b I X Y ε β β β = = = = = ) ( ) ( 1 0 d i d i d j d d j d i b X I I X Y ε η β β β 6
7 Y = β + β X 0 1 Linear Mixed Model Pinheiro and Bates-Springer, 2000 i β d I( j= d ) + ηd I( j= d ) X i + bi + d = 2 d = 2 Test of whether a covariate influences the menstrual cycle pattern of a particular hormone can be formulated as a likelihood ratio test of η = η = = η 3 8 ε 7
8 Advantages: Linear Mixed Model Pinheiro and Bates-Springer, 2000 Simple to implement in standard software Easily allows for flexible correlation structure (example: AR-1, ARIMA, exponential model, spherical model) Model diagnostics built in in standard software such as lme. 8
9 Advantages (Continued): Y Can be easily extended to multiple cycles per women = β + β X 0 + Disadvantages: 1 Linear Mixed Model Pinheiro and Bates-Springer, 2000 i 8 8 β d I( j= d ) + ηd I( j= d ) X i + bi + bik + d = 2 d = 2 Requires uniform measurement times on each subject. Global tests are not very powerful. Difficult to study effect of covariates on hormonal pattern. ε 9
10 Semi-parametric Stochastic Mixed Model Zhang, Lin, Raz, and Sowers - JASA 1998 Y = X ' β + f ( t ) + Z' b + U ( t ) + ε i i Fixed Effects Smoothed function of time Random Effects Random process Random process follows a nonhomogenous Ornstein-Uhlenbeck (NOU) process 10
11 Semi-parametric Stochastic Mixed Model Zhang, Lin, Raz, and Sowers - JASA 1998 NOU process: Var( Ui ( t)) = ζ ( t) with log{ ζ ( t)} = ζ + ζ t 1 0 Corr( U i ( t), U i ( s)) = ρ t s 11
12 Semi-parametric Stochastic Mixed Model Zhang, Lin, Raz, and Sowers - JASA 1998 Advantages: Allows for irregularly spaced measurements. Accounts for serial correlation and variances of measurements that change over the menstrual cycle. Allows for very flexible hormonal patterns. Disadvantages: Covariates affect cycle only by a constant shift in mean. Needs specialized software to implement. Does not easily extend to multiple cycles per women. 12
13 Semi-parametric Stochastic Mixed Model Zhang, Lin, and Sowers - Biometrics 2000 Comparing the non-parametric curves across two groups Y k = X ' k β + f ( t ) + Z' b + U ( t ) + ε k ki ki k k Groups K=1 or 2 Are the two curves different across groups? { T f = 1 (.), f2(.)} [ f1( t) f2( t)] 0 Chi-Squared test to compare the two curves 2 dt 13
14 Semi-parametric Stochastic Mixed Model Zhang, Lin, and Sowers - Biometrics 2000 Advantages: Very flexible and general test Disadvantages: Only compares two discrete groups Difficult to implement with standard software 14
15 Characterizing Hormonal Profiles Using Functional Data Analysis Meyer, Zeger, Harlow, et al. - AJE 2007 Fit separate curves to each subject s menstrual cycle data. Cubic B-splines with 20 basis functions were used with associated weights estimated from the data. Mean curves for each group were estimated by averaging estimated weights and multiplying average weights by each basis function. 15
16 Characterizing Hormonal Profiles Using Functional Data Analysis Meyer, Zeger, Harlow, et al. - AJE 2007 Interest is on comparing menstrual cycle function by cycle length 16
17 Characterizing Hormonal Profiles Using Functional Data Analysis Meyer, Zeger, Harlow, et al. - AJE 2007 Evaluating synchrony between FSH and LH: 17
18 Characterizing Hormonal Profiles Using Functional Data Analysis Meyer, Zeger, Harlow, et al. AJE 2007 Advantages: Summarizes menstrual cycle function in a sensible way. Easily allows for examining how different hormones relate to each other. Disadvantages: Needs many measurements on each cycle. Good for urine. May not work for serum measurements. Does not easily adapt to a regression framework. Focus is not on statistical testing, but is more descriptive. 18
19 y Nonlinear Mixed Models: Shape Invariant Model Wang, Ye, and Brown - Biometrics 2003 Albert and Hunsberger - Biometrics 2005 = φ + φ f t a it φ )} + ε 1i φ 1i = X 1i β + Z 1i b 1i φ 2i = X 2i β + Z 2i b 2i φ 3i = X 3i β + Z 3i b 3i exp( 2i ) { log ( 3i Mean level Amplitude: nadir to peak distance Phase shift a logit( x) = exp( x) /{1 + exp( x)} f(t) is a periodic function with t ranging from 0 to 1 19
20 Nonlinear Mixed Models: Shape Invariant Model Wang, Ye, and Brown- Biometrics 2003 Albert and Hunsberger - Biometrics 2005 y = φ + exp( φ ) f { t a logit( φ3 )} + ε 1i 2i i Mean level Mean Shift hormone value scaled time 20
21 Nonlinear Mixed Models: Shape Invariant Model Wang, Ye, and Brown- Biometrics 2003 Albert and Hunsberger - Biometrics 2005 y = φ + exp( φ ) f { t a logit( φ3 )} + ε 1i 2i i Amplitude: nadir to peak distance Amplitude Shift hormone value scaled time 21
22 Nonlinear Mixed Models: Shape Invariant Model Wang, Ye, and Brown- Biometrics 2003 Albert and Hunsberger - Biometrics 2005 y = φ + exp( φ ) f { t a logit( φ3 )} + ε 1i 2i i Phase Shift Phase shift hormone value scaled time 22
23 Nonlinear Mixed Models: Shape Invariant Model Wang, Ye, and Brown- Biometrics 2003 Albert and Hunsberger - Biometrics 2005 y = φ + φ f t a it φ )} + ε 1i exp( 2i ) { log ( 3i f is a periodic function which can be represented in a number of ways: Periodic splines (Wang, Ke, and Brown 2003) Harmonic models (Albert and Hunsberger 2005) 23
24 Nonlinear Mixed Models: Shape Invariant Model Characterizing f Using periodic splines: f and f are absolutely continuous f(0)=f(1), f (0)=f (1) 1 0 f ( t) dt = Using harmonic models: 0 f K ( t) = β cos(2πkt + θ ) k = 1 k k 24
25 Nonlinear Mixed Models: Shape Invariant Model Characterizing f with harmonic models K=1 (one harmonic) hormone fraction of standardized menstrual cycle 25
26 Nonlinear Mixed Models: Shape Invariant Model Characterizing f with harmonic models More than one harmonic 26
27 Nonlinear Mixed Models: Shape Invariant Model Wang, Ye, and Brown- Biometrics 2003 Albert and Hunsberger - Biometrics 2005 y = φ + φ f t a it φ )} + ε 1i exp( 2i ) { log ( 3i Mean level Amplitude: nadir to peak distance Phase shift The model allows for complex between subject variation in mean, amplitude, and phase shift 27
28 Complex Between-Subject Variation Variation in Amplitude & Phase Not Mean Variation in Mean & Phase Not Amplitude Variation in Mean & Amplitude Not Phase 28
29 How to Choose the Number of Harmonics (K)? Akaike s Information Criterion (AIC) and the Bayesian Information Criterion (BIC) for choosing K. AIC( p) = 2log L + 2 p BIC( p) = 2log L p log( n i= 1 Since f is very flexible with 4 harmonics, we recommend choosing 1,2, 3, or 4 harmonics. Inference on covariates not very sensitive to choice of number of harmonics. + n i ) 29
30 Testing the Shape Invariant Assumption Simple model for two groups (G=0 or 1): yg =τ G + b1 i + exp( b2i + d1g) f {( t a logit( b3i + d2g))} + ε Likelihood ratio test for shape invariance: f K ( t) = β1 k cos(2kπ ( t + θ1 k a logit( φ3i )) k = 1 + β G cos( 2k( t + θ2kg a logit( φ3 ))) 2k i k Test β 2 k = θ2 k = 0 for k = 2,3,..., K 30
31 Nonlinear Mixed Models: Shape Invariant Model Non-parametric versus harmonic modeling of f Non-parametric models must choose the smoothing parameter while harmonic models must choose the number of harmonics. Non-parametric models require specialized software. Wang et al. developed a package called ASSIST for SPLUS. Harmonic models are easy to fit in R or SAS with standard code. 31
32 Estimation in R using nlme Uses an estimation algorithm described by Lindstrom and Bates (1990) Sample R code: 32
33 Estimation in R using nlme: Choosing starting values Choose starting values for fixed effects by fitting a non-linear model without random effects. Sensitivity to starting values is more pronounced when the number of harmonics is large 33
34 Statistical Properties of the Procedure Even in small samples, estimation is unbiased and tests of coefficients have the correct type- I error (with and without estimating the number of harmonics). 34
35 Using the Shape Invariant Model for Menstrual Cycle Longitudinal Data y k = φ + )} + ε 1ik exp( φ2ik ) f { tk / Tik a logit( φ3ik k th cycle j th Standardized cycle length time point k i th subject a logit( x) = exp( x) /{1 + exp( x)} 35
36 Using the Shape Invariant Model for Menstrual Cycle Longitudinal Data Random Effects: φ φ φ 1 ik X1 i + b1 i + b1 ik = β 2 ik = X 2iβ + b2i + b2ik 3 ik = X 3iβ + b3i + b3ik ( 1i, b2i, b3i ) is b N 0, Σ ) ( b Between-subject variation ( b1 ik, b2ik, b3ik ) is N ( c 0, Σ ) Within-subject variation 36
37 A Comparison of Methods Method Advantages Disadvantages Linear Mixed Model Semi-parametric Stochastic Mixed Model Non-linear Mixed Model: Shape Invariant Model - Simple to implement - Built-in diagnostics - Easy to incorporate multiple cycles - Irregularly timed measurements - Flexible hormonal pattern - Irregularly timed measurements - Flexible way to examine covariates on hormonal process - Uniform measurements - Studying effect of covariates on hormonal process is limited - Studying effect of covariates on hormonal process is limited - Needs specialized software - Shape invariance assumption 37
38 Advanced Topics and New Methods 38
39 Menstrual Cycle Pattern May Vary By Cycle Length 39
40 Joint modeling of longitudinal menstrual cycle and menstrual cycle length data Measure hormones at multiple fixed time points over a 28 day period. Actual cycle length has a much wider range: days in BioCycle. Actual pattern may vary by menstrual cycle length y = φ + exp( φ ) f { t / T a logit( φ3 )} + ε 1i 2i i i T i is N ( + µ φ + µ φ + µ φ 2 µ 0 1 1i 2 2i 3 3i σ, ) 40
41 Joint modeling of longitudinal menstrual cycle and menstrual cycle length data In BioCycle, women enrolled with regular menstrual cycle length Modeling menstrual cycle length using a mixture distribution (Guo et al. Biostatistics 2002) 41
42 Relationship Between Menstrual Cycle Length and Number of Measurements Informative missing data problem 42
43 Estimation of Joint Model y T i = φ + exp( φ ) f { t / T a logit( φ3 )} + ε is 1i 2i N ( + µ φ + µ φ + µ φ 2 µ 0 1 1i 2 2i 3 3i σ i, i ) L = I J b i= 1 j= 1 i f ( y T i, b) f b) g( b) db Use trapezoidal rule for numerical integration Obtain maximum-likelihood estimates ( T i 43
44 Results of Analyses for Progesterone 3 harmonics model: Joint model and naïve model not considering menstrual cycle length Joint model Naïve model Parameters Estimate SE Estimate SE μ μ μ μ β β β σ b σ b σ b
45 Results for Progesterone 45
46 Simulation Studies Under the correct joint model - 3 harmonics, 50 individuals, 500 simulated datasets, and 5 time points: 2, 7, 13, 18, and 27. Parameters Truth Estimate MC SE Aympt SE μ μ μ μ β β β
47 Simulation Studies Naïve model - 3 harmonics, 50 individuals, 500 simulated datasets, and 5 time points: 2, 7, 13, 18, and 27. Parameters Truth Estimate MC SE Aympt SE β β β
48 Simulation Studies Naïve model -3 harmonics, 50 individuals, 1000 simulated datasets, and 5 time points: 2, 7, 13, 18, and 27 : µ 1 = µ 2 = 10 Parameters Truth Estimate MC SE Aympt SE β β β
49 Extensions to Modeling Multivariate Hormonal Profiles m th hormone y mk = φ + )} + ε 1imk exp( φ2imk ) f { tk / Tik a logit( φ3imk mk a logit( x) = exp( x) /{1 + exp( x)} φ φ φ 1 imk X1 i m + b1 im + b1 ikm 2 imk = X 2imβm + b2im + b2ikm = β 3 imk = X 3imβm + b3im + b3ikm 49
50 Extensions to Modeling Multivariate Hormonal Profiles Let b = im ( b im, b2im, b3 1 im )' and b = ikm ( b ikm, b2ikm, b3 1 ikm )' Assume that Cov(( b, b 2,..., b i1 i im )') = Σ b Between-Subject Variation Cov(( b, b 2,..., b ik1 ik ikm )') = Σ W Within-Subject Variation 50
51 Acknowledgements This work was supported by the Intramural Research Program of the Eunice Kennedy Shriver National Institute of Child Health & Human Development, National Institutes of Health. BioCycle working group: Sunni Mumford, Anna Pollack, Enrique Schisterman, Aun Ye, Jean Wactawski-Wende, Edwina Yeung, Cuilin Zhang
52 References 1. Albert, P.S. and Hunsberger, S. (2005). On analyzing circadian rhythms data using nonlinear mixed models with harmonic terms. Biometrics 61, Guo, Y., Manatunga, A.K., Chen, S., and Marcus, M. (2006). Modeling menstrual cycle length using a mixture model. Biostatistics 7, Meyer, P.M., Zeger, S.L., et al. (2007). Characterizing daily urinary hormone profiles for women at midlife using functional data analysis. American Journal of Epidemiology 165, Pinheiro, J.C. and Bates, D.M. Mixed-Effects Models in S and S-PLUS. Springer, New-York ISBN Wang Y., Ke C., and Brown, M.B. (2003). Shape-invariant modeling of circadian rhythms with random effects and smoothing spline ANOVA decompositions. Biometrics 59, Zhang, D., Raz, J., and Sowers, M. (1998). Semiparametric stochastic mixed models for longitudinal data. Journal of the American Statistical Association 93, Zhang, D. Lin, X., and Sowers, M. (2000). Semiparametric regression for periodic longitudinal hormone data with multiple menstrual cycles. Biometrics 56,
A Bayesian multi-dimensional couple-based latent risk model for infertility
A Bayesian multi-dimensional couple-based latent risk model for infertility Zhen Chen, Ph.D. Eunice Kennedy Shriver National Institute of Child Health and Human Development National Institutes of Health
More informationOnline Supplementary Material to: Perfluoroalkyl Chemicals, Menstrual Cycle Length and Fecundity: Findings from a. Prospective Pregnancy Study
Online Supplementary Material to: Perfluoroalkyl Chemicals, Menstrual Cycle Length and Fecundity: Findings from a Prospective Pregnancy Study by Kirsten J. Lum 1,a,b, Rajeshwari Sundaram a, Dana Boyd Barr
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 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 informationLinear, Generalized Linear, and Mixed-Effects Models in R. Linear and Generalized Linear Models in R Topics
Linear, Generalized Linear, and Mixed-Effects Models in R John Fox McMaster University ICPSR 2018 John Fox (McMaster University) Statistical Models in R ICPSR 2018 1 / 19 Linear and Generalized Linear
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 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 informationThe forecasting of menstruation based on a state-space modeling of basal body temperature time series
Research Article Received 8 July 2016, Accepted 27 April 2017 Published online 22 May 2017 in Wiley Online Library (wileyonlinelibrary.com) DOI: 10.1002/sim.7345 The forecasting of menstruation based on
More informationFractional Polynomials and Model Averaging
Fractional Polynomials and Model Averaging Paul C Lambert Center for Biostatistics and Genetic Epidemiology University of Leicester UK paul.lambert@le.ac.uk Nordic and Baltic Stata Users Group Meeting,
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 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 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 informationBayesian inference in semiparametric mixed models for longitudinal data
Bayesian inference in semiparametric mixed models for longitudinal data Yisheng Li 1, Xihong Lin 2 and Peter Müller 1 1 Department of Biostatistics, Division of Quantitative Sciences University of Texas
More informationTIME SERIES ANALYSIS AND FORECASTING USING THE STATISTICAL MODEL ARIMA
CHAPTER 6 TIME SERIES ANALYSIS AND FORECASTING USING THE STATISTICAL MODEL ARIMA 6.1. Introduction A time series is a sequence of observations ordered in time. A basic assumption in the time series analysis
More informationThe equivalence of the Maximum Likelihood and a modified Least Squares for a case of Generalized Linear Model
Applied and Computational Mathematics 2014; 3(5): 268-272 Published online November 10, 2014 (http://www.sciencepublishinggroup.com/j/acm) doi: 10.11648/j.acm.20140305.22 ISSN: 2328-5605 (Print); ISSN:
More informationRejoinder. 1 Phase I and Phase II Profile Monitoring. Peihua Qiu 1, Changliang Zou 2 and Zhaojun Wang 2
Rejoinder Peihua Qiu 1, Changliang Zou 2 and Zhaojun Wang 2 1 School of Statistics, University of Minnesota 2 LPMC and Department of Statistics, Nankai University, China We thank the editor Professor David
More informationSTAT 518 Intro Student Presentation
STAT 518 Intro Student Presentation Wen Wei Loh April 11, 2013 Title of paper Radford M. Neal [1999] Bayesian Statistics, 6: 475-501, 1999 What the paper is about Regression and Classification Flexible
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 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 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 informationPrerequisite: STATS 7 or STATS 8 or AP90 or (STATS 120A and STATS 120B and STATS 120C). AP90 with a minimum score of 3
University of California, Irvine 2017-2018 1 Statistics (STATS) Courses STATS 5. Seminar in Data Science. 1 Unit. An introduction to the field of Data Science; intended for entering freshman and transfers.
More 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 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 informationHST 583 FUNCTIONAL MAGNETIC RESONANCE IMAGING DATA ANALYSIS AND ACQUISITION A REVIEW OF STATISTICS FOR FMRI DATA ANALYSIS
HST 583 FUNCTIONAL MAGNETIC RESONANCE IMAGING DATA ANALYSIS AND ACQUISITION A REVIEW OF STATISTICS FOR FMRI DATA ANALYSIS EMERY N. BROWN AND CHRIS LONG NEUROSCIENCE STATISTICS RESEARCH LABORATORY DEPARTMENT
More informationBayesian inference on dependence in multivariate longitudinal data
Biometrics 000, 000 000 DOI: 000 000 0000 Bayesian inference on dependence in multivariate longitudinal data Hongxia Yang a, Fan Li a, Enrique Schisterman b and David Dunson a a Department of Statistical
More informationSemiparametric Mixed Effects Models with Flexible Random Effects Distribution
Semiparametric Mixed Effects Models with Flexible Random Effects Distribution Marie Davidian North Carolina State University davidian@stat.ncsu.edu www.stat.ncsu.edu/ davidian Joint work with A. Tsiatis,
More informationContrasting Models for Lactation Curve Analysis
J. Dairy Sci. 85:968 975 American Dairy Science Association, 2002. Contrasting Models for Lactation Curve Analysis F. Jaffrezic,*, I. M. S. White,* R. Thompson, and P. M. Visscher* *Institute of Cell,
More informationECONOMETRICS HONOR S EXAM REVIEW SESSION
ECONOMETRICS HONOR S EXAM REVIEW SESSION Eunice Han ehan@fas.harvard.edu March 26 th, 2013 Harvard University Information 2 Exam: April 3 rd 3-6pm @ Emerson 105 Bring a calculator and extra pens. Notes
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 Dynamic Linear Modelling for. Complex Computer Models
Bayesian Dynamic Linear Modelling for Complex Computer Models Fei Liu, Liang Zhang, Mike West Abstract Computer models may have functional outputs. With no loss of generality, we assume that a single computer
More informationA joint modeling approach for multivariate survival data with random length
A joint modeling approach for multivariate survival data with random length Shuling Liu, Emory University Amita Manatunga, Emory University Limin Peng, Emory University Michele Marcus, Emory 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 informationEconometrics Honor s Exam Review Session. Spring 2012 Eunice Han
Econometrics Honor s Exam Review Session Spring 2012 Eunice Han Topics 1. OLS The Assumptions Omitted Variable Bias Conditional Mean Independence Hypothesis Testing and Confidence Intervals Homoskedasticity
More informationHypothesis Testing in Smoothing Spline Models
Hypothesis Testing in Smoothing Spline Models Anna Liu and Yuedong Wang October 10, 2002 Abstract This article provides a unified and comparative review of some existing test methods for the hypothesis
More informationBIOSTATISTICAL METHODS
BIOSTATISTICAL METHODS FOR TRANSLATIONAL & CLINICAL RESEARCH Cross-over Designs #: DESIGNING CLINICAL RESEARCH The subtraction of measurements from the same subject will mostly cancel or minimize effects
More informationModeling Real Estate Data using Quantile Regression
Modeling Real Estate Data using Semiparametric Quantile Regression Department of Statistics University of Innsbruck September 9th, 2011 Overview 1 Application: 2 3 4 Hedonic regression data for house prices
More informationImpact of serial correlation structures on random effect misspecification with the linear mixed model.
Impact of serial correlation structures on random effect misspecification with the linear mixed model. Brandon LeBeau University of Iowa file:///c:/users/bleb/onedrive%20 %20University%20of%20Iowa%201/JournalArticlesInProgress/Diss/Study2/Pres/pres.html#(2)
More informationCourse topics (tentative) The role of random effects
Course topics (tentative) random effects linear mixed models analysis of variance frequentist likelihood-based inference (MLE and REML) prediction Bayesian inference The role of random effects Rasmus Waagepetersen
More informationLikelihood ratio testing for zero variance components in linear mixed models
Likelihood ratio testing for zero variance components in linear mixed models Sonja Greven 1,3, Ciprian Crainiceanu 2, Annette Peters 3 and Helmut Küchenhoff 1 1 Department of Statistics, LMU Munich University,
More informationSTK4900/ Lecture 10. Program
STK4900/9900 - Lecture 10 Program 1. Repeated measures and longitudinal data 2. Simple analysis approaches 3. Random effects models 4. Generalized estimating equations (GEE) 5. GEE for binary data (and
More informationMixed effects models
Mixed effects models The basic theory and application in R Mitchel van Loon Research Paper Business Analytics Mixed effects models The basic theory and application in R Author: Mitchel van Loon Research
More informationOn the Behavior of Marginal and Conditional Akaike Information Criteria in Linear Mixed Models
On the Behavior of Marginal and Conditional Akaike Information Criteria in Linear Mixed Models Thomas Kneib Department of Mathematics Carl von Ossietzky University Oldenburg Sonja Greven Department of
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 informationA general mixed model approach for spatio-temporal regression data
A general mixed model approach for spatio-temporal regression data Thomas Kneib, Ludwig Fahrmeir & Stefan Lang Department of Statistics, Ludwig-Maximilians-University Munich 1. Spatio-temporal regression
More informationLASR 2015 Geometry-Driven Statistics and its Cutting Edge Applications: Celebrating Four Decades of Leeds Statistics Workshops
LASR 2015 Geometry-Driven Statistics and its Cutting Edge Applications: Celebrating Four Decades of Leeds Statistics Workshops Programme and Abstracts Edited by K.V. Mardia, A. Gusnanto, C. Nooney & J.
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 informationCircle a single answer for each multiple choice question. Your choice should be made clearly.
TEST #1 STA 4853 March 4, 215 Name: Please read the following directions. DO NOT TURN THE PAGE UNTIL INSTRUCTED TO DO SO Directions This exam is closed book and closed notes. There are 31 questions. Circle
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 informationCharles E. McCulloch Biometrics Unit and Statistics Center Cornell University
A SURVEY OF VARIANCE COMPONENTS ESTIMATION FROM BINARY DATA by Charles E. McCulloch Biometrics Unit and Statistics Center Cornell University BU-1211-M May 1993 ABSTRACT The basic problem of variance components
More informationOn the Behavior of Marginal and Conditional Akaike Information Criteria in Linear Mixed Models
On the Behavior of Marginal and Conditional Akaike Information Criteria in Linear Mixed Models Thomas Kneib Institute of Statistics and Econometrics Georg-August-University Göttingen Department of Statistics
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 informationRegression of Time Series
Mahlerʼs Guide to Regression of Time Series CAS Exam S prepared by Howard C. Mahler, FCAS Copyright 2016 by Howard C. Mahler. Study Aid 2016F-S-9Supplement Howard Mahler hmahler@mac.com www.howardmahler.com/teaching
More informationSTAT Financial Time Series
STAT 6104 - Financial Time Series Chapter 4 - Estimation in the time Domain Chun Yip Yau (CUHK) STAT 6104:Financial Time Series 1 / 46 Agenda 1 Introduction 2 Moment Estimates 3 Autoregressive Models (AR
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 informationLesson 6 Maximum Likelihood: Part III. (plus a quick bestiary of models)
Lesson 6 Maximum Likelihood: Part III (plus a quick bestiary of models) Homework You want to know the density of fish in a set of experimental ponds You observe the following counts in ten ponds: 5,6,7,3,6,5,8,4,4,3
More informationLikelihood Ratio Tests. that Certain Variance Components Are Zero. Ciprian M. Crainiceanu. Department of Statistical Science
1 Likelihood Ratio Tests that Certain Variance Components Are Zero Ciprian M. Crainiceanu Department of Statistical Science www.people.cornell.edu/pages/cmc59 Work done jointly with David Ruppert, School
More informationHigh-dimensional regression modeling
High-dimensional regression modeling David Causeur Department of Statistics and Computer Science Agrocampus Ouest IRMAR CNRS UMR 6625 http://www.agrocampus-ouest.fr/math/causeur/ Course objectives Making
More informationLISA Short Course Series Generalized Linear Models (GLMs) & Categorical Data Analysis (CDA) in R. Liang (Sally) Shan Nov. 4, 2014
LISA Short Course Series Generalized Linear Models (GLMs) & Categorical Data Analysis (CDA) in R Liang (Sally) Shan Nov. 4, 2014 L Laboratory for Interdisciplinary Statistical Analysis LISA helps VT researchers
More informationGroup Sequential Tests for Delayed Responses. Christopher Jennison. Lisa Hampson. Workshop on Special Topics on Sequential Methodology
Group Sequential Tests for Delayed Responses Christopher Jennison Department of Mathematical Sciences, University of Bath, UK http://people.bath.ac.uk/mascj Lisa Hampson Department of Mathematics and Statistics,
More informationIntroduction to Within-Person Analysis and RM ANOVA
Introduction to Within-Person Analysis and RM ANOVA Today s Class: From between-person to within-person ANOVAs for longitudinal data Variance model comparisons using 2 LL CLP 944: Lecture 3 1 The Two Sides
More informationMixture modelling of recurrent event times with long-term survivors: Analysis of Hutterite birth intervals. John W. Mac McDonald & Alessandro Rosina
Mixture modelling of recurrent event times with long-term survivors: Analysis of Hutterite birth intervals John W. Mac McDonald & Alessandro Rosina Quantitative Methods in the Social Sciences Seminar -
More 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 informationarxiv: v3 [stat.ap] 8 Sep 2017
A self-excited threshold autoregressive state-space model for menstrual cycles: forecasting menstruation and identifying ovarian phases based on basal body temperature arxiv:177.6452v3 [stat.ap] 8 Sep
More informationAnalysing longitudinal data when the visit times are informative
Analysing longitudinal data when the visit times are informative Eleanor Pullenayegum, PhD Scientist, Hospital for Sick Children Associate Professor, University of Toronto eleanor.pullenayegum@sickkids.ca
More informationTime Series: Theory and Methods
Peter J. Brockwell Richard A. Davis Time Series: Theory and Methods Second Edition With 124 Illustrations Springer Contents Preface to the Second Edition Preface to the First Edition vn ix CHAPTER 1 Stationary
More informationEstimation of Optimally-Combined-Biomarker Accuracy in the Absence of a Gold-Standard Reference Test
Estimation of Optimally-Combined-Biomarker Accuracy in the Absence of a Gold-Standard Reference Test L. García Barrado 1 E. Coart 2 T. Burzykowski 1,2 1 Interuniversity Institute for Biostatistics and
More informationarxiv: v1 [stat.ap] 8 Jun 2016
The forecasting of menstruation based on a state-space modeling of basal body temperature time series Keiichi Fukaya, 1,4 Ai Kawamori, 1 Yutaka Osada, 2 Masumi Kitazawa, 3 and Makio Ishiguro 1 1 The Institute
More informationCovariance modelling for longitudinal randomised controlled trials
Covariance modelling for longitudinal randomised controlled trials G. MacKenzie 1,2 1 Centre of Biostatistics, University of Limerick, Ireland. www.staff.ul.ie/mackenzieg 2 CREST, ENSAI, Rennes, France.
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 informationCOLLABORATION OF STATISTICAL METHODS IN SELECTING THE CORRECT MULTIPLE LINEAR REGRESSIONS
American Journal of Biostatistics 4 (2): 29-33, 2014 ISSN: 1948-9889 2014 A.H. Al-Marshadi, This open access article is distributed under a Creative Commons Attribution (CC-BY) 3.0 license doi:10.3844/ajbssp.2014.29.33
More informationOn the Behavior of Marginal and Conditional Akaike Information Criteria in Linear Mixed Models
On the Behavior of Marginal and Conditional Akaike Information Criteria in Linear Mixed Models Thomas Kneib Department of Mathematics Carl von Ossietzky University Oldenburg Sonja Greven Department of
More informationA brief introduction to mixed models
A brief introduction to mixed models University of Gothenburg Gothenburg April 6, 2017 Outline An introduction to mixed models based on a few examples: Definition of standard mixed models. Parameter estimation.
More informationA Course in Applied Econometrics Lecture 14: Control Functions and Related Methods. Jeff Wooldridge IRP Lectures, UW Madison, August 2008
A Course in Applied Econometrics Lecture 14: Control Functions and Related Methods Jeff Wooldridge IRP Lectures, UW Madison, August 2008 1. Linear-in-Parameters Models: IV versus Control Functions 2. Correlated
More informationAnalysis of Longitudinal Data: Comparison between PROC GLM and PROC MIXED.
Analysis of Longitudinal Data: Comparison between PROC GLM and PROC MIXED. Maribeth Johnson, Medical College of Georgia, Augusta, GA ABSTRACT Longitudinal data refers to datasets with multiple measurements
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 informationSome general observations.
Modeling and analyzing data from computer experiments. Some general observations. 1. For simplicity, I assume that all factors (inputs) x1, x2,, xd are quantitative. 2. Because the code always produces
More informationSome methods for handling missing values in outcome variables. Roderick J. Little
Some methods for handling missing values in outcome variables Roderick J. Little Missing data principles Likelihood methods Outline ML, Bayes, Multiple Imputation (MI) Robust MAR methods Predictive mean
More informationAnalysis of Longitudinal Data: Comparison Between PROC GLM and PROC MIXED. Maribeth Johnson Medical College of Georgia Augusta, GA
Analysis of Longitudinal Data: Comparison Between PROC GLM and PROC MIXED Maribeth Johnson Medical College of Georgia Augusta, GA Overview Introduction to longitudinal data Describe the data for examples
More informationDirected acyclic graphs and the use of linear mixed models
Directed acyclic graphs and the use of linear mixed models Siem H. Heisterkamp 1,2 1 Groningen Bioinformatics Centre, University of Groningen 2 Biostatistics and Research Decision Sciences (BARDS), MSD,
More informationBayesian adaptive regression splines for hierarchical data
Bayesian adaptive regression splines for hierarchical data Jamie Lynn Bigelow and David B. Dunson Biostatistics Branch National Institute of Environmental Health Sciences Research Triangle Park, NC 27709
More informationBiostatistics 301A. Repeated measurement analysis (mixed models)
B a s i c S t a t i s t i c s F o r D o c t o r s Singapore Med J 2004 Vol 45(10) : 456 CME Article Biostatistics 301A. Repeated measurement analysis (mixed models) Y H Chan Faculty of Medicine National
More informationStatistics 262: Intermediate Biostatistics Model selection
Statistics 262: Intermediate Biostatistics Model selection Jonathan Taylor & Kristin Cobb Statistics 262: Intermediate Biostatistics p.1/?? Today s class Model selection. Strategies for model selection.
More informationStandard Errors & Confidence Intervals. N(0, I( β) 1 ), I( β) = [ 2 l(β, φ; y) β i β β= β j
Standard Errors & Confidence Intervals β β asy N(0, I( β) 1 ), where I( β) = [ 2 l(β, φ; y) ] β i β β= β j We can obtain asymptotic 100(1 α)% confidence intervals for β j using: β j ± Z 1 α/2 se( β j )
More informationOptimal Hypothesis Testing: From Semi to Fully Bayes Factors
Noname manuscript No. (will be inserted by the editor) Optimal Hypothesis Testing: From Semi to Fully Bayes Factors Optimal Hypothesis Testing Albert Vexler Chengqing Wu Kai Fun Yu Received: date / Accepted:
More informationTwo Applications of Nonparametric Regression in Survey Estimation
Two Applications of Nonparametric Regression in Survey Estimation 1/56 Jean Opsomer Iowa State University Joint work with Jay Breidt, Colorado State University Gerda Claeskens, Université Catholique de
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 informationChapter 3: Regression Methods for Trends
Chapter 3: Regression Methods for Trends Time series exhibiting trends over time have a mean function that is some simple function (not necessarily constant) of time. The example random walk graph from
More informationElements of Multivariate Time Series Analysis
Gregory C. Reinsel Elements of Multivariate Time Series Analysis Second Edition With 14 Figures Springer Contents Preface to the Second Edition Preface to the First Edition vii ix 1. Vector Time Series
More informationIntroduction to General Linear Models
Introduction to General Linear Models Moo K. Chung University of Wisconsin-Madison mkchung@wisc.edu September 27, 2014 In this chapter, we introduce general linear models (GLM) that have been widely used
More informationP -spline ANOVA-type interaction models for spatio-temporal smoothing
P -spline ANOVA-type interaction models for spatio-temporal smoothing Dae-Jin Lee 1 and María Durbán 1 1 Department of Statistics, Universidad Carlos III de Madrid, SPAIN. e-mail: dae-jin.lee@uc3m.es and
More informationCovariance Models (*) X i : (n i p) design matrix for fixed effects β : (p 1) regression coefficient for fixed effects
Covariance Models (*) Mixed Models Laird & Ware (1982) Y i = X i β + Z i b i + e i Y i : (n i 1) response vector X i : (n i p) design matrix for fixed effects β : (p 1) regression coefficient for fixed
More informationRANDOM and REPEATED statements - How to Use Them to Model the Covariance Structure in Proc Mixed. Charlie Liu, Dachuang Cao, Peiqi Chen, Tony Zagar
Paper S02-2007 RANDOM and REPEATED statements - How to Use Them to Model the Covariance Structure in Proc Mixed Charlie Liu, Dachuang Cao, Peiqi Chen, Tony Zagar Eli Lilly & Company, Indianapolis, IN ABSTRACT
More informationGOODNESS-OF-FIT TESTS FOR ARCHIMEDEAN COPULA MODELS
Statistica Sinica 20 (2010), 441-453 GOODNESS-OF-FIT TESTS FOR ARCHIMEDEAN COPULA MODELS Antai Wang Georgetown University Medical Center Abstract: In this paper, we propose two tests for parametric models
More informationNon-maximum likelihood estimation and statistical inference for linear and nonlinear mixed models
Optimum Design for Mixed Effects Non-Linear and generalized Linear Models Cambridge, August 9-12, 2011 Non-maximum likelihood estimation and statistical inference for linear and nonlinear mixed models
More informationMixed Effects Multivariate Adaptive Splines Model for the Analysis of Longitudinal and Growth Curve Data. Heping Zhang.
Mixed Effects Multivariate Adaptive Splines Model for the Analysis of Longitudinal and Growth Curve Data Heping Zhang Yale University NIA Workshop on Statistical Methods for Longitudinal Data on Aging
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 informationBayesian shrinkage approach in variable selection for mixed
Bayesian shrinkage approach in variable selection for mixed effects s GGI Statistics Conference, Florence, 2015 Bayesian Variable Selection June 22-26, 2015 Outline 1 Introduction 2 3 4 Outline Introduction
More informationLongitudinal Data Analysis. Michael L. Berbaum Institute for Health Research and Policy University of Illinois at Chicago
Longitudinal Data Analysis Michael L. Berbaum Institute for Health Research and Policy University of Illinois at Chicago Course description: Longitudinal analysis is the study of short series of observations
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 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 information