 Research
 Open access
 Published:
The estimation of long and short term survival time and associated factors of HIV patients using mixture cure rate models
BMC Medical Research Methodology volumeÂ 23, ArticleÂ number:Â 123 (2023)
Abstract
Background
HIV is one of the deadliest epidemics and one of the most critical global public health issues. Some are susceptible to die among people living with HIV and some survive longer. The aim of the present study is to use mixture cure models to estimate factors affecting short and longterm survival of HIV patients.
Methods
The total sample size was 2170 HIVinfected people referred to the disease counseling centers in Kermanshah Province, in the west of Iran, from 1998 to 2019. A Semiparametric PH mixture cure model and a mixture cure frailty model were fitted to the data. Also, a comparison between these two models was performed.
Results
Based on the results of the mixture cure frailty model, antiretroviral therapy, tuberculosis infection, history of imprisonment, and mode of HIV transmission influenced shortterm survival time (pvalueâ€‰<â€‰0.05). On the other hand, prison history, antiretroviral therapy, mode of HIV transmission, age, marital status, gender, and education were significantly associated with longterm survival (pvalueâ€‰<â€‰0.05). The concordance criteria (Kindex) value for the mixture cure frailty model was 0.65 whereas for the semiparametric PH mixture cure model was 0.62.
Conclusion
This study showed that the frailty mixture cure models is more suitable in the situation where the studied population consisted of two groups, susceptible and nonsusceptible to the event of death. The people with a prison history, who received ART treatment, and contracted HIV through injection drug users survive longer. Health professionals should pay more attention to these findings in HIV prevention and treatment.
Introduction
Human immunodeficiency virus (HIV) is still one of the most critical public health issues globally, as almost 33 million people have died from the virus so far. By taking HIV medicine, antiretroviral therapy (ART), infected people can live a long and healthy life. Nevertheless, without treatment a person with HIV is more likely to develop a severe condition called acquired immunodeficiency syndrome (AIDS). If left untreated, the life expectancy at the final stage of AIDS is about three years [1].
At the end of 2019, about 38 million people were living with HIV, of whom 81% knew their status. Due to the gap in HIV services, 690,000 people died at the end of 2019 due to HIVrelated causes, and 1.7 million people have recently been infected with the virus [1,2,3]. Iran has been one of the active countries in fighting against HIV/AIDS in the Middle East for decades [4]. However, despite the efforts made, there are 59,000 (33,000â€“130,000) people living with HIV in this country, and 4,100 (1200â€“12,000) people have recently been infected with the virus, of which about 2,500 (1200â€“5600) have died of AIDS [5].
In the contemporary world, extensive research is being done to increase the lifetime of patients with a variety of diseases, including deadly diseases such as cancer and HIV / AIDS. As a result, many patients are expected to be cured permanently. They are longterm survivors and do not experience socalled death. These people are called cured patients [6]. In such a situation, where the followup period is long and a proportion of patients do not experience the event of interest (nonsusceptible), the rate of censorship increases, and as a result, the estimated survival rate of conventional methods such as the Cox model is overestimated because Cox model regards all the patients as susceptible [7, 8].
When a nonnegligible proportion of patients long survive (nonsusceptible), mixture cure models describe to survival process more appropriately [8]. For several reasons, mixture cure models can be an excellent alternative to the conventional Cox proportional hazard models for highcensorship data. First, proportional hazards may not be assumed when the survival curves in their sequence are flat (horizontal line). Second, survival curves with a long horizontal line in the sequence may indicate heterogeneity in patientsâ€™ populations, which can be helpful for explicitly describing mixture cure models [8]. Mixture cure models make it possible to estimate the proportion of longterm survivors and investigate the effect of covariates on short and longterm survival time [9]. In several studies, the cured models have been theoretically developed and in the proposed models, HIV data have been used [10,11,12,13]. Varshney et al.â€™s study estimated the proportion of cured people living with HIV/AIDS treated with ART using Bayesian methods [8].
Price and Manatunga have pointed out that in medical and epidemiological studies, there may be some unobserved covariates between individuals that observed variables in survival analysis cannot explain. If these unobserved heterogeneities are not accounted for, the results will be biased [14]. In mixture cure models, random effects can be used to compensate for the lack of information on other important covariates and controlling some individual characteristics that effect on survival time. Such models are referred to as mixture cure frailty models [15].
Many studies have been conducted on the mixture cure frailty models [16,17,18]. For example, Calsavara introduced the longterm frailty model using a nonproportional risk model and applied this model to melanoma datasets [19]. Other studies have used the mixture cure model for competing risks [20,21,22]. Ghavami used the parametric mixture cure gamma frailty model and vertical modeling of competing risks for breast cancer data in her study [23].
To our knowledge, no study in Iran has performed the factors affecting the short and longterm survival time of HIV patients using mixture cure models. For this reason, in the present work, to estimate the probability of longterm survival of each HIVinfected individual and the factors affecting the mortality of HIVinfected individuals, a semiparametric mixture cure model was applied. Also, to compensate for the lack of information on other important covariates, a mixture frailty model was used.
Materials and methods
Data
In this study, the data of 2170 HIVinfected people referring to the Kermanshah Province Disease Counseling Centers from 1998 to 2019 were examined. The study was a retrospective cohort study. The information on baseline demographic and clinical variables, including diagnosis age, gender, history of imprisonment, history of addiction, education, marital status, occupation status, ART, modes of HIV transmission, and tuberculosis infection was collected from medical records. We included those people who were diagnosed as HIVâ€‰+â€‰/ADIS based on two positive Elisa tests and one positive Western blot test or a rapid HIV test and two positive Elisa tests (that one the fourth generation of HIV ELISA test), which are used now simultaneously irrespective of nationality, age, and gender. We had no exclusion criteria. The outcome of the interest was time to diagnosis until death (survival time). All the patients who were alive at the end of the study were censored for death.
Semiparametric PH mixture cure model
The mixture cure model is a specific type of survival model. In these models, the studied population is assumed to be a mixture of uncured, susceptible, individuals who may experience the event of interest and cured individuals, nonsusceptible, who will never experience the event. The advantage of using this model over conventional survival models is that the estimated probability of longterm survival and its effective factors can also be analyzed in this model simultaneously.
Let T denote the failure time of interest, \(1  \pi ({\mathbf{Z}})\) be the probability of a patient.being cured and \(S(t{\mathbf{X}})\) be the survival probability of an uncuredpatient. Where, X and Z are observed values of two covariate vectors that may affect the survival function and probability of the uncured, respectively. Boag suggests a mixture cure model \(S_{pop} (t{\mathbf{X}},{\mathbf{Z}})\), in which \(S(t{\mathbf{Z}}) = 1\) for the cured group in the population [24]. A mixture function was used as follows:
Usually, \(\pi ({\mathbf{Z}})\) is referred to as â€œincidence" and \(S(t{\mathbf{X}})\) is â€œlatency". The proportional hazard (PH) model is used to model the latency part, the cure model is called the PH mixture.cure model. A logit link function \(\pi ({\mathbf{Z}}) = \frac{{\exp ({\mathbf{bZ}})}}{{1 + \exp ({\mathbf{bZ}})}}\), where b is a vector of unknown parameters, is used to model the effects of Z covariate vector. Proportional hazard models describe the event time distributions given to uncured patients. The proportional hazards model is defined as \(h(t) = h_{0} (t)\exp ({\mathbf{\beta X}})\), where \(h_{0} (t)\) is the arbitrary baseline hazard at time t, where \({{\varvec{\upbeta}}}\) is the regression coefficient of the covariate vector X.
The EM algorithm was used to estimate the unknown parameters in the mixture cure model [20]. To estimate the mixture cure model parameters, the followup time should be sufficiently long and the sample size should be large [25].
The Kaplan Meier curve and the Maller and Zhou test were used to test the hypothesis of curing (presence of patients with longterm survival time). If there were cured patients, the Kaplanâ€“Meier curve becomes horizontal before reaching zero [26, 27].
Mixture cure frailty model
Suppose there is a latent binary variable \(v_{i}\), \(v_{i} = 0\) if subject i is cured, and \(v_{i} = 1\) if subject i is not cured. We denote the probability of not being cured by \(\pi_{i}\) and define it as follows:
where g (Â·) is a logit link function, \({\mathbf{Z}}\) is a vector of covariates, and \({{\varvec{\uptheta}}}\) is the regression coefficient vector. If person i is uncured \(v_{i} = 1\), the hazard function at time t is defined as follows.
where \({{\varvec{\upbeta}}}\) is the regression coefficient of the covariate X, \(w_{i}\) is the subjectspecific frailty and follows a Gamma distribution, \(\Gamma (\Psi^{  1} ,\Psi )\), where \(\Psi\) is the variance of the frailty distribution. The probability density function for \(w_{i}\) is
In addition, \(h_{0} (t)\) is the arbitrary baseline hazard function. Let \(\Omega_{i} = \{ {\mathbf{X}},{\mathbf{Z}},t_{i} ,\delta_{i} :i = 1,...,n\}\) denote the observed data in subject i. The complete data likelihood for subject i can be written as:
where
and
The complete data likelihood of all the study subjects is \(L_{c} = \Pi_{i = 1}^{n} L_{i,c}\). It was used to estimate the unknown parameters of the EM algorithm for this model [28].
The probability of cure for each person was calculated using the formula \(1\widehat{p}=1(\mathrm{exp}({\varvec{Z}}\widehat{{\varvec{\theta}}})/(1+\mathrm{exp}(Z\widehat{\theta })))\), with the coefficients resulting from the fitted model and putting each person's characteristics into the formula.
The missing data were imputed using the multiple regression method [29, 30]. Also, in order to check collinearity between the variables, the variance inflation factorÂ (VIF) was used [31]. We used concordance criteria (Kindex) to evaluate the goodness of fit of the two models.
In this study, R 4.1. software (the smcure package) and STATA 16.0 (the strmcure command) were used for analysis [28, 32]. We used concordance criteria in the eva_cure package to check the Kindex [33]. The significance level was considered 0.05.
Results
The mean and median followup times were 173.75 and 119.13Â months, respectively. About 50% of the patients had more than 33Â years old. Of these, 57.5% were censored. Among patients, 83.5% were men, 52.1% were married, and 58.0% were unemployed. In addition, 97.1% had a low level of education, and 62.6% had a history of imprisonment. The frequency of drug users was 74.2%. Also, 54.3% of people did not receive ART. About 67.1% were infected through HIVinjecting drug users. The number of infected with tuberculosis was 67.9% were. Further information is provided in Table 1.
As shown in Fig.Â 1, the Kaplan Meier curve becomes flat above zero. The Kaplanâ€“Meier plot shows a clear plateau of about 141Â months for death which could justify the use of the cure models. Also, the results of Maller and Zhou's test show that there are patients with longterm survival in the present study. Table 2 shows the results of the semiparametric mixture cure model and mixture cure frailty model in the presence of covariates. The Kindex value of the frailty model was better than that of the mixture cure model and was 0.65 against 0.62. The results of the mixture cure frailty model showed that variables of ART, tuberculosis infection, history of imprisonment and mode of HIV transmission affected death (pvalueâ€‰<â€‰0.05). The variables of age, ART, education, gender, marital status, history of imprisonment and mode of HIV transmission are significant in the probability of being cured of death (pvalueâ€‰<â€‰0.05). The results of this model also show that the frailty component is significant (pvalueâ€‰=â€‰0.001). It means that there are some other important variables and some individual characteristics that were not considered in the present study.
According to the results of the mixture cure frailty model in the first part of Table 2, the cure odds in men were 3.367 times higher than in women (Pâ€‰=â€‰0.001). On the other hand, antiretroviral therapy significantly reduced the risk of death compared to those who did not receive antiretroviral therapy (ORâ€‰=â€‰0.131; Pâ€‰=â€‰0.001). The cure probability was significantly higher among higheducated people than among loweducated individuals (ORâ€‰=â€‰0.289; Pâ€‰=â€‰0.013). The cure probability was much higher in prisoners than in nonprisoners (ORâ€‰=â€‰0.351; Pâ€‰=â€‰0.001). The cure odds for persons over 33Â years old is 1.645 times higher than for those under 33Â years of age (ORâ€‰=â€‰1.645; pâ€‰=â€‰0.001). Based on the second part of Table 2, ART reduced the risk of death by almost 50% (HRâ€‰=â€‰0.458; Pâ€‰=â€‰0.001), while coinfection with tuberculosis increased the risk by about 40% (HRâ€‰=â€‰1.391; Pâ€‰=â€‰0.001). Prisoners were at lower risk of death than nonprisoners (HRâ€‰=â€‰0.628; Pâ€‰=â€‰0.001). The risk of death was more than two times among injection drug users than in other modes of transmission. As shown in Fig.Â 2, People who received ART had higher survival than those who did not. People who had a history of imprisonment were more likely to survive longer. Also, people infected with HIV through drug injection use had higher survival than those infected through other ways.
Discussion
Various studies have been performed in Iran and the world on the survival of HIVinfected patients. The conventional Cox method have used in some studies [34,35,36,37] and few works in some parts of the world have performed mixture cure models [8, 11]. In some investigations, the parameters estimation method was the EM algorithm, but in Davies's study [38], the stochastic EM algorithm was used to estimate the parameters. In the studies that have used the EM and stochastic EM algorithm, the goodness of fit of the expressed models is graphical methods. In some studies, the Markov Chain Monte Carlo sampling algorithm has been employed for model fitting and comparison, and the deviance information criteria (DIC) and the logarithm of pseudomarginal likelihood (LPML) indices have been used for evaluating the goodness of fit [39,40,41,42]. However, in the present study, although the parameters estimation method was the EM algorithm, the goodness of fit index was a concordance index, which is a more common index for survival studies.
The presence of shortterm and longterm survivors in the mixture cure models and promotion time cure models, may lead to a violation of the PH assumption [43]. When the PH is not established in these models, the Schoenfeld residual cannot be used to test for this feature [44]. So in this study, we checked the PH assumption for all covariates using the ln(ln(survival)). Examining this assumption showed that the PH is not valid for the TB variable. Some authors have pointed out that the mixture cure model can be used if there is evidence of a cure fraction in the study, even if the PH assumption is not met; while the Cox PH and the promotion time models should not be used [43, 45]. Therefore, in the present study, despite the nonestablishment of the PH assumption for TB, a mixture cure model was used.
Zhang presented a suitable alternative for modeling survival data with cure fraction and no proportional treatment effects as the accelerated hazard mixture cure model [46]. Another approach to nonproportional hazard alternatives is a lagtime model proposed by Zhang and Quan [47]. Also, Liu et al. suggested that generalized accelerated hazard cure models can be applied when the PH assumption is violated [48].
Beretta introduces a package that implements an extension of the semiparametric proportional hazards cure model with timevarying covariates and a variable selection technique based on SCAD penalized likelihood [49].
In accordance with the above points, in future study, it is suggested the proposed models apply to these data to check whether there is a significant difference between the results of the present study or not.
We tried to follow both the mixture cure frailty model and semiparametric PH mixture cure model to check which model was better. For the mixture cure frailty model, the goodness of fit index (Kindex) value was 0.65, but in the semiparametric PH mixture cure model, this value was equal to 0.62. So, this index showed that the mixture cure frailty model was better. Therefore, we reported the results based on the mixture cure frailty model. In the mixture cure models, the population was classified into two subpopulations composed: patients with longterm survival time, the probability of whom was estimated using a logistic model, and the patients with shortterm survival that were estimated using a conventional Cox model. The applied mixture cure frailty model showed that the variables of gender, prison history, age, marital status, education, mode of HIV transmission, and ART were statistically significant in longterm survival time. Also, variables of prison history, mode of HIV transmission, ART, and tuberculosis affected the â€‹shortterm survival time.
According to the results, the risk of death was higher for older ages, consistent with some studies and inconsistent with others [34,35,36,37]. Men were also at higher risk of death than women, consistent with other studies, probably because most people living with HIV were men [34, 36, 50]. In the present study, people who did not receive ART treatment had lower survival than those who received treatment, which was consistent with the results of other studies [51,52,53]. Many studies have also reported a significant positive association between HIV transmission through injecting drug use and death, which is in line with the present study [54, 55]. In the present study, in line with Rubaihayoâ€™s study [50], people with lower levels of education were more at risk of death. Tabarsi et al. [56] examined the effect of different covariates on the mortality of people living with HIV/TB. They found that tuberculosis infection had a significant effect on death, consistent with the present study. Also, in their study, unlike the present work, no significant relationship existed between prison history and death. In the present study, in line with other studies, having a prison history reduced mortality and increased the lifetime of HIV patients [57]. As Dolan [58] noted in his study, healthcare facilities were provided for prisoners without any cost, which could be an opportunity to prevent and treat HIV. In addition, nonprisoners are more prone to highrisk behaviors such as injecting drug use, having abnormal sex, etc., which increase the risk of death. In the present study, in agreement with Kposowa [59], a significant association was found between marital status and the probability of death, as single individuals had a higher risk of death than married individuals.
The total number of variables was 10. Since the interest of health policymakers was to evaluate the effects of all variables on the longterm survival of HIV patients, we did not use any variable selection in this study, while Sun applied variable selection under a nonmixturecure model for the prostate cancer [60].
One of the limitations of this study was due to using retrospectively collected data by health centers; it was not possible to examine the impact of other influential variables. Our suggestion for future work is to perform prospective cohort studies using the applied model in order to evaluate the survival of HIV patients. Moreover, a prospective study could be performed in Iran to examine the effect of genetic factors [61,62,63] and behavioral variables on the long and shortterm survival time of HIV patients.
Conclusion
This study showed that the frailty mixture cure models were more appropriate in a situation where the study population was composed of two groups of susceptible and nonsusceptible to the event of death. Also, using mixture cure models, the variables affecting short and longterm survival time can be identified. The results of this study show that people with a prison history, people who received ART treatment, and people who contracted HIV through injection drug use survive longer. Health professionals should pay more attention to these findings in the field of HIV prevention and treatment.
Availability of data and materials
The datasets analyzed during the current study are not publicly available due [Because the data is related to HIV and belongs to the ViceChancellor for Health of Kermanshah University of Medical Sciences, it must be with the permission of this ViceChancellor] but are available from the corresponding author on reasonable request.
Abbreviations
 TB:

Tuberculosis
 HIV:

Human immunodeficiency virus
 ART:

Antiretroviral therapy
 ELISA:

Enzymelinked immunosorbent assay
 AIDS:

Immunodeficiency syndrome
 EM algorithm:

Expectationâ€“maximization algorithm
 DIC:

Deviance information criteria
 LPML:

Logarithm of pseudomarginal likelihood
 VIF:

Variance inflation factor
References
World Health Organization 2020 [Available from: https://www.who.int/newsroom/factsheets/detail/hivaids.
Global HIV & AIDS statistics 2020 [Available from: https://www.unaids.org/en/resources/factsheet.
The Global HIV/AIDS Epidemic 2021 [Available from: https://www.kff.org/globalhealthpolicy/factsheet/theglobalhivaidsepidemic/.
Najafi Z, Taj L, Dadras O, Ghadimi F, Moradmand B, SeyedAlinaghi S. Epidemiology of HIV in Iran. Curr HIV Res. 2020;18(4):228â€“36.
Islamic Republic of Iran 2019 [Available from: https://www.unaids.org/en/regionscountries/countries/islamicrepublicofiran.
LeÃ£o J, Leiva V, Saulo H, Tomazella V. Incorporation of frailties into a cure rate regression model and its diagnostics and application to melanoma data. Stat Med. 2018;37(29):4421â€“40.
Guo L, Xiong Y, Joan HuX. Estimation in the Cox cure model with covariates missing not at random, with application to disease screening/prediction. Canadian Journal of Statistics. 2020;48(4):608â€“32.
Varshney M, Grover G, Ravi V, Thakur A. Cure fraction model for the estimation of longterm survivors of HIV/AIDS patients under antiretroviral therapy. J Commun Dis. 2018;5(3):1â€“10.
Taweab F, Ibrahim NA. Cure rate models: a review of recent progress with a study of changepoint cure models when cured is partially known. J Appl Sci. 2014;1(14):609â€“16.
Kim YJ. Cure rate model with bivariate interval censored data. Communications in StatisticsSimulation and Computation. 2017;46(9):7116â€“24.
Barbieri A, Legrand C. Joint longitudinal and timetoevent cure models for the assessment of being cured. Stat Methods Med Res. 2020;29(4):1256â€“70.
Wang P, Tong X, Sun J. A semiparametric regression cure model for doubly censored data. Lifetime Data Anal. 2018;24(3):492â€“508.
Kattumannil SK. Semiparametric transformation model for competing risks data with cure fraction. arXiv preprint arXiv:200702305. 2020.
Price DL, Manatunga AK. Modelling survival data with a cured fraction using frailty models. Stat Med. 2001;20(9â€“10):1515â€“27.
Vaupel JW, Manton KG, Stallard E. The impact of heterogeneity in individual frailty on the dynamics of mortality. Demography. 1979;16(3):439â€“54.
Tawiah R, McLachlan GJ, Ng SK. A bivariate joint frailty model with mixture framework for survival analysis of recurrent events with dependent censoring and cure fraction. Biometrics. 2020;76(3):753â€“66.
Peng Y, Zhang J. Estimation method of the semiparametric mixture cure gamma frailty model. Stat Med. 2008;27(25):5177â€“94.
Liu K, Balakrishnan N, He M. Generalized Birnbaumâ€“Saunders mixture cure frailty model: inferential method and an application to bone marrow transplant data. Communications in StatisticsSimulation and Computation. 2021:1â€“25.
Calsavara VF, Milani EA, Bertolli E, Tomazella V. Longterm frailty modeling using a nonproportional hazards model: Application with a melanoma dataset. Stat Methods Med Res. 2020;29(8):2100â€“18.
Cai C. Advanced methodology developments in mixture cure models. 2013.
Basu S, Tiwari RC. Breast cancer survival, competing risks and mixture cure model: a Bayesian analysis. J R Stat Soc A Stat Soc. 2010;173(2):307â€“29.
Gee C. Cure rate estimation in the analysis of survival data with competing risks: University of Southern California; 2004.
Ghavami V, Mahmoudi M, Foroushani AR, Baghishani H, Yaseri M, Shandiz FH. A competing risks cure frailty model: An application to relapsefree survival of breast cancer patients. Pakistan Journal of Statistics and Operation Research. 2021:591â€“605.
Boag JW. Maximum likelihood estimates of the proportion of patients cured by cancer therapy. J Roy Stat Soc: Ser B (Methodol). 1949;11(1):15â€“44.
Farewell VT. Mixture models in survival analysis: Are they worth the risk? Canadian Journal of Statistics. 1986;14(3):257â€“62.
Zhou S, Maller RA. The likelihood ratio test for the presence of immunes in a censored sample. Statistics. J Theoretical Appl Stat. 1995;27(1â€“2):181â€“201.
Maller RA, Zhou X. Survival analysis with longterm survivors: Wiley New York; 1996.
Xu Y, Cheung YB. Frailty models and frailtymixture models for recurrent event times. Stand Genomic Sci. 2015;15(1):135â€“54.
Stuart EA, Azur M, Frangakis C, Leaf P. Multiple imputation with large data sets: a case study of the Childrenâ€™s Mental Health Initiative. Am J Epidemiol. 2009;169(9):1133â€“9.
Klebanoff MA, Cole SR. Use of multiple imputation in the epidemiologic literature. Am J Epidemiol. 2008;168(4):355â€“7.
Johnston R, Jones K, Manley D. Confounding and collinearity in regression analysis: a cautionary tale and an alternative procedure, illustrated by studies of British voting behaviour. Qual Quant. 2018;52(4):1957â€“76.
Cai C, Zou Y, Peng Y, Zhang J. smcure: An RPackage for estimating semiparametric mixture cure models. Comput Methods Programs Biomed. 2012;108(3):1255â€“60.
Zhang Y, Shao Y. Concordance measure and discriminatory accuracy in transformation cure models. Biostatistics. 2018;19(1):14â€“26.
Rezaianzadeh A, Abbastabar H, Rajaeefard A, Ghaem H, Abdollahi M. Determinant factors of survival time in a cohort study on HIV patient using by timevarying cox model: Fars province, south of Iran. International Journal of Epidemiologic Research. 2017;4(2):145â€“55.
Mangal TD, Meireles MV, Pascom ARP, de Almeida CR, Benzaken AS, Hallett TB. Determinants of survival of people living with HIV/AIDS on antiretroviral therapy in Brazil 2006â€“2015. BMC Infect Dis. 2019;19(1):1â€“9.
Poorolajal J, Molaeipoor L, Mohraz M, Mahjub H, Ardekani MT, Mirzapour P, et al. Predictors of progression to AIDS and mortality postHIV infection: a longterm retrospective cohort study. AIDS Care. 2015;27(10):1205â€“12.
Seyoum D, Degryse JM, Kifle YG, Taye A, Tadesse M, Birlie B, et al. Risk factors for mortality among adult HIV/AIDS patients following antiretroviral therapy in Southwestern Ethiopia: an assessment through survival models. Int J Environ Res Public Health. 2017;14(3):296.
Davies K, Pal S, Siddiqua JA. Stochastic EM algorithm for generalized exponential cure rate model and an empirical study. J Appl Stat. 2021;48(12):2112â€“35.
Kannan N, Kundu D, Nair P, Tripathi R. The generalized exponential cure rate model with covariates. J Appl Stat. 2010;37(10):1625â€“36.
Swain PK, Grover G, Goel K. Mixture and nonmixture cure fraction models based on generalized gompertz distribution under Bayesian Approach. Tatra Mt Math Publ. 2016;66(1):121â€“35.
Martinez EZ, Achcar JA, JÃ¡come AA, Santos JS. Mixture and nonmixture cure fraction models based on the generalized modified Weibull distribution with an application to gastric cancer data. Comput Methods Programs Biomed. 2013;112(3):343â€“55.
Kim S, Chen MH, Dey DK. A new threshold regression model for survival data with a cure fraction. Lifetime Data Anal. 2011;17(1):101â€“22.
Legrand C, Bertrand A. Cure models in cancer clinical trials. Textbook of Clinical Trials in Oncology: Chapman and Hall/CRC; 2019. p. 465â€“92.
Amico M, Van Keilegom I. Cure models in survival analysis. Annual Review of Statistics and Its Application. 2018;5:311â€“42.
Amico M. Cure models in survival analysis: from modelling to prediction assessment of the cure fraction 2018.
Zhang J, Peng Y. Accelerated hazards mixture cure model. Lifetime Data Anal. 2009;15(4):455â€“67.
Kim HT, Gray R. Threecomponent cure rate model for nonproportional hazards alternative in the design of randomized clinical trials. Clin Trials. 2012;9(2):155â€“63.
Liu X, Xiang L. Generalized accelerated hazards mixture cure models with intervalcensored data. Comput Stat Data Anal. 2021;161: 107248.
Beretta A, Heuchenne C. penphcure: Variable selection in proportional hazards cure model with timevarying covariates. The R Journal. 2021;13(1):116.
Rubaihayo J, Tumwesigye NM, KondeLule J, Makumbi F, Nakku EJ, Wamani H, et al. Trends and predictors of mortality among HIV positive patients in the era of highly active antiretroviral therapy in Uganda. Infectious Disease Reports. 2015;7(3):66â€“73.
Samson OO, Nyandwaki M, Mutwiri P, Wandera O, Oriema O, Otumba E. The application of logistic regression in modeling of survival chances of HIVpositive patients under highly active antiretroviral therapy (HAART): A Case of Nyakach District, Kenya. Basic Research Journal of Medicine and Clinical Sciences. 2014;3(3):14â€“20.
Shahrin L, Leung DT, Matin N, Pervez MM, Azim T, Bardhan PK, et al. Characteristics and predictors of death among hospitalized HIVinfected patients in a low HIV prevalence country: Bangladesh. PLoS ONE. 2014;9(12): e113095.
Bisson GP, Bastos M, Campbell JR, Bang D, Brust JC, Isaakadis P, et al. Mortality in adults with multidrugresistant tuberculosis and HIV by antiretroviral therapy and tuberculosis drug use: an individual patient data metaanalysis. The Lancet. 2020;396(10248):402â€“11.
Qian HZ, Schumacher JE, Chen HT, Ruan YH. Injection drug use and HIV/AIDS in China: review of current situation, prevention and policy implications. Harm Reduct J. 2006;3(1):1â€“8.
Etemad K, Sharhani A, Higgs P, Gheibipour H. Prevalence of HIV infection and associated risk behaviors among people who inject drugs in Kermanshah, Iran: 2017. Journal of Substance Use. 2020;25(4):421â€“4.
Tabarsi P, Chitsaz E, Moradi A, Baghaei P, Farnia P, Marjani M, et al. Treatment outcome, mortality and their predictors among HIVassociated tuberculosis patients. Int J STD AIDS. 2012;23(9):e1â€“4.
Braithwaite RL, Arriola KR. Male prisoners and HIV prevention: a call for action ignored. Am J Public Health. 2003;93(5):759â€“63.
Dolan K, Wirtz AL, Moazen B, NdeffoMbah M, Galvani A, Kinner SA, et al. Global burden of HIV, viral hepatitis, and tuberculosis in prisoners and detainees. The Lancet. 2016;388(10049):1089â€“102.
Kposowa AJ. Marital status and HIV/AIDS mortality: evidence from the US National Longitudinal Mortality Study. Int J Infect Dis. 2013;17(10):e868â€“74.
Sun L, Li S, Wang L, Song X. Variable selection in semiparametric nonmixture cure model with intervalcensored failure time data: An application to the prostate cancer screening study. Stat Med. 2019;38(16):3026â€“39.
Laplana M, Caruz A, Pineda JA, Puig T, Fibla J. Association of BST2 gene variants with HIV disease progression underscores the role of BST2 in HIV type 1 infection. J Infect Dis. 2013;207(3):411â€“9.
Singh R, Ramsuran V, Naranbhai V, YendeZuma N, Garrett N, Mlisana K, et al. Epigenetic Regulation of BST2 Expression Levels and the Effect on HIV1 Pathogenesis. Front Immunol. 2021;12: 669241.
Dyer WB, Zaunders JJ, Yuan FF, Wang B, Learmont JC, Geczy AF, et al. Mechanisms of HIV nonprogression; robust and sustained CD4+ Tcell proliferative responses to p24 antigen correlate with control of viraemia and lack of disease progression after longterm transfusionacquired HIV1 infection. Retrovirology. 2008;5(1):1â€“14.
Acknowledgements
This study is adapted from a part of the first authorâ€™s Ph.D. We hereby thank and appreciate the Research and Technology ViceChancellor of Hamadan University of Medical Sciences for the approval and support of the study.
Funding
The study was funded by Vicechancellor for Research and Technology, Hamadan University of Medical Sciences (Grant No.9812209740).
Author information
Authors and Affiliations
Contributions
Kh. NGh., H. M., Gh. R. and J. P. conceived the research topic, explored that idea, performed the statistical analysis and drafted the manuscript. E. Sh and K. KH provided the data. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
We obtained written informed consent from all the participants; for the illiterates and participants under the age of 16, it was obtained from the parents/legally authorized representatives with confidentiality regarding the patientsâ€™ names and surnames. This study was approved by the Research Ethics Committee of Hamadan University (no. IR.UMSHA.REC.1398.1066). The study adhered to the relevant guidelines and regulations.
Consent for publication
Not applicable.
Competing interests
The authors declare no conflict of interest.
Additional information
Publisherâ€™s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Ghobadi, K.N., Roshanaei, G., Poorolajal, J. et al. The estimation of long and short term survival time and associated factors of HIV patients using mixture cure rate models. BMC Med Res Methodol 23, 123 (2023). https://doi.org/10.1186/s1287402301949x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s1287402301949x