Premium
Modelling the random effects covariance matrix in longitudinal data
Author(s) -
Daniels Michael J.,
Zhao Yan D.
Publication year - 2003
Publication title -
statistics in medicine
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.996
H-Index - 183
eISSN - 1097-0258
pISSN - 0277-6715
DOI - 10.1002/sim.1470
Subject(s) - cholesky decomposition , estimation of covariance matrices , covariance matrix , mathematics , covariance , gibbs sampling , random effects model , estimator , statistics , matrix (chemical analysis) , bayesian probability , eigenvalues and eigenvectors , medicine , physics , meta analysis , materials science , quantum mechanics , composite material
A common class of models for longitudinal data are random effects (mixed) models. In these models, the random effects covariance matrix is typically assumed constant across subject. However, in many situations this matrix may differ by measured covariates. In this paper, we propose an approach to model the random effects covariance matrix by using a special Cholesky decomposition of the matrix. In particular, we will allow the parameters that result from this decomposition to depend on subject‐specific covariates and also explore ways to parsimoniously model these parameters. An advantage of this parameterization is that there is no concern about the positive definiteness of the resulting estimator of the covariance matrix. In addition, the parameters resulting from this decomposition have a sensible interpretation. We propose fully Bayesian modelling for which a simple Gibbs sampler can be implemented to sample from the posterior distribution of the parameters. We illustrate these models on data from depression studies and examine the impact of heterogeneity in the covariance matrix on estimation of both fixed and random effects. Copyright © 2003 John Wiley & Sons, Ltd.