Mixture models for quantitative HIV RNA data

Abstract

Clinical investigators are increasing their use of quantitative determinations of HIV viral load in their study populations. The distributions of these measures may be highly skewed, left-censored, and with an extra spike below the detection limit of the assay. We recommend use of a mixture model in this situation, with two sets of explanatory covariates. We extend this model to incorporate multiple measures across time, and to employ shared parameters as a way of increasing model efficiency and parsimony. Data from a cohort of HIV-infected men are used to illustrate these features, and simulations are performed to assess the utility of shared parameters.

Get full access to this article

View all access options for this article.

References

Paxton WB , Coombs RW , McElrath MJ et al. Longitudinal analysis of quantitative virologic measures in human immunodeficiency virus-infected subjects with ≥ 400 CD4 lymphocytes: implications for applying measurements to individual patients. National Institute of Allergy and Infectious Diseases AIDS Vaccine Evaluation Group . Journal of Infectious Diseases 1997; 175: 247-254 .

Hughes JP. Mixed effects models with censored data with application to HIV RNA levels . Biometrics 1999; 55: 625-629 .

Ding AA , Wu H. A comparison study of models and fitting procedures for biphasic viral dynamics in HIV-1 infected patients treated with antiviral therapies . Biometrics 2000; 56: 293-300 .

Wu H , Ding AA , De Gruttola, V. Estimation of HIV dynamic parameters . Statistics in Medicine 1998; 17: 2463-2485 .

Farewell VT. The use of mixture models for the analysis of survival data with long-term survivors . Biometrics 1982; 38: 1041-1046 .

Moulton LH , Halsey NA. A mixture model with detection limits for regression analyses of antibody response to vaccine . Biometrics 1995; 51: 1570-1578 .

Moulton LH , Halsey NA. A mixed gamma model for regression analyses of quantitative assay data . Vaccine 1996; 14: 1154-1158 .

Hall DB. Zero-inflated Poisson and binomial regression with random effects: a case study . Biometrics 2000; 56: 1030-1039 .

Zhou X-H , Tu W. Comparison of several independent population means when their samples contain log-normal and possibly zero observations . Biometrics 1999; 55: 645-651 .

10.

Prentice RL. A log gamma model and its maximum likelihood estimation . Biometrika 1974; 61: 539-544 .

11.

Lachenbruch PA. Comparisons of two-part models with competitors . Abstract, Joint Statistical Meetings, Indianapolis, 2000.

12.

Wu MC , Carroll RJ. Estimation and comparison of changes in the presence of informative right censoring by modeling the censoring process . Biometrics 1988; 44: 175-188 .

13.

Ten Have T , Kunselman AR , Pulkstenis EP , Landis JR. Mixed effects logistic regression models for longitudinal binary response data with informative drop-out . Biometrics 1998; 54: 367-383 .

14.

Albert PS , Follmann DA , Barnhart HX. A generalized estimating equation approach for modeling random length binary vector data . Biometrics 1997; 53: 1116-1124 .

15.

Rotnitzky A , Jewell NP. Hypothesis testing of regression parameters in semiparametric generalized linear models for cluster correlated data . Biometrika 1990; 77: 485-497 .

16.

Huber PJ. The behavior of maximum likelihood estimates under non-standard conditions . In: Proceedings of the Fifth Berkeley Symposium on Mathematical Statistics and Probability. Berkeley: University of California Press, 1967.

17.

Royall RM. Model robust confidence intervals using maximum likelihood estimators . International Statistical Review 1986; 54: 221-226 .

18.

Aptech Systems, Inc. The GAUSS System; Applications: Maximum Likelihood. (Version 4.0). Maple Valley, WA: Aptech Systems, Inc ., 2000.

19.

Barroso PF , Schechter M , Gupta P et al. Effect of antiretroviral therapy on HIV shedding in semen . Annals of Internal Medicine 2000; 133: 280-284 .

20.

Berk K , Lachenbruch PA. Repeated measures with zeroes . Abstract, Joint Statistical Meetings, Indianapolis, 2000.