IntroductionCombination antiretroviral therapy has become the standard care of human immunodeficiency virus (HIV)-infected patients and has further led to a dramatically decreased progression probability to acquired immune deficiency syndrome (AIDS) for patients under such a therapy. However, responses of the patients to this therapy have recorded heterogeneous complexity and high dynamism. In this paper, we simultaneously model long-term viral suppression, viral rebound, and state-specific duration of HIV-infected patients.MethodsFull-parametric and semi-parametric Markov multistate models were applied to assess the effects of covariates namely TB co-infection, educational status, marital status, age, quality of life (QoL) scores, white and red blood cell parameters, and liver enzyme abnormality on long-term viral suppression, viral rebound and state-specific duration for HIV-infected individuals before and after treatment. Furthermore, two models, one including and another excluding the effect of the frailty, were presented and compared in this study.ResultsResults from the diagnostic plots, Akaike information criterion (AIC) and likelihood ratio test showed that the Weibull multistate frailty model fitted significantly better than the exponential and semi-parametric multistate models. Viral rebound was found to be significantly associated with many sex partners, higher eosinophils count, younger age, lower educational level, higher monocyte counts, having abnormal neutrophils count, and higher liver enzyme abnormality. Furthermore, viral suppression was also found to be significantly associated with higher QoL scores, and having a stable sex partner. The analysis result also showed that patients with a stable sex partner, higher educational levels, higher QoL scores, lower eosinophils count, lower monocyte counts, and higher RBC indices were more likely to spend more time in undetectable viral load state.ConclusionsTo achieve and maintain the UNAIDS 90% suppression targets, additional interventions are required to optimize antiretroviral therapy outcomes, specifically targeting those with poor clinical characteristics, lower education, younger age, and those with many sex partners. From a methodological perspective, the parametric multistate approach with frailty is a flexible approach for modeling time-varying variables, allowing for dealing with heterogeneity between the sequence of transitions, as well as allowing for a reasonable degree of flexibility with a few additional parameters, which then aids in gaining a better insight into how factors change over time.