Advanced search
Publication Cover
Global Health Action Volume 7, 2014 - Issue 1
Submit an article Journal homepage
712
Views
12
CrossRef citations to date
0
Altmetric
Original Articles

Factors associated with mortality in HIV-infected people in rural and urban South Africa

Kennedy N. OtwombePerinatal HIV Research UnitFaculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa;School of Public HealthFaculty of Health Sciences, University of the Witwatersrand, Johannesburg, South AfricaCorrespondence[email protected]
,
Max PetzoldCentre for Applied Biostatistics, Occupational and Environmental Medicine, University of Gothenburg, Gothenburg, Sweden
,
Tebogo ModisenyanePerinatal HIV Research UnitFaculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa
,
Neil A. MartinsonPerinatal HIV Research UnitFaculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa;Center for TB ResearchJohns Hopkins University, Baltimore, MD, USA
&
Tobias ChirwaSchool of Public HealthFaculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa
Article: 25488 | Received 18 Jul 2014, Accepted 03 Sep 2014, Published online: 29 Sep 2014

Abstract

Background

Factors associated with mortality in HIV-infected people in sub-Saharan Africa are widely reported. However rural–urban disparities and their association with all-cause mortality remain unclear. Furthermore, commonly used classical Cox regression ignores unmeasured variables and frailty.

Objective

To incorporate frailty in assessing factors associated with mortality in HIV-infected people in rural and urban South Africa.

Design

Using data from a prospective cohort following 6,690 HIV-infected participants from Soweto (urban) and Mpumalanga (rural) enrolled from 2003 to 2010; covariates of mortality were assessed by the integrated nested Laplace approximation method.

Results

We enrolled 2,221 (33%) rural and 4,469 (67%) urban participants of whom 1,555 (70%) and 3,480 (78%) were females respectively. Median age (IQR) was 36.4 (31.0–44.1) in rural and 32.7 (28.2–38.1) in the urban participants. The mortality rate per 100 person-years was 11 (9.7–12.5) and 4 (3.6–4.5) in the rural and urban participants, respectively. Compared to those not on HAART, rural participants had a reduced risk of mortality if on HAART for 6–12 (HR: 0.20, 95% CI: 0.10–0.39) and >12 months (HR: 0.10, 95% CI: 0.05–0.18). Relative to those not on HAART, urban participants had a lower risk if on HAART >12 months (HR: 0.35, 95% CI: 0.27–0.46).

The frailty variance was significant and >1 in rural participants indicating more heterogeneity. Similarly it was significant but <1 in the urban participants indicating less heterogeneity.

Conclusion

The frailty model findings suggest an elevated risk of mortality in rural participants relative to the urban participants potentially due to unmeasured variables that could be biological, socio–economic, or healthcare related. Use of robust methods that optimise data and account for unmeasured variables could be helpful in assessing the effect of unknown risk factors thus improving patient management and care in South Africa and elsewhere.

South Africa has the highest number of people living with HIV and the largest antiretroviral treatment (ART) programme globally (Citation1). Driven by the need to achieve the millennium development goal of increasing the number of HIV-infected people accessing treatment (Citation2), over 2 million people have been started on ART with another half million planned to be started each year. Mortality in HIV-infected individuals is high (Citation3, Citation4), especially prior to ART treatment and while awaiting ART (Citation5). This is associated with highly immunosuppressed patients with low CD4 counts and high viral loads. While survival on ART is markedly improved, mortality rates in the early stages of treatment remain relatively high (Citation4, Citation6) signifying the importance of immunosuppression at the time of ART initiation. Knowledge of risk factors for mortality is important to improve ART programmes and to prolong life. Factors associated with mortality in HIV-infected people in sub-Saharan Africa are widely reported with variations. They include CD4 and total lymphocyte counts, HIV RNA level, WHO and CDC staging, anaemia, malnutrition, and gender amongst others (Citation6–Citation11). However, prior studies do not specifically report and compare mortality in rural and urban settings. Furthermore time-to-event methods remain the most commonly used for the investigation of mortality risk factors particularly in HIV-infected cohorts (Citation12). Cox proportional hazards regression analysis, the most popular approach assumes that the risk of mortality is homogeneous (Citation13), an assumption that does not often hold (Citation14). In epidemiological studies where some participants are related either genetically or are from the same family, share some unmeasured covariates, or have multiple recurrence events, their survival times are not independent due to within -subject dependency (Citation15). Analysis of such data requires use of a method that accounts for dependence leading to heterogeneity such as frailty modelling (Citation14–Citation16), an extension of the Cox regression. Frailty modelling incorporates within-subject dependency by considering that the individual risk of events differs between subjects through accounting for unmeasured variables (Citation15–Citation17). Adjusting for within-subject correlation and accounting for unmeasured variables are some of the advantages of frailty modelling. Disadvantages include handling of competing risks and assuming the same risk within a cluster of unobserved variables.

Understanding the dynamics of mortality in HIV-infected people in sub-Saharan Africa is important as it influences patient management and care. As such, assessment of mortality surveillance patterns requires the use of robust statistical analysis methods that will optimise findings leading to appropriate policy decisions from evidence-based research. While current methods are well established, recent approaches, though more complex, provide further insight through their ability to model dependency. Detailed research from sub-Saharan Africa comparing factors associated with mortality in HIV-infected people between rural and urban areas using recently developed advanced methods are limited. Furthermore, those employing methods that account for unmeasured variables in this setting are not properly described. And, more importantly, we hypothesise that they inadequately capture uncertainty and are less precise. We therefore sought to assess and compare factors associated with mortality in HIV-infected adults in a large cohort of participants in a wellness programme in rural and urban South Africa while accounting for unmeasured variables in both sites.

Methods

Setting and study population

Data from two South African study sites were used: Perinatal HIV Research Unit (PHRU) in Soweto and Tintswalo hospital in Mpumalanga. In this study, participants enrolled in PHRU were classified as urban while those enrolled in Tintswalo were classified as rural. While both sites provided treatment using the same guideline from the department of health and similar resources such as nurses with similar qualifications, Tintswalo hospital is a primary health centre while the PHRU clinic is a research centre.

Soweto is an urban area in the city of Johannesburg in Gauteng, South Africa, with a population greater than a million people (Citation18). It is in the southwestern part of Johannesburg where approximately 40% of the population resides. While many households are poor with high unemployment rates, there is a mix of poor and wealthy residents (Citation18).

Tintswalo hospital is situated in the Bushbuckridge district of rural South Africa's northeast with a population of 541,247 people. Bushbuckridge is among the poorest areas of South Africa with up to 75% of the population living in poverty (Citation19).

A nurse-based wellness programme was established in PHRU and Tinstwalo to provide pre- and post-treatment care for HIV-infected adults between the years 2003 and 2010. Participants in the study were referred from voluntary counselling and testing centres, surrounding clinics, hospital wards, or research programmes. Referrals for ART were made for those eligible, based on the South African ART treatment guidelines at the time (CD4 cell count <200 cells/mm3) (Citation20). Inclusion criteria were diagnosis of HIV, 18 years of age or older, receiving primary clinical care at PHRU or Tintswalo clinics and providing written informed consent. Further study procedures are presented elsewhere (Citation21, Citation22).

Study questionnaires were administered at baseline and follow-up visits that occurred every 6 months. Losses to follow-up were minimised through actively following up participants by way of telephone calls, letters, and then home visits. Weight, CD4 counts, creatinine, and haemoglobin estimations were measured at scheduled six monthly visits while other investigations were done according to clinical presentation. To be considered positive for tuberculosis (TB), subjects had to meet at least one of the following criteria (Citation21): initiation of multidrug TB therapy, presence of acid fast bacilli on microscopy or a biopsy suggestive of TB, mycobacterium culture positive for acid fast bacilli or Mycobacterium tuberculosis, hospitalisation due to TB or cause of death ascribed to TB.

Measures

For this study, the outcome measure was all-cause mortality, primarily through making contact with the next of kin. Socio-demographic measures included age, household income, number of people and rooms in a household, ever smoked and ever employed. Others were height and weight for calculating body mass index (BMI), HAART status, ever had TB and CD4 counts. BMI was categorised into three groups using the CDC classification: underweight (<18.5), normal (18.5–24.9), or overweight/obese (≥25) while CD4 count was categorised into 0–200, 201–350, 351–500, and >500 cells/mm3.

Ethics

Ethical approval was provided by the Human Research Ethics Committee of the University of the Witwatersrand.

Statistical analysis

Data were stratified by site. Continuous measures at enrolment such as age, BMI, and CD4 count were assessed descriptively by median and interquartile ranges and compared between the two groups using the Kruskal–Wallis non-parametric test. Frequencies and associated proportions were determined for categorical measures and compared using the Fishers Exact and Chi--square analysis as appropriate.

Mortality rates (95% confidence interval) per 100 person-years were determined using Poisson regression modelling by site for different lengths of time pre- and post-HAART and presented graphically. The Kaplan–Meier test was used to determine differences in surveillance patterns by site.

Factors associated with mortality by site were assessed using the integrated nested Laplace approximation (INLA) (Citation23, Citation24), where univariate and multivariate shared gamma frailty models were fitted. INLA is a Bayesian inference technique, computationally efficient due to its speed relative to existing methods and more accurate parameter estimates. For the purpose of interpretation, posterior Bayesian estimates generated from INLA were exponentiated to determine their approximate hazard ratios and 95% confidence intervals (Citation25).

INLA accounted for unmeasured variables and allowed for the comparison of frailty variances between sites. A significant frailty variance with a value >1 suggests a higher rate for the event (or shorter survival times) than would be predicted under the basic Cox model while <1 suggests a lower rate (longer survival times). Hence, failure to account for frailty may overestimate or underestimate the hazard rate.

Backward selection together with inclusion of plausible socio-demographic and clinical factors was used to determine the final multivariate model. Socio-demographic variables included gender, ever employed, and ever smoked. Clinical variables included HAART status, BMI, CD4 cell count, and ever had TB. The most recent information was used for HAART status while BMI and CD4 cell count values were considered over time.

A 95% confidence interval excluding one was determined as significant. Model fit was assessed using the deviance information criterion.

Multiple imputation (MI) under the missing-at-random mechanism was used to impute 18% (n=1,195) and 11% (n=704) of missing CD4 count and BMI data, respectively (Citation26–Citation28). Five imputed datasets reflecting different scenarios were generated and analysed following the MI procedure.

Several sensitivity analyses were performed using INLA. This included assessing factors associated with mortality overall, pre- and post-HAART. Kaplan–Meier plots of time to death by gender, HAART use, CD4 count, WHO staging, and BMI were also fitted. For participants that were lost to follow-up, their time in the study up to their last visit was included in the analysis.

All statistical analyses were performed using R and SAS Enterprise Guide 5.1 under the assumption of a two-sided test at 5% significance level.

Results

Overall, 6,690 participants were enrolled and followed up for a cumulative total of 2225.8 and 8185.9 person-years in the rural (33%) and urban (67%) sites, respectively. The median (IQR) pre- and post-HAART time in years was 0.6 (0–1.9) and 1.1 (0.5–2.6), respectively. Overall, 75% were women: rural (70%) and urban (78%). The median age at enrolment was 36.4 (IQR: 31.0–44.1) and 32.7 (IQR: 28.2–38.1) years for rural and urban sites, respectively.

Rural participants were more likely to report a household income of ≤R1,000 per month compared to those in the urban areas (72.8% vs. 48.5%; p<0.0001). The proportion not on HAART at the end of the study was significantly higher in the urban area compared to the rural area (63.2% vs. 30.0%; p<0.0001) while unemployment was significantly higher in the rural area (86.7% vs. 70.7%; p<0.0001) relative to the urban areas.

At enrolment, 54% of the rural and 27% of the urban participants had either started HAART or were already on HAART. On follow-up, 17 and 10% of rural and urban participants, respectively, initiated HAART.

The median BMI of urban participants was higher than that of the rural participants (24.9 vs. 21.0; p<0.0001). Similarly, the median CD4 count at entry of rural participants was lower than that in the urban participants (178 cells/mm3 vs. 303 cells/mm3; p<0.0001, ). The proportion on HAART for <6 months and with a most recent CD4 count <200 cells/mm3 was significantly higher in the urban compared to the rural participants (76% vs. 59%; p<0.0001). Similarly, urban participants on HAART for 6–12 months and with a most recent CD4 count <200 cells/mm3 were significantly more than in the rural participants (69% vs. 53%; p=0.0004). The proportion unemployed and with a CD4 count <200 cells/mm3 was significantly higher in the rural participants compared to the urban participants (40% vs. 24%; p<0.0001). The proportion ever with TB in the rural participants was significantly higher than in the urban participants (44.8% vs. 25.4%; p<0.0001, ).

Table 1 Participant characteristics at enrolment by site

Table 2 Distribution of HAART status, employment, smoking, and TB by site

Overall, 566 (8%) of the participants died during follow-up: 11 and 7% in rural and urban sites, respectively (p<0.0001). Of these, majority were females (63% rural vs. 69% urban; p=0.1952). Overall, the mortality rate per 100 person-years in rural sites was higher at 11 (95% CI: 9.7–12.5) compared to urban sites at 4 (95% CI: 3.6–4.5). In the first 6 months, 48% of rural participants had died while 51% of the deaths in the urban site occurred after 12 months. Mortality rates per 100 person-years while on HAART were 11 (95% CI: 8.7–14.0) and 7 (95% CI: 5.9–8.3) for rural and urban participants, respectively. When taking HAART, majority of the deaths in the rural site occurred within the first 3 months of treatment initiation (49%) whereas they occurred after 12 months in the urban site (44%). Further description of those who died pre- and post-treatment in the two sites by CD4 category and BMI is presented in . Relative to the urban participants, mortality rates in the rural participants was higher ().

Fig. 1 Survival and mortality rate plots by site. a) displays the cumulative probability of death by site in the cohort of HIV-infected people in South Africa. b) is a plot of mortality rates during follow-up by site showing the period prior to initiating HAART to post-HAART.

Fig. 1 Survival and mortality rate plots by site. a) displays the cumulative probability of death by site in the cohort of HIV-infected people in South Africa. b) is a plot of mortality rates during follow-up by site showing the period prior to initiating HAART to post-HAART.

Table 3 Characteristics of those who died

In the rural participants, being on HAART for 6–12 months (HR: 0.20, 95% CI: 0.10–0.39) and >12 months (HR: 0.10, 95% CI: 0.05–0.18) relative to no HAART was protective after adjusting for CD4 count (). However, underweight BMI compared to normal BMI was associated with a higher risk of mortality (HR: 1.74, 95% CI: 1.06–2.85).

Table 4 Factors associated with mortality in rural and urban South Africa

Being on HAART for <6 months (HR: 4.48, 95% C: 3.0–6.56) and 6–12 months (HR: 2.39, 95% CI: 1.60–3.47) compared to no HAART was associated with an elevated risk of mortality in the urban participants after adjusting for CD4 count, employment, smoking, and TB history. However, those with overweight/obese BMI (HR: 0.73, 95% CI: 0.56–0.90) had a significantly lower risk of mortality compared to those with normal BMI.

The estimated frailty variances were statistically significant for rural and urban sites; 5.65 (95% CI: 4.50–7.19) and 5.4×10−5 (95% CI: 1.5×10−5 – 8.0×10−4) respectively.

Discussion

In sub-Saharan Africa where the burden of HIV is the greatest, research on factors associated with all-cause mortality in HIV-infected cohorts has widely been done. However, research on risk factors using advanced Bayesian methods such as INLA that incorporate unmeasured variables is limited (Citation23, Citation24). Common protective factors against mortality included HAART use for >12 months and CD4 count above 200 cells/mm3 while underweight BMI was positively associated with mortality. Even with differences between the rural and urban participants, the hazard rates were often in the same direction.

While both frailty variances were significant, our findings show that the frailty variance for participants in the rural site was higher than that of the urban. This suggests important unmeasured variables relevant to participants in the rural site may not have been considered (as the variance was >1) signifying more heterogeneity (variability). In contrast, the variance for participants in the urban site was <1 signifying less heterogeneity and probably that the model accounted for most of the important variables. Available literature on covariates of mortality relies on classical statistical methods that do not account for unmeasured variables (Citation12). Potential unmeasured variables include 1) biological factors such as undiagnosed TB and/or other opportunistic infections such as cryptococcosis, 2) socio-economic factors that were not included in the data collection process, 3) healthcare-related factors, 4) adherence to HAART medication, and 5) cultural and religious beliefs. A study from Portugal assessed predictors of mortality in HIV-associated hospitalisations, and through frailty modelling, it was shown how unmeasured variables in the form of quality in healthcare in different hospitals affected mortality (Citation17). Hence, an assessment that allows for the use of a mixed methods approach may provide an opportunity for further insight of the unmeasured variables that may qualify as risk factors.

Underweight BMI was associated with mortality in both rural and urban areas with a slightly higher hazard in the urban group. Our findings concur with others from previous studies (Citation5, Citation21) (Citation29). While overweight or obese BMI was not associated with mortality in rural participants, it was associated with decreased mortality in the urban participants. Previous studies have shown that overweight or obese BMI reduces the risk of mortality (Citation21). The proportion of participants with underweight BMI was higher in the rural group while that of overweight or obese participants was higher in the urban group suggesting differences in nutrition between the two groups. It may also be that the use of HAART mediated in the lack of association between mortality and overweight or obese BMI in the rural participants.

Higher, most recent CD4 count was protective in both groups as similarly reported in previous studies (Citation30, Citation31). However at enrolment, rural participants had a lower median CD4 count that qualified them to initiate therapy as per the treatment guidelines at the time. As in many treatment programmes across sub-Saharan Africa, first time testers often present for care with low CD4 count. Previous studies have underscored the relationship between treatment initiation, low CD4 count, and early mortality (Citation32–Citation35).

Those on treatment for >6 months had a lower risk of mortality in the rural participants compared to >12 months in the urban participants. Among those on treatment for 6–12 months in the urban participants, a greater proportion had a most recent CD4 count <200 cells/mm3. Findings from prior research show an elevated risk of early mortality especially in immunocompromised patients in the early period of treatment (Citation4, Citation9) (Citation34, Citation36) (Citation37). The high mortality rate in urban participants may be due to late initiation of treatment leading to advanced HIV disease. Furthermore, the association between late initiation of treatment and low CD4 count is well established. It may also be that mortality was accelerated in this group by opportunistic infections that were not identified in time including undiagnosed and untreated TB. TB remains the most common opportunistic infection diagnosed in HIV-infected people.

Unemployed urban participants had a higher risk of mortality, a finding that concurs with previous research (Citation38–Citation41). HIV care was provided at no cost in both sites irrespective of their employment status. It maybe that mortality is affected by differential access to care in urban participants. Other studies have shown that mortality could be affected by poor social situations such as unemployment, poor housing, and social isolation. Furthermore, these factors may influence one's health negatively (Citation38, Citation39).

There was an association between ever suffering from TB and mortality in participants from the urban site. Contrary to our findings, a previous study from Durban in South Africa found no association (Citation36). TB is a common opportunistic infection in HIV-infected individuals and its association with mortality is indisputable (Citation21, Citation37) (Citation42). Since urban participants were significantly more likely to have CD4 counts <200 cells/mm3 up to 12 months on HAART, it may be that they are immune-compromised and more susceptible to opportunistic infections, such as TB, and early mortality.

Some variables such as employment were self-reported and may have been influenced by social desirability bias. Using cumulative variables such as ‘ever smoking’ or ‘ever employed’ may lead to misclassification given the fact that it is the cumulative, rather than the instantaneous. While all attempts were made to collect and report all measures, there were some missing values that were estimated by MIs which may have introduced some bias. Some deaths may have been misclassified since some participants were lost to follow-up and reasons for loss to follow-up were not systematically recorded. However, an advantage of our study was the large sample size and the use of a Bayesian frailty survival model that accounts for unmeasured factors.

Conclusion

Our findings suggest that using robust advanced statistical analysis methods may provide further insight to risk factors for mortality. Furthermore, and where appropriate, a mixture of methods involving both quantitative and qualitative approaches could be employed allowing for further understanding of the unmeasured variables. Such an approach is likely to inform HIV care managers and policy makers of important factors for which interventions could be developed or existing ones enhanced to improve patient management and care. Any such interventions should alleviate rural–urban healthcare disparities not only in South Africa but also in other countries with similar experiences. Interventions should specifically help in alleviating the challenges faced by rural communities in accessing care.

Financial support

Patient care was funded by a US PEPFAR grant through USAID South Africa (674-A-00-08-00009-00).  Analysis was supported by the Consortium for Advanced Research Training in Africa (CARTA). CARTA is jointly led by the African Population and Health Research Center and the University of the Witwatersrand and funded by the Wellcome Trust (UK) (Grant No: 087547/Z/08/Z), the Department for International Development (DfID) under the Development Partnerships in Higher Education (DelPHE), the Carnegie Corporation of New York (Grant No: B 8606), the Ford Foundation (Grant No: 1100-0399), Google.Org (Grant No: 191994), Sida (Grant No: 54100029), and MacArthur Foundation Grant No: 10-95915-000-INP.  KO received an academic scholarship from Canadian--Africa Prevention Trials Network and the Forgaty International Center (grant no. TW007370/3).

Conflict of interest and funding

The authors have not received any funding or benefits from industry or elsewhere to conduct this study.

Acknowledgements

We are grateful to all the participants who enrolled in this study and allowed their data to be used for research purpose, the clinical staff for providing care and treatment, and the Statistical and Data Management Center at the Perinatal HIV Research Unit.

Related Research Data

Risk Factors for Infectious Diseases in Urban Environments of Sub-Saharan Africa: A Systematic Review and Critical Appraisal of Evidence.
Source: MDPI AG

Linking provided by 

References

  • UNAIDS. Outlook. 2010. Available from: http://www.unaids.org/en/resources/presscentre/featurestories/2010/july/20100713outloklaunch/ [cited 27 February 2014].
  • UN. The millenium development goals report. 2013. Available from: http://www.un.org/millenniumgoals/pdf/report-2013/mdg-report-2013-english.pdf [cited 27 February 2014].
  • Oliveira I, Andersen A, Furtado A, Medina C, da Silva D, da Silva ZJ, etal. Assessment of simple risk markers for early mortality among HIV-infected patients in Guinea-Bissau: a cohort study. BMJ Open. 2012; 2: 1–10.
  • Otwombe KN, Laher F, Tutu-Gxashe T, Gray G, Mohapi L. The effect of a maturing antiretroviral program on early mortality for patients with advanced immune-suppression in Soweto, South Africa. PLoS One. 2013; 8: e81538.
  • Masiira B, Baisley K, Mayanja BN, Kazooba P, Maher D, Kaleebu P. Mortality and its predictors among antiretroviral therapy naive HIV-infected individuals with CD4 cell count (350 cells/mm3 compared to the general population: data from a population-based prospective HIV cohort in Uganda. Glob Health Action. 2014; 7: 21843. [PubMed Abstract].
  • Fetzer BC, Hosseinipour MC, Kamthuzi P, Hyde L, Bramson B, Jobarteh K, etal. Predictors for mortality and loss to follow-up among children receiving anti-retroviral therapy in Lilongwe, Malawi. Trop Med Int Health. 2009; 14: 862–9.
  • Badri M, Lawn SD, Wood R. Short-term risk of AIDS or death in people infected with HIV-1 before antiretroviral therapy in South Africa: a longitudinal study. Lancet. 2006; 368: 1254–9.
  • Erikstrup C, Kallestrup P, Zinyama R, Gomo E, Mudenge B, Gerstoft J, etal. Predictors of mortality in a cohort of HIV-1-infected adults in rural Africa. J Acquir Immune Defic Syndr. 2007; 44: 478–83.
  • Jerene D, Endale A, Hailu Y, Lindtjorn B. Predictors of early death in a cohort of Ethiopian patients treated with HAART. BMC Infect Dis. 2006; 6: 136.
  • Johannessen A, Naman E, Ngowi BJ, Sandvik L, Matee MI, Aglen HE, etal. Predictors of mortality in HIV-infected patients starting antiretroviral therapy in a rural hospital in Tanzania. BMC Infect Dis. 2008; 8: 52.
  • Kouanda S, Meda IB, Nikiema L, Tiendrebeogo S, Doulougou B, Kabore I, etal. Determinants and causes of mortality in HIV-infected patients receiving antiretroviral therapy in Burkina Faso: a five-year retrospective cohort study. AIDS Care. 2012; 24: 478–90. [PubMed Abstract].
  • Otwombe KN, Petzold M, Martinson N, Chirwa T. A review of the study designs and statistical methods used in the determination of predictors of all-cause mortality in HIV-infected cohorts: 2002–2011. PLoS One. 2014; 9: e87356.
  • Cox D. Regression models and life-tables. J Roy Stat Soc B. 1972; 43: 187–220.
  • Chuang SK, Cai T, Douglass CW, Wei LJ, Dodson TB. Frailty approach for the analysis of clustered failure time observations in dental research. J Dent Res. 2005; 84: 54–8.
  • Wienke A. Frailty models in survival analysis. 2011; Florida, USA: Chapman and Hall/CRC.
  • Vaupel JW, Manton KG, Stallard E. The impact of heterogeneity in individual frailty on the dynamics of mortality. Demography. 1979; 16: 439–54.
  • Dias SS, Andreozzi V, Martins MO, Torgal J. Predictors of mortality in HIV-associated hospitalizations in Portugal: a hierarchical survival model. BMC Health Serv Res. 2009; 9: 125.
  • City of Johannesburg. 2012/16 Integrated Development Plan. 2012. Available from: http://www.joburg-archive.co.za/2012/pdfs/idp/idp201216.pdf [cited 6 August 2014].
  • Bushbuckridge Municipality. Bushbuckridge local municipality local development strategy. 2010. Available from: http://bushbuckridge.gov.za/wp-content/uploads/2014/03/Bushbuckridge-LED-Strategy-2010–20142.pdf [cited 11 June 2014].
  • Southern African HIV Clinicians Society: antiretroviral therapy in adults. South Afr J HIV Med. 2008; 9: 18–31.
  • Hanrahan CF, Golub JE, Mohapi L, Tshabangu N, Modisenyane T, Chaisson RE, etal. Body mass index and risk of tuberculosis and death. AIDS. 2010; 24: 1501–8.
  • Venkatesh KK, de Bruyn G, Lurie MN, Mohapi L, Pronyk P, Moshabela M, etal. Decreased sexual risk behavior in the era of HAART among HIV-infected urban and rural South Africans attending primary care clinics. AIDS. 2010; 24: 2687–96.
  • Martino S, Akerkar R, Rue H. Approximate Bayesian inference for survival models. Scand J Stat. 2010; 38: 514–28.
  • Rue H, Martino S, Chopin N. Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations. J Roy Stat Soc B. 2009; 71: 319–92.
  • Fairley L, Forman D, West R, Manda S. Spatial variation in prostate cancer survival in the Northern and Yorkshire region of England using Bayesian relative survival smoothing. Br J Canc. 2008; 99: 1786–93.
  • Little RJA. Regression with missing X's: a review. J Am Stat Assoc. 1992; 87: 1227–37.
  • Rubin DB. Multiple imputation after 18+ years. J Am Stat Assoc. 1996; 91: 473–89.
  • Rubin DB, Schenker N. Multiple imputation in health-care databases: an overview and some applications. Stat Med. 1991; 10: 585–98.
  • Koethe JR, Blevins M, Nyirenda CK, Kabagambe EK, Chiasera JM, Shepherd BE, etal. Serum phosphate predicts early mortality among underweight adults starting ART in Zambia: a novel context for refeeding syndrome?. J Nutr Metabol. 2013; 2013: 545439.
  • Poka-Mayap V, Pefura-Yone EW, Kengne AP, Kuaban C. Mortality and its determinants among patients infected with HIV-1 on antiretroviral therapy in a referral centre in Yaounde, Cameroon: a retrospective cohort study. BMJ Open. 2013; 3: 1–7.
  • Tadesse K, Haile F, Hiruy N. Predictors of mortality among patients enrolled on antiretroviral therapy in Aksum Hospital, Northern Ethiopia: a retrospective cohort study. PLoS One. 2014; 9: e87392.
  • Barth RE, Tempelman HA, Moraba R, Hoepelman AIM. Long-term outcome of an HIV-treatment programme in rural Africa: viral suppression despite early mortality. AIDS Res Treat. 2011; 2011: 434375. [PubMed Abstract] [PubMed CentralFull Text].
  • Mageda K, Leyna GH, Mmbaga EJ. High initial HIV/AIDS-related mortality and -its predictors among patients on antiretroviral therapy in the Kagera Region of Tanzania: a five-year retrospective cohort study. AIDS Res Treat. 2012; 2012: 843598. [PubMed Abstract] [PubMed CentralFull Text].
  • Schmaltz CA, Santoro-Lopes G, Lourenco MC, Morgado MG, Velasque Lde S, Rolla VC. Factors impacting early mortality in tuberculosis/HIV patients: differences between subjects naive to and previously started on HAART. PLoS One. 2012; 7: e45704.
  • Wubshet M, Berhane Y, Worku A, Kebede Y, Diro E. High loss to followup and early mortality create substantial reduction in patient retention at antiretroviral treatment program in north-west Ethiopia. ISRN AIDS. 2012; 2012: 721720.
  • Ojikutu BO, Zheng H, Walensky RP, Lu Z, Losina E, Giddy J, etal. Predictors of mortality in patients initiating antiretroviral therapy in Durban, South Africa. S Afr Med J. 2008; 98: 204–8. [PubMed Abstract] [PubMed CentralFull Text].
  • Russell EC, Charalambous S, Pemba L, Churchyard GJ, Grant AD, Fielding KL. Low haemoglobin predicts early mortality among adults starting antiretroviral therapy in an HIV care programme in South Africa: a cohort study. BMC Public Health. 2010; 10: 433.
  • Burkey MD, Weiser SD, Fehmie D, Alamo-Talisuna S, Sunday P, Nannyunja J, etal. Socioeconomic determinants of mortality in HIV: evidence from a clinical cohort in Uganda. J Acquir Immune Defic Syndr. 2014; 66: 41–47. [PubMed Abstract].
  • Delpierre C, Cuzin L, Lauwers-Cances V, Datta GD, Berkman L, Lang T. Unemployment as a risk factor for AIDS and death for HIV-infected patients in the era of highly active antiretroviral therapy. Sex Transm Infect. 2008; 84: 183–6.
  • DeSilva MB, Merry SP, Fischer PR, Rohrer JE, Isichei CO, Cha SS. Youth, unemployment, and male gender predict mortality in AIDS patients started on HAART in Nigeria. AIDS Care. 2009; 21: 70–7.
  • Sartorius B, Kahn K, Collinson MA, Sartorius K, Tollman S. Dying in their prime: determinants and space-time risk of adult mortality in rural South Africa. Geospat Health. 2013; 7: 237–49. [PubMed Abstract] [PubMed CentralFull Text].
  • MacPherson P, Moshabela M, Martinson N, Pronyk P. Mortality and loss to follow-up among HAART initiators in rural South Africa. Trans Roy Soc Trop Med Hyg. 2009; 103: 588–93.
Download PDF
Supercharge Your Next Research Paper: Research and write your next paper with Jenni AI.

Related research

People also read lists articles that other readers of this article have read.

Recommended articles lists articles that we recommend and is powered by our AI driven recommendation engine.

Cited by lists all citing articles based on Crossref citations.
Articles with the Crossref icon will open in a new tab.