An assessment of adverse drug reactions among HIV positive patients receiving antiretroviral treatment in South Africa
© Masenyetse et al.; licensee BioMed Central. 2015
Received: 26 August 2014
Accepted: 2 February 2015
Published: 5 March 2015
Antiretroviral treatment (ART) has been effective in reducing HIV/AIDS related morbidity and mortality. However, the use and uptake of ART has resulted in adverse reactions, due mainly to the medicine’s toxicity and interactions with other medicines. The timing of adverse drug reactions (ADRs) among these patients is a critical public health issue for antiretroviral (ARV) treatment adherence and retention. Reliable monitoring of HIV patients on ART is through a structured pharmacovigilance surveillance system. However, recurrent nature of these data pose challenges in their analyses. This study aimed at modelling the timing of ADR events in HIV patients on ART using correlated time-to-event models.
The data concern 590 HIV patients registered onto the Medunsa National ARV Pharmacovigilance Surveillance System within 6 months of ART initiation between February 2007 and July 2011. Recurrent times of ADRs and baseline characteristics: patient gender, and age, ART regimen, clinic and initiation period were extracted from the data. The recurrent ADR events data were modelled using both shared frailty and marginal models on the five patients’ characteristics as covariates.
Out of 590 patients, 67% were female, 68% started on regimen: Stavudine, Lamivudine and Efavirenz; 37% had experienced at least one ADR and 67% started ART in 2009–2011. Age (p-value = 0.0210), clinic (p-value < 0.0001) and period of ART initiation (p-value = 0.0002) were significantly associated with timing of first ADR. There was a significantly higher rates of ADR recurrences in patients aged 38–44 years [HR = 2.45; 95% CI = (1.47; 4.10)] vs. 30 years and less, patients taking regimen: Zidovudine, Lamivudine and Nevarapine) vs. regimen: Stavudine, Lamivudine and Efavirenz [HR = 2.09; 95% CI = (1.35; 3.22)], while the rate was lower among those who started ART in 2009–2011 vs. those who initiated in 2007–2008 [HR = 0.55; 95% CI = (0.40; 0.76)].
More realistic time-to-event models for recurrent events data have been used to analyse timing of ADR events in HIV patients taking ARV treatment. Age, antiretroviral regimen type and period of initiation of ART were associated with the timing of HIV/AIDS drug related adverse reactions regardless of the analysis model used. This study has public health policy implications in addressing the added morbidity among HIV patients taking ARV treatment in the context of universal scaling up of ARV treatment.
The Human Immunodeficiency Virus (HIV) has changed from life threatening to chronic condition due to the almost universal use and accessibility of antiretroviral treatment (ART) among HIV patients . Antiretroviral (ARV) treatment works by providing suppression of viral load and restoring the immune system. It is estimated that out of the 35.3 million people living with HIV worldwide, 10.6 million were receiving ART in 2012 . Nearly, 6.6 million HIV/AIDS related deaths worldwide have been prevented as a result of ART . Despite these gains, adverse reactions to these medicines remain a significant public health concern and may compromise the effectiveness of the ART programmes [3,4].
The risk of adverse drug reactions (ADRs) arises because of the effect of the disease on the immune systems and the safety profiles of the complex ART drugs . There are a number of ADRs related to ART that have been documented, and may be mild to severe; and short to long term depending on the environment [1,5-12]. ADRs in developing countries may differ from those in developed countries because of high prevalence of conditions such as malnutrition, tuberculosis and patients presenting with advanced HIV disease . For instance, it has been found that in Africa, neuropathy, neutropenia and lipodystrophy are the predominant ADRs . Short term ADRs are a potential threat to successful initiation and adherence to ART . The timing of ADRs may also depend on the type of drugs. Studies have shown that patients on Efavirenz, Lamivudine and Zidivudine or Indinavir, Zidovudine and Lamivudine may present with ADRs within the first 12 or 24 weeks, respectively [16,17]. ADRs may be common or specific to class of drugs [1,8,15]. Drugs classified as non-nucleoside reverse transcriptase inhibitors (NNRTIs) which include Efavirenz (EFZ) and Nevirapine (NVP) are known to cause rashes and hepatotoxicity. On the other hand drugs classified as nucleoside reverse transcriptase inhibitors (NRTIs) including Zidovudine (AZT) and Stavudine (d4T) are known to cause anemia, nausea, rashes, lipoatrophy and lactic acidosis .
Apart from ADR depending on the environment and the type of ART regimen, a number of other risk factors have been identified, that include patient age, gender, duration on treatment, disease biomarkers such as CD4 count and viral load and body mass index (BMI) [7,10,18,19]. These risk factors have been found to interact with type of ADR. For instance females are more likely to develop rashes and hepatotoxicity [7,18]; and patients aged 40 years and above are at a higher risk of developing peripheral neuropathy when taking d4T . The longer a patient is on ART the less likely they would experience ADRs; possibly as a result of stability in ARV regimen, coming after many changes and eventually settling on an acceptable regimen .
Monitoring safety and toxicity related to ART remains a challenge facing the public health sector. Monitoring is usually done using spontaneous surveillance of HIV patients on treatment. Spontaneous reporting of ADRs is a very inefficient system in detecting drug-related conditions, leading to underestimation of the burden due to ADRs [3,20,21]. Thus, more systematic and more robust surveillance methods including structured surveillance pharmacovigilance systems, which assesse and monitor safety profile and impact of antiretroviral medicines have been advocated . Structured surveillance tracks HIV positive patients who are on ART to assess drug related morbidity and mortality over time. South Africa, a country heavily hit by the HIV epidemic, uses spontaneous surveillance of HIV patients on ART to assess ART-related adverse effects. Though these data are routinely available, the coverage of important patient data may not be adequate. Thus, for the purposes of this study, data from a structured surveillance system in South Africa are used.
The adverse drug reaction events in patients often are of recurrent nature, such that the repetitions tend to cluster more in some patients than in others. Analyses of these data are complicated due to the fact that independence between the recurrent event times cannot be assumed in a subject. In medical studies, time-to-event models have been developed to account for possible dependence between recurrent events data [22,23]. The aim of this paper was to provide a unified analysis of recurrent ADR events data from a structured antiretroviral pharmacovigilance surveillance system. The authors are not aware of any such study that has comprehensively analysed data from this kind of structured surveillance system for HIV positive patients. A previous study investigated ADRs in a sample of adult inpatients at a local hospital in South Africa, and compared the distribution between HIV and non-HIV patients . This present study is very specific and unique by assessing ARV-related ADRs in a cohort of only HIV patients in a country with one of the highest HIV burden of the disease and with over 75% of the 6.1 million HIV patients on free ART .
Data are from the Medical University of South Africa (MEDUNSA) National ARV Pharmacovigilance Centre, which is based at the University of Limpopo MEDUNSA campus. Recruitment of patients started in January 2007 from four clinics, two in Gauteng, one in Limpopo and one in Mpumalanga provinces. Patients were followed up to obtain information on the ADRs they were experiencing and these data were collected during patients’ clinic visitations. Other information contained in the data were the type of regimens patients were taking together with the dates when they switched or collected the same regimens. A descriptive of the data generated by the system can be found in Dube et al. .
The study included HIV positive patients aged 15 years and above, receiving ART in the South Africa public health sector . Patients who met the inclusion criteria were recruited when visiting public health sector clinics using systematic random sampling. Patients who agreed to participate in the study were asked to sign consent forms. Patients who were antiretroviral naïve or on ART at private clinics were excluded from the study. Further details about the study can be found in Dube et al. . For purposes of this paper, data were restricted only to 590 patients enrolled within 6 months of their ART initiation between January, 2007 and August, 2011. The covariates considered in this analysis were age, gender, clinics and antiretroviral medication patients are taking.
In previous analyses of these types of data, standard logistic and survival models, assuming ADR events as independent, were fitted [20,21]. For this study ADR recurrent events in a patient could not be taken to be independent of each other: some patients could have been more prone to experience ADRs than others. Thus analyses of these data must take into account; possible dependence in the occurrence of ADRs. This study used more robust modelling techniques based on time-to-event models for correlated survival data. In particular, the marginal and shared frailty models were used to model the recurring ADR events data in order to take account of possible correlation between the data [22,23].
There are some analytical assumptions differentiating the two proposed correlated survival models. The marginal model fits separate failure times using survival models, while the dependence structure unspecified but inflates the estimated variance-covariance parameters of regression coefficients . Under the shared frailty model subject-specific random effects, which capture dependence in the failure times within a subject, are assumed to be time constant and independent and identically distributed from a known distribution function. This can be relaxed to have the subject-specific random effects to be time-varying as in Manda and Renate  and to be distributed nonparametrically as in Manda . The fitting of the models were implemented in STATA version 12 as described in Cleves and StataCorp .
Patients demographic and baseline characteristics, HIV Patients on ART, South Africa 2007-2012
Patients with adverse drug effects (%)
0 ADR → 1 st ADR
1 st ADR → 2 nd ADR
2 st ADR → 3 nd ADR
Age at clinic visit
30 and less
31 – 37
38 – 44
1a (d4T, 3TC, EFZ)*
1b (d4T, 3TC, NVP)*
1c (AZT, 3TC, EFZ)*
Year started ART
2009 and after
Adjusted Hazard Ratios (HR) of various characteristics on the occurrence of ADRs, HIV Patients on ART, South Africa 2007-2012
1 st ADR
1 st ADR → 2 nd ADR
HR (95% CI)
HR (95% CI)
HR (95% CI)
HR (95% CI)
1.13 (0.82, 1.55)
1.32 (0.72, 2.42)
1.23 (0.88, 1.72)
1.23 (0.94, 1.62)
30 and less
31 – 37
1.33 (0.81, 2.19)
2.41 (0.54, 10.69)
1.59 (0.95, 2.66)
1.51 (0.94, 2.43)
38 – 44
1.83 (1.12, 3.01)*
3.19 (0.73, 13.93)
2.45 (1.47, 4.10)*
2.07 (1.29, 3.32)*
1.43 (0.85, 2.42)
3.23 (0.73, 14.29)
2.02 (1.16, 3.51)*
1.80 (1.11, 2.93)*
0.77 (0.27, 2.20)
0.73 (0.25, 2.11)
0.78 (0.28, 2.22)
3.74 (2.54, 5.50)*
2.36 (0.89, 6.22)
4.03 (2.75, 5.91)*
3.20 (2.24, 4.57)*
1.99 (1.10, 3.58)*
2.76 (0.85, 8.99)
2.14 (1.12, 4.09)*
1.96 (1.18, 3.25)*
1a (d4T, 3TC, EFZ)
1b (d4T, 3TC, NVP)
1.01 (0.61, 1.66)
0.65 (0.19, 2.21)
1.08 (0.65, 1.80)
0.99 (0.64, 1.55)
1c (AZT, 3TC, EFZ)
1.86 (1.20, 2.87)*
1.83 (0.91, 3.69)
2.09 (1.35, 3.22)*
1.83 (1.27, 2.62)*
1.00 (0.65, 1.54)
1.18 (0.50; 2.82)
1.35 (0.90, 2.01)
1.43 (0.99, 2.04)
Year started ART
2009 and after
0.78 (0.58, 1.06)
0.41 (0.23, 0.74)*
0.55 (0.40, 0.76)*
0.90 (0.70, 1.16)
The rate of experiencing first and second ADRs were higher for females compared to males though not significant, [Hazard ratio = 1.13; 95% Confidence Interval = (0.82; 1.55)] and [HR = 1.32; 95% CI = (0.72; 2.42)] respectively. Older age groups had higher rate of experiencing ADR compared to the younger age group of 30 years and less, with significant risks being for patients in aged 38 – 44 years [HR = 1.83; 95% CI = (1.12; 3.01)]. Patients who were in regimen 1(c) had significantly higher risks of having an ADR compared to those in regimen 1(a) [HR = 1.86; 95% CI = (1.20; 2.87)]. Patients on regimen 1b had a lesser rate of experiencing recurrent ADRs compared to those in regimen 1a though this was not significant [HR = 0.65; 95% CI = (0.19; 2.21)].
The marginal model results indicate that patients in age 38 – 44 years and from 45 years had significantly higher rate of recurrence of ADRs compared to those 30 years and younger [HR = 2.45; 95% CI = (1.47; 4.10)] and [HR = 2.02; 95% CI = (1.16; 3.51)] respectively. Patients on regimen 1(c) had significantly higher hazards of recurrence of ADRs compared to patients on regimen 1(a) [HR = 2.09; 95% CI = (1.35; 3.22)]. Patients who initiated ART from 2009 had significantly lower hazards of recurrence of ADRs compared to patients who initiated ART before 2009 [HR = 0.55; 95% CI = (0.40; 0.76)].
Frailty model results show that 31 years and older patients experienced higher rates of ADRs compared to those aged 30 years and less. Patients who were initiated ART from 2009 had significantly lower risks of recurrence of ADRs compared to patients who were initiated before 2009 [HR = 0.90; 95% CI = (0.70; 1.16)], though this finding was not significant. Patients attending clinic C and D had significantly higher rate of recurrence of ADRs compared to patients in clinic A [HR = 2.07; 95% CI = (1.29; 3.32)] and [HR = 1.80; 95% CI = (1.11; 2.93)] respectively.
Timing of first ADR was also analysed using logistic regression where the outcome was categorized as 1 if a patient had an ADR and 0 otherwise. The substantive results were the same as those obtained using survival models. Using logistic regressions models this way does not account for possible correlation between recurrent ADRs, Additionally, exposure time is not taken into account in the logistics regression models, in contrast to time-to-event models described in the preceding sections.
This study has shown the utility of using correlated survival models to the analysis of recurrent adverse drug reactions events to account for possible dependence in the events. The models have been applied to a very important data based on structured surveillance of HIV patients on antiretroviral treatment in South Africa. The country has the largest number of people living with the HIV virus and on ART in the world. Age, gender, ARV regimen type, period of initiation of ART were found to be associated with the timing of HIV/AIDS drug related adverse reactions regardless of the analysis model used, though gender was not statistically significant. The highly prevalent ADRs included neuropathy, rash, cough and lipodystrophy, comparable with results found in Nwokike . As in Reddenna et al. , neuropathy was the most prevalent adverse drug reaction condition.
Patients older that ages 38 years experienced significantly higher recurrence of ADRs compared to patients aged 30 years and less, which supports the findings previous studies [7,19]. Gender differentials were also found with females having higher risks of ADRs than males. Though not statistically significant, the finding of gender differences in experiencing ARV-related ADRs has been observed in other African settings [18,19,27].
Patients taking AZT + 3TC + NVP combination had higher rates of ADRs compared to patients on d4T + 3TC + EFZ. It has also been found that patients taking d4T + 3TC + NVP experienced higher rates of ADRs compared to patients taking d4T + 3TC + EFZ. These findings confirm previous ARV-related ADR data in Spaulding et al.  and Webster et al. . Due to high rates of ARV-related ADRs, South Africa patients are now using Tenofovir (TDF) containing regimen as a first line ARV treatment. Depending on the prognosis, AZT + 3TC + NVP and then d4T + 3TC + EFZ/NVP are offered . More pharmacovigilance studies are needed to compare ADRs among patients on TDF containing regimen with patients on AZT + 3TC + NVP, a regimen found with lower ADR rates in this study. Patients who started ART after 2009 had a lesser rate of occurrence of ADRs; this may be due better care and management of HIV treatment.
This study showing high rates of ARV-related adverse drug reactions has highlighted the added morbidity among HIV patients taking ARV treatment. As indicated in Mehta et al.  in contract to other ADRs, almost all ARV-related ADRs are often inevitable and unpredictable, which makes treatment of these ADRs problematic. This creates extra cost burden to most public health systems that are stretched over many important health problems. In terms of policy implication, a better understanding, timely and proactive pharmacovigilance surveillance and reporting, especially of problematic regimen and patients subgroups, of ARV-related ADRs is advisable .
Limitations of the study
The findings from this study should be interpreted within caution. The surveillance data that was used did not have grading of the severity of ADRs. Some studies have graded the ADRs using the World Health Organization grading or the Hartwig scale to determine the level of severity to patients’ morbidity [10,19,31,32]. Not all the data in the database were used as some patients were registered into the surveillance system study long after they had already initiated on ART. Thus, information on ADRs prior to being enrolled into the study was missing; using these patients in the analyses would have biased timing of first ADR. Other relevant covariates including CD4 count, viral load and body mass index were not used as these data had a high rate of missing values.
Further limitation of the analyses in this study was the lack of adherence data, which has been found to affect the rate of ADRs due to inadequate treatment uptake . There were very few patients in the study who were started on ARTs based on the 2010 South African guidelines. These patients were included under “others” category, therefore the Tenofovir (TDF) containing combination was not individually assessed in the analysis. Lastly, several factors that may predispose patients to adverse reactions of antiretroviral medications including non-antiretroviral co-administered medication and treatment; alcoholism and viral hepatitis co-infection were also not explored [34-36].
Regarding statistical models for recurrent events data analysis, other models that could have been used include Anderson and Gill  and conditional models of Prentice, William and Peterson  instead of the WLW marginal model. The specification of the frailty model that was used for analysis in this paper could be improved by included time dependant random effects . In terms of missing data, imputation methods could have been used on the variables with missing values and included in the analysis.
This study has shown differential effects of patient’s age and antiretroviral regimen on the risks of adverse drug reactions in HIV-infection populations. It has highlighted the important of optimal drug selection and monitoring of vulnerable patients in mitigating the clinical severity of the adverse reactions. Further analyses using much larger and more complete data from different settings would be needed to solidify the findings in this study. Using correlated survival models is encouraged; however most of these data are observational in nature, which makes ascertaining casualty and associations problematic. Thus, to fully understand the determinants and dynamics of ADRs in HIV patients on ART, statistical methods including matching, propensity scores and instrumental variables techniques and where it is ethically possible, pharmacovigilance clinical trials should be encouraged. ADRs in HIV patients on ARV treatment are a huge public health problem in the free-delivery of ARV, especially in the high HIV epidemic countries of the sub-Saharan, where health resources are very limited.
The project proposal was submitted to and approved by the Medunsa Research and Ethics Committee at the University of Limpopo in 2006 (Project number MP119/2006).
We acknowledge members of the MEDUNSA National ARV Pharmacovigilance Centre for providing the data.
The original study on which the present article is based received financial support from the South African Field Epidemiology and Laboratory Training Programme (SAFELTP) funded by the Centers for Disease Control and Prevention (CDC). The Foundation for Professional Development (FPD) contributed to the development and maintenance of the database. SAFELTP, CDC and FPD had no role in the analysis, interpretation and preparation of the present article.
- Hawkins T. Understanding and managing the adverse effects of antiretroviral therapy. Antivir Res. 2010;85:201–9.View ArticlePubMedGoogle Scholar
- Joint United Nations Programme on HIV/AIDS (UNAIDS). Global report: UNAIDS report on the global AIDS epidemic 2013. Available at: http://www.unaids.org/sites/default/files/en/media/unaids/contentassets/documents/epidemiology/2013/gr2013/UNAIDS_Global_Report_2013_en.pdf.
- Mehta U. Pharmacovigilance: the devastating consequences of not thinking about adverse drug reactions. Contin Med Educ. 2011;29(6):247–2251.Google Scholar
- World Health Organization. The safety of medicines in public health programmes: pharmacovigilance an essential tool. 2006. Available at http://www.who.int/medicines/areas/quality_safety/safety_efficacy/Pharmacovigilance_B.pdf.
- Shubber Z, Calmy A, Andrieux-Meyer I, Vitoria M, Renaud_Thery F, Shaffer N, et al. Adverse events associated with nevirapine and efavirenz-based first-line antiretroviral therapy: a systematic review and meta-analysis. AIDS. 2013;27:1403–12.View ArticlePubMedGoogle Scholar
- Reddenna L, Basha SA, Gopal DV, Krishna TR. Highly Active antiretroviral Therapy: Incidence of Adverse Drug Reactions. International Journal of Allied Medicinal Sciences and Clinical Research. 2013 (1): 25–30.Google Scholar
- Srikanth AB, Babu SC, Yadav HN, Sk J. Incidence of adverse drug reactions in human immune deficiency virus-positive patients using highly active antiretroviral therapy. J Adv Pharm Technol Res. 2012;3:62–7.PubMed CentralPubMedGoogle Scholar
- Montessori V, Press N, Harris M, Akagi L, Montaner JSG. Adverse effects of antiretroviral therapy for HIV infection. Can Med Assoc J. 2004;170:229–38.Google Scholar
- Dube MN, Summers R, Tint K-S, Mayayise G. A pharmacovigilance study of adults on highly active antiretroviral therapy, South Africa: 2007–2011. Pan Afr Med J. 2012;11:39.PubMed CentralPubMedGoogle Scholar
- Nemaura T, Dhoro M, Nhachi C, Kadzirange G, Chonzi P, Masemirembwa C. Evaluation of the Prevalence, Progression and Severity of Common Adverse Reactions (Lipodystrophy, CNS, Peripheral Neuropathy, and Hypersensitivity Reactions) Associated with Anti-Retrovirial Therapy (ART) and Anti-Tuberculosis Treatment in Outpatients in Zimbabwe. J AIDS Clin Res. 2013;4:203.View ArticleGoogle Scholar
- Breckenridge A. Pharmacology of drugs for HIV. Medicine. 2009;37:374–7.View ArticleGoogle Scholar
- Meintjes G, Maartens G, Boulle A, Conradie F, Goemaere E, Hefer E, et al. Guidelines for antiretroviral therapy in adults. South Afr J HIV Med. 2012;13:114–33.Google Scholar
- Subbaraman R, Chaguturu SK, Mayer KH, Flanigan TP, Kurarasamy N. Adverse effects of Highly Active Antiretroviral Therapy in developing countries. Clin Infect Dis. 2007;45:1093–101.View ArticlePubMedGoogle Scholar
- Nwokike J. Monitoring Adverse Drug Reactions in public health programs: the case of the Nigeria TB program. Submitted to the U.S. Agency for International Development by the TBCAP Project. 2008. Available at: http://apps.who.int/medicinedocs/documents/s18400en/s18400en.pdf.
- Max B, Sherer R. Management of the adverse effects of antiretroviral therapy and medication adherence. Clin Infect Dis. 2000;30 Suppl 2:S96–116.View ArticlePubMedGoogle Scholar
- Eron JJ, Benoit SL, Jemsek J, MacArthur RD, Santana J, Quinn JB, et al. Treatment with Lamivudine, Zidovudine, or Both in HIV positive patients with 200 to 500 CD4 cells per cubic millimeter. N Engl J Med. 1995;333:1662–9.View ArticlePubMedGoogle Scholar
- Staszewski S, Morales-Ramirez J, Tashima KT, Rachlis A, Skiest D, Stanford J, et al. Efavirenz plus Zidovudine and Lamivudine, Efavirenz plus Indinavir, Indinavir plus Zidovudine and lamivudine. N Engl J Med. 1999;341:1865–73.View ArticlePubMedGoogle Scholar
- Luma NH, Doualla M-S, Choukem S-P, Temfack E, Ashuntantang G, Joko HA, et al. Adverse drug reactions of Highly Active Antiretroviral Therapy (HAART) in HIV infected patients at the General Hospital, Douala, Cameroon: a cross sectional study. Pan Afr Med J. 2012;12:87.Google Scholar
- Eluwa GI, Badru T, Akpoigbe KJ. Adverse drug reactions to antiretroviral therapy (ARVs): incidence, type and risk factors in Nigeria. BMC Cen Clin Pharmacol. 2012;12:7.View ArticleGoogle Scholar
- Agu KA, Isah MA, Oqua D, Habeeb MA, Agada PO, Ohiaeri SI, et al. Incidence of Adverse Drug Reactions in Patients on Antiretroviral Therapy: A study of Pharmaceutical Care in HIV Interventions in Nigeria. West Afr J Pharm. 2013;24:30–42.Google Scholar
- Mehta U, Durrheim DN, Blockman M, Kredo T, Gounden R, Barnes KI. Adverse drug reactions in adult medical inpatients in a South African hospital serving a community with a high HIV/AIDS prevalence: prospective observational study. Br J Clin Pharmacol. 2008;65:396–406.View ArticlePubMed CentralPubMedGoogle Scholar
- Manda SOM, Meyer R. Bayesian inference for recurrent events data using time-dependent frailty. Stat Med. 2005;24:1263–74.View ArticlePubMedGoogle Scholar
- Wei LJ, Lin DY, Weissfeld L. Regression analysis of multivariate failure time data by modelling marginal distributions. J Am Stat Assoc. 1989;48:1065–73.View ArticleGoogle Scholar
- Joint United Nations Programme on HIV/AIDS (UNAIDS). Global AIDS response progress reporting 2013: Construction of core indicators for monitoring the 2011 UN Political Declaration on HIV/AIDS.2013. Available at: http://www.unaids.org/sites/default/files/media_asset/GARPR_2013_guidelines_en_0.pdf.
- Manda SOM. A nonparametric frailty model for clustered survival data. Commun Stat: Methods Theory. 2011;40:863–75.View ArticleGoogle Scholar
- Cleves M. Analysis of multiple failure-time data with Stata. Stata Technical Bulletin 1999;49:30–39.Google Scholar
- Ofotokun I, Pomeroy C. Sex differences in adverse reactions to antiretroviral drugs. Top HIV Med. 2003;11:55–9.PubMedGoogle Scholar
- Spaulding A, Rutherford GW, Siegfried N. Stavudine or zidovudine in three-drug combination therapy for initial treatment of HIV infection in antiretroviral-naïve individuals. The Cochrane Database Syst Rev. 2010 (8): CD008651.Google Scholar
- Wester WC, Thomas AM, Bussmann H, Moyo S, Makhema JM, Gaolath T, et al. Non-Nucleoside Reverse Transcriptase Inhibitor Outcomes Among cART-Treated Adults in Botswana. AIDS. 2010;24 Suppl 1:S27–36.View ArticlePubMed CentralPubMedGoogle Scholar
- Department Health: Republic of South Africa. The South African Antiretroviral Treatment Guidelines. 2013. Available at: http://www.sahivsoc.org/upload/documents/2013%20ART%20Guidelines-Short%20Combined%20FINAL%20draft%20guidelines%2014%20March%202013.pdf.
- Anwikar SR, Bandekar MS, Smrati B, Pazare AP, Tatke PA, Kshirsagar NA. HAART induced adverse drug reactions: a retrospective analysis at a tertiary referral health care center in India. Int J Risk Saf Med. 2011;23:163–9.PubMedGoogle Scholar
- Lorio M, Colasanti J, Moreira S, Gutierrez G, Quant C. Adverse Drug Reactions to Antiretroviral Therapy in HIV-Infected Patients at the Largest Public Hospital in Nicaragua. J Int Assoc Provid AIDS Care. 2014;13:466–70.View ArticlePubMedGoogle Scholar
- Mitiku H, Abdosh T, Teklemariam Z. Factors Affecting Adherence to Antiretroviral Treatment in Harari National Regional State, Eastern Ethiopia. International Scholarly Research Notices AIDS. 2013. Article ID 960954, 7 pages, 2013. doi:10:1155/2013/960954Google Scholar
- Este JA, Cihlar T. Current status and challenges of antiretroviral research and therapy. Antivir Res. 2010;85:25–33.View ArticlePubMedGoogle Scholar
- Mallal S, Phillips E, Carosi G, Molina J-M, Workman C, Tomazic J, et al. HLA-B*5701 Screening for Hypersensitivity to Abacavir. N Engl J Med. 2008;358:568–79.View ArticlePubMedGoogle Scholar
- Dieterich DT, Robinson PA, Love J, Stern JO. Drug-induced liver injury associated with the use of nonnucleoside reverse-transcriptase inhibitors. Clin Infect Dis. 2004;38 Suppl 2:S80–9.View ArticlePubMedGoogle Scholar
- Anderson PK, Gill RD. Cox’s regression model for counting processes: a large sample study. Ann Stat. 1981;10:1100–20.View ArticleGoogle Scholar
- Prentice RL, Williams BJ, Peterson AV. On the regression analysis of multivariate failure time data. Biometrika. 1981;68:373–9.View ArticleGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.