Abstract

Longitudinal HIV-1 RNA viral load measures are often subject to censoring due to upper and lower detection limits depending on the quantification assays. A complication arises when these continuous measures present a heavy-tailed behavior because inference can be seriously affected by the misspecification of their parametric distribution. For such data structures, we propose a robust nonlinear censored regression model based on the scale mixtures of normal distributions. By taking into account the autocorrelation existing among irregularly observed measures, a damped exponential correlation structure is considered. A stochastic approximation of the EM algorithm is developed to obtain the maximum likelihood estimates of the model parameters. The main advantage of this new procedure os to allow estimating the parameters of interest and evaluating the log-likelihood function easily and quickly. Furthermore, the standard errors of the fixed effects and predictions of unobservable values of the response can be obtained as a byproduct. The practical utility of the proposed method is exemplified using both simulated and real data.