A joint modeling and estimation method for multivariate longitudinal data with mixed types of responses to analyze physical activity data generated by accelerometers

A mixed effect model is proposed to jointly analyze multivariate longitudinal data with continuous, proportion, count, and binary responses. The association of the variables is modeled through the correlation of random effects. We use a quasi-likelihood type approximation for nonlinear variables and transform the proposed model into a multivariate linear mixed model framework for estimation and inference. Via an extension to the EM approach, an efficient algorithm is developed to fit the model. The method is applied to physical activity data, which uses a wearable accelerometer device to measure daily movement and energy expenditure information. Our approach is also evaluated by a simulation study.

[1]  D. Rubin,et al.  The ECME algorithm: A simple extension of EM and ECM with faster monotone convergence , 1994 .

[2]  R. Wolfinger,et al.  Generalized linear mixed models a pseudo-likelihood approach , 1993 .

[3]  G. Molenberghs,et al.  Models for Discrete Longitudinal Data , 2005 .

[4]  Henrik Holmberg,et al.  Generalized linear models with clustered data: Fixed and random effects models , 2011, Comput. Stat. Data Anal..

[5]  Harvey Goldstein,et al.  Improved Approximations for Multilevel Models with Binary Responses , 1996 .

[6]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[7]  Runze Li,et al.  Statistical models for longitudinal zero‐inflated count data with applications to the substance abuse field , 2012, Statistics in medicine.

[8]  D. Bates,et al.  Fitting Linear Mixed-Effects Models Using lme4 , 2014, 1406.5823.

[9]  S. Ferrari,et al.  Beta Regression for Modelling Rates and Proportions , 2004 .

[10]  N. Breslow,et al.  Approximate inference in generalized linear mixed models , 1993 .

[11]  Geert Verbeke,et al.  Pairwise Fitting of Mixed Models for the Joint Modeling of Multivariate Longitudinal Profiles , 2006, Biometrics.

[12]  D. Bates,et al.  Newton-Raphson and EM Algorithms for Linear Mixed-Effects Models for Repeated-Measures Data , 1988 .

[13]  Geert Verbeke,et al.  High dimensional multivariate mixed models for binary questionnaire data , 2006 .

[14]  D. Bates,et al.  Nonlinear mixed effects models for repeated measures data. , 1990, Biometrics.

[15]  John Staudenmayer,et al.  Changes in sedentary time and physical activity in response to an exercise training and/or lifestyle intervention. , 2014, Journal of physical activity & health.

[16]  N. Goldman,et al.  Improved estimation procedures for multilevel models with binary response: a case‐study , 2001 .

[17]  A. Agresti,et al.  A Correlated Probit Model for Joint Modeling of Clustered Binary and Continuous Responses , 2001 .

[18]  Geert Molenberghs,et al.  The analysis of multivariate longitudinal data: A review , 2014, Statistical methods in medical research.

[19]  Raymond J Carroll,et al.  Hierarchical functional data with mixed continuous and binary measurements , 2014, Biometrics.

[20]  R. Carroll,et al.  An Evaluation of Accelerometer-derived Metrics to Assess Daily Behavioral Patterns , 2017, Medicine and science in sports and exercise.

[21]  Dibyen Majumdar,et al.  Conditional Second-Order Generalized Estimating Equations for Generalized Linear and Nonlinear Mixed-Effects Models , 2002 .