Joint modeling of HIV data in multicenter observational studies: A comparison among different approaches

Abstract

Disease process over time results from the combination of event history information and longitudinal process. Commonly, separate analyses of longitudinal and survival outcomes are performed. However, discharging the dependence between these components may cause misleading results. Separate analyses are difficult to interpret whenever one deals with observational retrospective multicenter cohort studies where the biomarkers are poorly monitored over time, while the survival component may be affected by several sources of bias, such as multiple endpoints, multiple time-scales, and informative censoring. We discuss how joint modeling of longitudinal and survival data represents an effective strategy to incorporate all information simultaneously and to provide valid and efficient inferences, thus allowing to produce a better insight into the biological mechanisms underlying the phenomenon under study. Accounting for the whole dynamics of the disease process is crucial in retrospective longitudinal studies. In this work, we present different approaches for modeling longitudinal and time-to-event data, retrieved from 648 HIV-infected patients enrolled in the Italian cohort of the CASCADE (Concerted Action on SeroConversion to AIDS and Death in Europe) study, one of the largest AIDS collaborative cohort studies. In particular, we evaluate CD4 lymphocyte evolution over time (from the date of seroconversion) and overall survival, CD4 being one of the most important immunologic biomarker for HIV progression. Besides a standard separate modeling approach, we consider two alternative joint models: the traditional joint model and the joint latent class mixed model. Advantages and disadvantages of the different approaches are discussed. To compare the performances of these models, cross-validation procedures are also performed.

Keywords

CASCADE study Cox regression cross-validation joint modeling latent class mixed model linear mixed-effect model

Get full access to this article

View all access options for this article.

References

Jewell

Kalbfleisch

. Marker processes in survival analysis. Lifetime Data Anal 1996; 2: 15–29.

Henderson

Diggle

Dobson

. Joint modeling of longitudinal measurements and event time data. Biostatistics 2000; 4: 465–480.

Hogan

Daniels

. A hierarchical modelling approach to analysing longitudinal data with drop-out and non-compliance, with application to an equivalence trial in paediatric acquired immune defficiency syndrome. Appl Stat 2002; 51: 1–21.

Dupuy

Mesbah

. Joint modeling of event time and nonignorable missing longitudinal data. Lifetime Data Anal 2002; 8: 99–115.

Verzilli

Carpenter

. A Monte-Carlo EM algorithm for random-coefficient-based dropout models. J Appl Stat 2002; 29: 1011–1021.

Dobson

Henderson

. Diagnostics for joint longitudinal and dropout time modelling. Biometrics 2003; 59: 741–751.

Wood

White

Hillsdon

et al.

Comparison of imputation and modelling methods in the analysis of a physical activity trial with missing outcomes. Int J Epidemiol 2005; 34: 89–99.

Roy

Daniels

. A general class of pattern mixture models for nonignorable dropout with many possible dropout times. Biometrics 2008; 64: 538–545.

Wulfsohn

Tsiatis

. A joint model for survival and longitudinal data measured with error. Biometrics 1997; 53: 330–339.

10.

Zeger

. Joint analysis of longitudinal data comprising repeated measures and times to events. Appl Stat 2001; 50: 375–387.

11.

Wang

Taylor

JMG

. Joint modeling longitudinal and event time data with application to acquired immunodefi-ciency syndrome. J Am Stat Assoc 2001; 96: 895–905.

12.

Song

Davidian

Tsiatis

. A semiparametric likelihood approach to joint modeling of longitudinal and time-to-event data. Biometrics 2002; 58: 742–753.

13.

Brown

Ibrahim

. A bayesian semiparametric joint hierarchical model for longitudinal and survival data. Biometrics 2003; 59: 221–228.

14.

Tseng

Hsieh

Wang

. Joint modelling of accelerated failure time and longitudinal data. Biometrika 2005; 92: 587–603.

15.

Lin

Taylor

JMG

. Semiparametric modeling of longitudinal measurements and time-to-event data–a two-stage regression calibration approach. Biometrics 2008; 64: 1238–1246.

16.

Zeng

Cai

. Simultaneous modelling of survival and longitudinal data with an application to repeated quality of life measures. Lifetime Data Anal 2005; 11: 151–174.

17.

Elashoff

. An approach to joint analysis of longitudinal measurements and competing risks failure time data. Stat Med 2007; 26: 2813–2835.

18.

Elashoff

. A joint model for longitudinal measurements and survival data in the presence of multiple failure types. Biometrics 2008; 64: 762–771.

19.

Huang

Elashoff

. A joint model of longitudinal and competing risks survival data with heterogeneous random effects and outlying longitudinal measurements. Stat Interface 2010; 3: 185–195.

20.

Liu

et al.

Analysis of longitudinal and survival data: joint modeling, inference methods, and issues. J Prob Stat 2012; 2012: 17 pp–17 pp. DOI: 10.1155/2012/640153.

21.

Lin

Turnbull

McCulloch

et al.

Latent class models for joint analysis of longitudinal biomarker and event process data. J Am Stat Assoc 2002; 97: 53–65.

22.

Proust-Lima

Taylor

JMG

. Development and validation of a dynamic prognostic tool for prostate cancer recurrence using repeated measures of post-treatment PSA: a joint modelling approach. Biostatistics 2009; 10: 535–549.

23.

Proust-Lima

Joly

Jacqmin-Gadda

. Joint modelling of multivariate longitudinal outcomes and a time-to-event: a nonlinear latent class approach. Comput Stat Data Anal 2009; 53: 1142–1154.

24.

The CASCADE Collaboration. Survival after introduction of HAART in people with known duration of HIV-1 infection. Lancet 2000; 355: 1158–1159.

25.

Zugna

Geskus

De Stavola

et al.

Time to virological failure, treatment change and interruption for individuals treated within 12 months of HIV seroconversion and in chronic infection. Antivir Ther 2012; 17: 1039–1048.

26.

Thiébaut

Jacqmin-Gadda

Babiker

et al.

Joint modelling of bivariate longitudinal data with informative dropout and left-censoring, with application to the evolution of CD4+ cell count and HIV RNA viral load in response to treatment of HIV infection. Stat Med 2005; 24: 65–82.

27.

Laird

Ware

. Random effects models for longitudinal data. Biometrics 1982; 38: 963–974.

28.

Pinheiro

Bates

. Mixed-effects models in S and S-PLUS, New York: Springer, 2000.

29.

Mellors

Muñoz

Giorgi

et al.

Plasma viral load and CD4+ lymphocytes as prognostic markers of HIV-1 infection. Ann Intern Med 1997; 126: 946–954.

30.

Lichtenstein

Armon

Buchacz

et al.

Initiation of antiretroviral therapy at CD4 cell counts ≥350 cells/mm³ does not increase incidence or risk of peripheral neuropathy, anemia, or renal insufficiency. J Acquir Immune Defic Syndr 2008; 47: 27–35.

31.

Balter

. How does HIV overcome the body's T-cell bodyguards? Science 1997; 278: 1399–1400.

32.

McKinnon

Kimani

Wachihi

et al.

Effect of baseline HIV disease parameters on CD4+ T cell recovery after antiretroviral therapy initiation in Kenyan women. PLoS One 2010; 5: e11434–e11434.

33.

Rizopoulos

. JM: an R package for the joint modelling of longitudinal and time-to-event data. J Stat Softw 2010; 35: 1–33.

34.

Zhang

Lin

et al.

Interleukin-17-producing CD4(+) T cells increase with severity of liver damage in patients with chronic hepatitis B. Hepatology 2010; 51: 81–91.

35.

Sajadi

Pulijala

Redfield

et al.

Chronic immune activation and decreased CD4 cell counts associated with hepatitis C infection in HIV-1 natural viral suppressors. AIDS 2012; 26: 1879–1884.

36.

Branch

Van Natta

Vachon

et al.

Mortality in HCV-infected patient with a diagnosis of AIDS in the era of combination anti-retroviral therapy. Clin Infect Dis 2012; 55: 137–144.

37.

Chun

Roediger

Hullsiek

et al.

Hepatitis B virus coinfection negatively impacts on HIV outcomes in seroconverters. J Infect Dis 2012; 205: 185–193.

38.

Peters

Marston

. Preventing deaths in person with HIV/hepatitis B virus coinfection: a call to accelerate prevention and treatment efforts. J Infect Dis 2012; 205: 166–168.

39.

Proust-Lima

Sene

Taylor

JMG

et al.

Joint latent class models for longitudinal and time-to-event data: a review. Stat Methods Med Res 2012; 0: 1–17.

40.

Jacqmin-Gadda

Proust-Lima

Taylor

JMG

et al.

Score test for conditional independence between longitudinal outcome and time to event given the classes in the joint latent class model. Biometrics 2010; 66: 11–19.

41.

Graf

Schmoor

Sauerbrei

et al.

Assessment and comparison of prognostic classification schemes for survival data. Stat Med 1999; 18: 2529–2545.