title: Fast Implementation for Normal Mixed Effects Models With Censored Response.
abstract: We propose an EM algorithm for computing the maximum likelihood and restricted maximum likelihood for linear and nonlinear mixed effects models with censored response. In contrast with previous developments, this algorithm uses closed-form expressions at the E-step, as opposed to Monte Carlo simulation. These expressions rely on formulas for the mean and variance of a truncated multinormal distribution, and can be computed using available software. This leads to an improvement in the speed of computation of up to an order of magnitude. A wide class of mixed effects models is considered, including the Laird-Ware model, and extensions to different structures for the variance components, heteroscedastic and autocorrelated errors, and multilevel models. We apply the methodology to two case studies from our own biostatistical practice, involving the analysis of longitudinal HIV viral load in two recent AIDS studies. The proposed algorithm is implemented in the R package lmec. An appendix which includes further mathematical details, the R code, and datasets for examples and simulations are available as the online supplements.
authors.affiliation: Department of Family and Preventive Medicine, UC San Diego School of Medicine, La Jolla, CA 92093-0717.
authors.forename: Florin
authors.initials: F
authors.lastname: Vaida
authors.affiliation: Department of Family and Preventive Medicine, UC San Diego School of Medicine, La Jolla, CA 92093-0717.
authors.forename: Lin
authors.initials: L
authors.lastname: Liu
keywords: Detection limit; EM algorithm; HIV viral load; Maximum likelihood; Truncated multinormal distribution
references.citation: Neurology. 2009 Aug 4;73(5):342-8
references.pmid: 19474412
references.citation: Biometrics. 1999 Jun;55(2):625-9
references.pmid: 11318225
references.citation: Comput Stat Data Anal. 2007 Aug 15;51(12):5718-5730
references.pmid: 19578533
references.citation: Biometrics. 1990 Sep;46(3):673-87
references.pmid: 2242409
references.citation: Pediatrics. 2008 Mar;121(3):e513-21
references.pmid: 18310171