Shared random effects analysis of multi‐state Markov models: application to a longitudinal study of transitions to dementia

Multi-state models are appealing tools for analysing data about the progression of a disease over time. In this paper, we consider a multi-state Markov chain with two competing absorbing states: dementia and death and three transient non-demented states: cognitively normal, amnestic mild cognitive impairment (amnestic MCI), and non-amnestic mild cognitive impairment (non-amnestic MCI). The likelihood function for the data is derived and estimates for the effects of the covariates on transitions are determined when the process can be viewed as a polytomous logistic regression model with shared random effects. The presence of a shared random effect not only complicates the formulation of the likelihood but also its evaluation and maximization. Three approaches for maximizing the likelihood are compared using a simulation study; the first method is based on the Gauss-quadrature technique, the second method is based on importance sampling ideas, and the third method is based on an expansion by Taylor series. The best approach is illustrated using a longitudinal study on a cohort of cognitively normal subjects, followed annually for conversion to mild cognitive impairment (MCI) and/or dementia, conducted at the Sanders Brown Center on Aging at the University of Kentucky.

[1]  Marco Alfò,et al.  Regression models for binary longitudinal responses , 1998, Stat. Comput..

[2]  B. Reboussin,et al.  Mixed Effects Logistic Regression Models for Longitudinal Ordinal Functional Response Data with Multiple‐Cause Drop‐Out from the Longitudinal Study of Aging , 2000, Biometrics.

[3]  Sujuan Gao,et al.  A shared random effect parameter approach for longitudinal dementia data with non‐ignorable missing data , 2004, Statistics in medicine.

[4]  Milton Abramowitz,et al.  Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables , 1964 .

[5]  L. V. Rubinstein,et al.  Markov models for covariate dependence of binary sequences. , 1985, Biometrics.

[6]  Adele A. Barrett,et al.  Elements of Applied Stochastic Processes (2nd Edition) , 1984 .

[7]  S. Raudenbush,et al.  Maximum Likelihood for Generalized Linear Models with Nested Random Effects via High-Order, Multivariate Laplace Approximation , 2000 .

[8]  David R. Cox,et al.  Nonlinear component of variance models , 1992 .

[9]  Jong-Ling Fuh,et al.  Measuring Alzheimer's disease progression with transition probabilities in the Taiwanese population , 2004, International journal of geriatric psychiatry.

[10]  S. DeKosky,et al.  Mild cognitive impairment, amnestic type , 2004, Neurology.

[11]  J. Ashford,et al.  “Preclinical” AD revisited , 2000, Neurology.

[12]  P. Diggle,et al.  Analysis of Longitudinal Data , 2003 .

[13]  Sophia Rabe-Hesketh,et al.  Multilevel logistic regression for polytomous data and rankings , 2003 .

[14]  F B Hu,et al.  Random-effects regression analysis of correlated grouped-time survival data , 2000, Statistical methods in medical research.

[15]  M. Tanner.  Tools for statistical inference: methods for the exploration of posterior distributions and likeliho , 1994 .

[16]  R. Petersen.  Mild cognitive impairment as a diagnostic entity , 2004, Journal of internal medicine.

[17]  C. McCulloch.  Maximum Likelihood Algorithms for Generalized Linear Mixed Models , 1997 .

[18]  F. Schmitt,et al.  Models of progression in AD , 2001, Neurology.

[19]  L. Beckett,et al.  Analysis of change in self-reported physical function among older persons in four population studies. , 1996, American journal of epidemiology.

[20]  F. Schmitt,et al.  Brain donation in normal aging: procedures, motivations, and donor characteristics from the Biologically Resilient Adults in Neurological Studies (BRAiNS) Project. , 2001, The Gerontologist.