A mechanistic nonlinear model for censored and mismeasured covariates in longitudinal models, with application in AIDS studies.

When modeling longitudinal data, the true values of time-varying covariates may be unknown because of detection-limit censoring or measurement error. A common approach in the literature is to empirically model the covariate process based on observed data and then predict the censored values or mismeasured values based on this empirical model. Such an empirical model can be misleading, especially for censored values since the (unobserved) censored values may behave very differently than observed values due to the underlying data-generation mechanisms or disease status. In this paper, we propose a mechanistic nonlinear covariate model based on the underlying data-generation mechanisms to address censored values and mismeasured values. Such a mechanistic model is based on solid scientific or biological arguments, so the predicted censored or mismeasured values are more reasonable. We use a Monte Carlo EM algorithm for likelihood inference and apply the methods to an AIDS dataset, where viral load is censored by a lower detection limit. Simulation results confirm that the proposed models and methods offer substantial advantages over existing empirical covariate models for censored and mismeasured covariates.

[1]  Hulin Wu,et al.  Comparison of Two Indinavir/Ritonavir Regimens in the Treatment of HIV-Infected Individuals , 2004, Journal of acquired immune deficiency syndromes.

[2]  R. Grant,et al.  Sustained CD4+ T cell response after virologic failure of protease inhibitor-based regimens in patients with human immunodeficiency virus infection. , 2000, The Journal of infectious diseases.

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

[4]  M. Lederman,et al.  Immunologic responses associated with 12 weeks of combination antiretroviral therapy consisting of zidovudine, lamivudine, and ritonavir: results of AIDS Clinical Trials Group Protocol 315. , 1998, The Journal of infectious diseases.

[5]  Harry Joe,et al.  Accuracy of Laplace approximation for discrete response mixed models , 2008, Comput. Stat. Data Anal..

[6]  Lin Liu,et al.  Fast Implementation for Normal Mixed Effects Models With Censored Response , 2009, Journal of computational and graphical statistics : a joint publication of American Statistical Association, Institute of Mathematical Statistics, Interface Foundation of North America.

[7]  J. Hughes,et al.  Mixed Effects Models with Censored Data with Application to HIV RNA Levels , 1999, Biometrics.

[8]  J. Booth,et al.  Maximizing generalized linear mixed model likelihoods with an automated Monte Carlo EM algorithm , 1999 .

[9]  Joseph G. Ibrahim,et al.  Missing covariates in generalized linear models when the missing data mechanism is non‐ignorable , 1999 .

[10]  V De Gruttola,et al.  Estimation of HIV dynamic parameters. , 1998, Statistics in medicine.

[11]  Russell D. Wolfinger,et al.  Laplace's approximation for nonlinear mixed models. , 1993 .

[12]  Lang Wu,et al.  A Joint Model for Nonlinear Mixed-Effects Models With Censoring and Covariates Measured With Error, With Application to AIDS Studies , 2002 .

[13]  L. Palmisano,et al.  A brief history of antiretroviral therapy of HIV infection: success and challenges. , 2011, Annali dell'Istituto superiore di sanita.

[14]  H Wu,et al.  Population HIV‐1 Dynamics In Vivo: Applicable Models and Inferential Tools for Virological Data from AIDS Clinical Trials , 1999, Biometrics.

[15]  G. C. Wei,et al.  A Monte Carlo Implementation of the EM Algorithm and the Poor Man's Data Augmentation Algorithms , 1990 .

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

[17]  Paul W. Bernhardt,et al.  Flexible modeling of survival data with covariates subject to detection limits via multiple imputation , 2014, Comput. Stat. Data Anal..