Joint analysis of repeatedly observed continuous and ordinal measures of disease severity

In biomedical studies often multiple measures of disease severity are recorded over time. Although correlated, such measures are frequently analysed separately of one another. Joint analysis of the outcomes variables has several potential advantages over separate analyses. However, models for response variables of different types (discrete and continuous) are challenging to define and to fit. Herein we propose correlated probit models for joint analysis of repeated measurements on ordinal and continuous variables measuring the same underlying disease severity over time. We demonstrate how to rewrite the models so that maximum-likelihood estimation and inference can be performed with standard software. Simulation studies are performed to assess efficiency gains in fitting the responses together rather than separately and to guide response variable selection for future studies. Data from a depression clinical trial are used for illustration.

[1]  Xihong Lin,et al.  Latent Variable Models for Longitudinal Data with Multiple Continuous Outcomes , 2000, Biometrics.

[2]  D. Charney,et al.  Addition of the α2-Antagonist Yohimbine to Fluoxetine: Effects on Rate of Antidepressant Response , 2004, Neuropsychopharmacology.

[3]  Y. Matsuyama,et al.  Mixed models for bivariate response repeated measures data using Gibbs sampling. , 1997, Statistics in medicine.

[4]  David B Dunson,et al.  A Bayesian Approach for Joint Modeling of Cluster Size and Subunit‐Specific Outcomes , 2003, Biometrics.

[5]  Alan Agresti.  A Model for Repeated Measurements of a Multivariate Binary Response , 1997 .

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

[7]  Nan M. Laird,et al.  Regression Models for a Bivariate Discrete and Continuous Outcome with Clustering , 1995 .

[8]  J. Markowitz,et al.  Self-Reported Depressive Symptom Measures: Sensitivity to Detecting Change in a Randomized, Controlled Trial of Chronically Depressed, Nonpsychotic Outpatients , 2005, Neuropsychopharmacology.

[9]  Jason Roy,et al.  Analysis of Multivariate Longitudinal Outcomes With Nonignorable Dropouts and Missing Covariates , 2002 .

[10]  R. Gueorguieva,et al.  A multivariate generalized linear mixed model for joint modelling of clustered outcomes in the exponential family , 2001 .

[11]  Qing Liu,et al.  A note on Gauss—Hermite quadrature , 1994 .

[12]  Geert Molenberghs,et al.  Two latent variable risk assessment approaches for mixed continuous and discrete outcomes from developmental toxicity data , 2001 .

[13]  P J Catalano,et al.  Bivariate modelling of clustered continuous and ordered categorical outcomes. , 1997, Statistics in medicine.

[14]  J. Booth,et al.  Standard Errors of Prediction in Generalized Linear Mixed Models , 1998 .

[15]  G Molenberghs,et al.  Multivariate probit analysis: a neglected procedure in medical statistics. , 1991, Statistics in medicine.

[16]  A. Agresti,et al.  Measures of relative model fit , 2002 .

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

[18]  N M Laird,et al.  Regression models for mixed discrete and continuous responses with potentially missing values. , 1997, Biometrics.

[19]  M M Regan,et al.  Likelihood Models for Clustered Binary and Continuous Out comes: Application to Developmental Toxicology , 1999, Biometrics.

[20]  Louise Ryan,et al.  Bivariate Latent Variable Models for Clustered Discrete and Continuous Outcomes , 1992 .

[21]  D. Dunson,et al.  Bayesian latent variable models for clustered mixed outcomes , 2000 .

[22]  D. Heitjan,et al.  Modelling repeated-series longitudinal data. , 1997, Statistics in medicine.

[23]  L. Ryan,et al.  Latent Variable Models for Mixed Discrete and Continuous Outcomes , 1997 .

[24]  L Ryan,et al.  Multivariate linear mixed models for multiple outcomes. , 1999, Statistics in medicine.

[25]  Scott L. Zeger,et al.  A Class of Logistic Regression Models for Multivariate Binary Time Series , 1989 .

[26]  Gregory C. Reinsel,et al.  Estimation and Prediction in a Multivariate Random Effects Generalized Linear Model , 1984 .

[27]  P. Catalano,et al.  Regression Models and Risk Estimation for Mixed Discrete and Continuous Outcomes in Developmental Toxicology , 2000, Risk analysis : an official publication of the Society for Risk Analysis.

[28]  L Ryan,et al.  A Scaled Linear Mixed Model for Multiple Outcomes , 2000, Biometrics.

[29]  David A. Schoenfeld,et al.  A Random-Effects Model for Multiple Characteristics with Possibly Missing Data , 1997 .

[30]  J Rochon,et al.  Analyzing bivariate repeated measures for discrete and continuous outcome variables. , 1996, Biometrics.

[31]  Carla Rampichini,et al.  Alternative Specifications of Multivariate Multilevel Probit Ordinal Response Models , 2003 .

[32]  D. Stram,et al.  Variance components testing in the longitudinal mixed effects model. , 1994, Biometrics.

[33]  Geert Molenberghs,et al.  The Use of Score Tests for Inference on Variance Components , 2003, Biometrics.

[34]  B A Rosner,et al.  Analysis of longitudinally observed irregularly timed multivariate outcomes: regression with focus on cross-component correlation. , 2001, Statistics in medicine.

[35]  S. Rabe-Hesketh,et al.  Reliable Estimation of Generalized Linear Mixed Models using Adaptive Quadrature , 2002 .