Forecast analysis of any opportunistic infection among HIV positive individuals on antiretroviral therapy in Uganda
© The Author(s). 2016
Received: 2 April 2016
Accepted: 5 August 2016
Published: 11 August 2016
Predicting future prevalence of any opportunistic infection (OI) among persons infected with the human immunodeficiency virus (HIV) on highly active antiretroviral therapy (HAART) in resource poor settings is important for proper planning, advocacy and resource allocation. We conducted a study to forecast 5-years prevalence of any OI among HIV-infected individuals on HAART in Uganda.
Monthly observational data collected over a 10-years period (2004–2013) by the AIDS support organization (TASO) in Uganda were used to forecast 5-years annual prevalence of any OI covering the period 2014–2018. The OIs considered include 14 AIDS-defining OIs, two non-AIDS defining OIs (malaria & geohelminths) and HIV-associated Kaposi’s sarcoma. Box-Jenkins autoregressive integrated moving average (ARIMA) forecasting methodology was used.
Between 2004 and 2013, a total of 36,133 HIV patients were enrolled on HAART of which two thirds (66 %) were female. Mean annual prevalence for any OI in 2004 was 57.6 % and in 2013 was 27.5 % (X2trend = 122, b = −0.0283, p <0.0001). ARIMA (1, 1, 1) model was the most parsimonious and best fit for the data. The forecasted mean annual prevalence of any OI was 26.1 % (95 % CI 21.1–31.0 %) in 2014 and 15.3 % (95 % CI 10.4–20.3 %) in 2018.
While the prevalence of any OI among HIV positive individuals on HAART in Uganda is expected to decrease overall, it’s unlikely that OIs will be completely eliminated in the foreseeable future. There is therefore need for continued efforts in prevention and control of opportunistic infections in all HIV/AIDS care programmes in these settings.
KeywordsHIV/AIDS Opportunistic infections Antiretroviral therapy Prevalence ARIMA Forecast TASO Uganda
Since the introduction of highly active antiretroviral therapy (HAART) in resource poor settings, there has been tremendous reduction in the incidence and overall prevalence OIs among persons living with HIV/AIDS . In 2014, UNAIDS came up with very ambitious programme code named “90-90-90” meaning 90 % of all the people living with HIV to be diagnosed, 90 % of all diagnosed to receive ART and 90 % of all on ART to have a fully suppressed viral load by 2020 . Basing on these targets and past achievements, there is global consensus that the HIV/AIDS epidemic could be halted by 2030 if the planned actions are fully implemented . These targets have been adopted by many developing countries including Uganda and form the basis for national policy and strategic planning. Uganda 5 year National HIV strategic plan (2015/16-2019/2020) (NSP) targets to reduce HIV-associated morbidity and mortality by 70 % and maintain 90 % viral suppression through a strategy of increased HAART coverage, early HAART initiation and improved retention in HIV care [4, 5].
HAART rollout in public funded HIV/AIDS care facilities in Uganda started in 2004 as part of the Global HAART roll out strategy . The Uganda government employed a public-private sector partnership strategy to accelerate the HAART rollout exercise. By the end of 2014, a cumulative total of 750,896 HIV positive patients out of the estimated 1.6 million HIV positive individuals in Uganda (50 %) had access to HAART . However, despite these impressive gains, there are still disparities in HAART usage by geographical area, sex and age  and the effect of increasing HAART coverage on frequency and patterns of OIs among HIV positive individuals over time and in different geographical areas is not well known.
Effective strategies for prevention of opportunistic infections and mortality associated with HIV/AIDS is of great importance given the high cost and uncertainty surrounding sustainability of lifelong ART in resource poor settings. There is therefore need to plan ahead of time in order to make current treatment/preventive interventions more effective. An accurate morbidity prediction of any OI among HIV-infected patients on HAART provides a good scientific basis for planning and formulating appropriate control or preventive policies.
The Box-Jenkins autoregressive integrated moving average (ARIMA) model is the most widely used method in forecasting of infectious diseases [7–14]. ARIMA fits univariate models for a time series where time dependent disturbances are allowed to follow a linear autoregressive moving average specification [15, 16]. Univariate ARIMA forecasting is a procedure for computing a forecast based on past and present values of a given series possibly augmented with a function of time (e.g. a linear trend) . However, univariate ARIMA forecasting is appropriate where there is large number of series to forecast and the time series are stationary (i.e. constant variance or rate of change). ARIMA models have been used successfully to monitor and forecast different infectious diseases including new annual HIV infections in Korea , monthly dengue fever cases in Thailand , monthly dengue incidence in Brazil , hantavirus outbreak in Chile , SARs outbreak in Singapore , hemorrhagic fever outbreak in China , monthly malaria cases in Burundi , monthly TB morbidity in China  and monthly mortality rate in under-five years old children in Iran .
However to the best of our knowledge, there has not been any published work on prediction of OIs associated with HIV/AIDS among HIV positive patients on HAART in Uganda. The purpose of the current study therefore was to model and forecast the prevalence of any OI among HIV-infected individuals on highly active antiretroviral therapy (HAART) in Uganda.
Data for this study was obtained from the AIDS Support Organization (TASO) in Uganda. TASO is the oldest and largest non-governmental organization (NGO) providing HIV/AIDS care and treatment in Uganda and sub-Saharan Africa . TASO started in 1987 and has since established 11 HIV clinics spread across Uganda serving mainly the rural population . TASO HIV/AIDS care and treatment programme offers a wide range of services include HIV testing and counseling, antiretroviral therapy, primary and secondary prophylaxis, OI diagnosis and treatment, home based care and psycho-social support.
With close to 3 decades of experience in HIV/AIDS care and treatment, TASO provided a good opportunity for assessing HAART effect on trends of opportunistic infections in real life programmatic settings in which HAART roll out is taking place in sub-Saharan Africa . Initially, access to HAART was based on the Ugandan ministry of Health and the World Health Organization (WHO) 2006 guide lines as indicated above [22, 23]. However, after 2010 access to HAART was expanded to include HIV patients with CD4 cell count ≤350 or WHO clinical stage3 or 4 irrespective of CD4 cell count as was recommended by Ugandan Ministry of Health and WHO/UNAIDS [24, 25]. By end of 2013, a total of 91,218 were actively in care of which 58,051 clients(64 %) were on HAART and 87,903(96.4 %) were on cotrimoxazole/dapson prophylaxis throughout the 11 HIV clinics across Uganda . All services are free of charge including anti-retroviral drugs (ARVs) for those who are eligible .
Sampling and sample size
Data collection procedure was according to an established protocol across all TASO HIV clinics. In brief, each client is expected to attend the clinic at least once a month. At each clinic visit, data per client is collected on a standardized case report form (CRF) detailing the client’s demographic characteristics, clinical condition, medical history, prophylaxis use, OI diagnosis and treatment given. The information was then compiled and entered into the TASO electronic data base using EPIINFO vs3 in Access format. Data on 14 AIDS-defining OIs and 2 non-AIDS defining OIs (malaria & geohelminths) plus HIV-associated Kaposi’s sarcoma were compiled and handed over for analysis.
To establish whether there is a trend, data on prevalence of any OI was summarized by month and calendar year. Then a monthly time plot was drawn. Monthly prevalence was calculated as a proportion of clients diagnosed with any OI in a given month. Mean annual prevalence was calculated as an average of the 12 monthly prevalences in a year. To examine long term trends in monthly prevalence, the moving average smoothing technique was used to filter out short-term fluctuations and random variation within the monthly series . The moving average smoothing technique achieves this by replacing each element of the time series by n neighbouring elements, where n is the width of the smoothing window. A centered moving average including three observations before and two observations after the current observation inclusive was used.
To forecast, the Box-Jenkins ARIMA methodology  was used. This method involves an interactive procedure including: model identification, evaluation and diagnostic checking before forecasting. Since ARIMA requires a stationary process, the Augmented Dickey-Fuller (ADF) unit-root test was used to test for stationarity of the monthly series. Non-stationary series were transformed by first order differencing to stabilize the variance. Autocorrelation function (ACF) and partial autocorrelation function (PACF) plots were used to determine possible values for the autoregressive and moving average orders. Akaike information criterion (AIC) and Bayesian information criterion (BIC) were used to identify the most parsimonious model. Parameter estimation was by maximum likelihood method. Diagnostic checks involved plotting ACF and PACF for autocorrelation structure and the Portmanteau test for white noise in the model residuals.
In order to evaluate the model forecast accuracy, data were split into two groups: one for the fitting process (2004–2010) and the other for validation (2011–2013). Forecast accuracy was assessed by computing the mean absolute percentage error (MAPE) . Finally, the fitted ARIMA model was used to forecast 5 year mean prevalence of any OI among HIV positive patients on HAART for the period 2014–2018. Root mean squared error (RMSE) was used to estimate lower and upper forecast limits . All analyses were conducted using Stata 13 (Stata Corp, TX) with p <0.05 considered significant.
Baseline characteristics of study participants who were started on HAART in the period between 2004 and 2013
Gender (36,133), n(%)
Median age(36,120), n(IQR)
Age (36,120), n(%)
Location (36,133), n(%)
Education (n = 34,337), n(%)
Marital status (n = 34,338), n(%)
Occupation (n = 33,454), n(%)
WHO clinical stage (n = 36, 133), n(%)
Weight (23,848), n(%)
Year of enrolment on HAART
Mean annual prevalence of any OI (2004–2013)
Mean annual prevalence
AIC and BIC model selection criteria
ARIMA (1 1 3)
ARIMA (2 1 1)
ARIMA (3, 2, 4)
ARIMA (1 1 2)
ARIMA (3 1 3)
ARIMA (2 1 2)
ARIMA (3 1 1)
ARIMA (4 1 3)
ARIMA (2 1 4)
Estimation of model parameters for any OI
Average annual prevalence rates declined from 58 % in 2004 when HAART was first rolled out to 28 % in 2013. The model estimated a decline in OI rate of −0.028 per year (95 % CI −0.034 to −0.023) or −0.0023 per month (95 % CI: −0.0028 to −0.0019).
t-test of the difference between observed and predicted mean annual any OI prevalence in the validation period (2011–2013)
Mean annual any OI prevalence
Mean annual any OI prevalence
t = 1.5421, p = 0.151
t = −0.620, p = 0.547
t = −0.035, p = 0.973
t = −0.519, p = 0.614
One-step ahead forecast of mean annual prevalence of any OI among HIV positive patients on HAART (2014–2018)
Predicted value (per 100)
95 % Confidence limits
Lower predicted limit
Upper predicted limit
ARIMA models have been widely used to predict morbidity and mortality due to infectious diseases with fairly good accuracy [7, 9, 10, 12, 17]. Early warning based on forecasts is very important for advocacy and resource allocation. For example ARIMA models have been successfully used to monitor and forecast new annual HIV infections in Korea , monthly dengue fever cases in Thailand , monthly dengue incidence in Brazil , hantavirus outbreak in Chile , SARs outbreak in Singapore , hemorrhagic fever outbreak in China , monthly malaria cases in Burundi , monthly TB morbidity in China  and monthly mortality rate in under-five years old children in Iran . Yu et al.  used observation data for the period 1985–2012 to forecast new HIV infections for the period 2013–2017. ARIMA models can take into account changing trends, periodic changes and random disturbances in the time series .
This study established that ARIMA (1,1,1) can be used to predict monthly OI morbidity among HIV positive patients on HAART in Uganda. The fitted model suggests the need for first-order differencing of the series to remove a stochastic decreasing trend, then a first order autoregressive (1) term and a first order moving average (1) term to cater for serial correlation in the data. This means that prevalence of any OI each month is directly influenced by the presence of any OI in the preceding month and the prediction errors of the current and preceding month. The model can also be used to estimate the cost of OIs prevention/treatment in the next 5 years of HAART after 2013 as well as informing policy, planning and advocacy for increased resource allocation towards HIV/AIDS care in Uganda.
In this study, the forecasted prevalence of any OI was 15.3 % in 2018 which show that though the prevalence of any OI is likely to reduce but will not vanish as was reported in previous studies elsewhere . Available reports suggest that OIs are likely to remain an important aspect of care even in the era of HAART especially in low income countries due to several factors including late HIV diagnosis, sub-optimal HAART use, poor adherence, dug resistance, poverty, poor nutrition, high exposure to infectious agents, just to mention a few [29–35]. Early warning based on forecasts is very important for future planning and rational resource allocation. There is therefore need for more efforts in prevention and control of opportunistic infections in all HIV/AIDS care programmes in these settings.
In the current study, results show that the ARIMA (1,1,1) model gives a good prediction of OIs prevalence in the period January 2014 - December 2018 (60 months) within 8 % margin of error i.e. MAPE = 8.2. The prediction error was close to that found in another previous study that used ARIMA (3,1,4) model to forecast dengue incidence in North-Eastern Thailand with a margin of error of 7 % (i.e. MAPE = 7.0) . The difference could be attribute to the fact that our predictions covered a slightly longer period of 60 months compared to 40 months in the previous study.
One of the limitation of this study was that the ARIMA model was fitted with a constant term (1st order differences) hence forcing the assumption that past trends (a reduction of 0.2 % in OI prevalence per month) will be maintained over time. This could only be possible if the current OI interventions are maintained and all other factors remain constant which may not be the case. Future studies should explore the effect of other time changing variables (HAART access policy, OI treatment/ prophylaxis, adherence, drug resistance, etc.) on the trends. Projections using alternative models such as ARIMA models without a constant term can also be explored in the future. Another limitation was that, although the model adequacy was evaluated, ARIMA models can only be used for short term forecasts and the accuracy of the predicted values gets less the longer the prediction period.
Despite the above limitations, the greatest strength of the present study is that it is the first in Uganda to model and predict out-of-sample prevalence of OIs among HIV positive patients on HAART. However, additional studies from other ART programmes in the country are needed in order to generalize the results at the national level.
While the prevalence of any OI among HIV positive individuals on HAART in Uganda is expected to decrease overall, it’s unlikely that OIs will be completely eliminated in the foreseeable future. There is urgent need to integrate control/prevention measures for all OIs in HIV/AIDS care programmes and to address all factors that may undermine HAART effect including poverty if OIs are to be eliminated in these settings. The forecast model in the current study could be the first attempt to predict morbidity due to OIs among HIV positive individuals on HAART in Uganda. However, ARIMA models are designed for short term forecasts. Future observational series should explore other predictive models that provide longer forecasts taking into account time changing variables and random disturbances in the time series.
ACF, autocorrelation function; AIC, Akaike information criteria; AIDS, acquired immuno-deficiency syndrome; AR, autoregressive; ARIMA, autoregressive integrated moving average; ART, antiretroviral therapy; ARV, antiretroviral; BIC, Bayesian information criteria; HAART, highly active antiretroviral therapy; HIV, human immunodeficiency virus; IQR, interquartile range; IRB, institutional review board; MA, moving average; MAPE, mean absolute percentage error; OI, opportunistic infection; PA F, patial autocorrelation function; SD, standard deviation; TASO, the AIDS support organisation; UNAIDS, joint United Nations programme on HIV/AIDS; WHO, World Health Organisation
The authors are grateful to the AIDS Support Organization (TASO) in Uganda for providing the data and Mountains of the Moon University for all the support towards this study.
Availability of data and materials
All relevant data generated or analyzed during this study are included in this manuscript and its supplementary information files (Additional file 1).
JR: conceived the study design, wrote the study protocol, collected the data, analyzed and drafted the manuscript; NMT: contributed to the study design, data analysis and review of the manuscript; JKL: contributed to the study design, data analysis and review of the manuscript. FM: contributed to the study design, data analysis and review of the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
The study obtained ethical approval from Makerere University School of Public Health Higher Degrees Research and Ethics committee (No. IRB00011353) and the Uganda National Council for Science and Technology (No. HS895). The above ethical committees waived off the need for consent to participate since this was secondary data analysis.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- WHO. HIV/AIDS facts sheet. Geneva: WHO; 2015.Google Scholar
- UNAIDS. 90-90-90 : an ambitious treatment target to help end the AIDS epidemic. Geneva: WHO; 2014.Google Scholar
- UNAIDS. Fast-Track: ending the AIDS epidemic by 2030. Geneva: UNAIDS; 2015.Google Scholar
- UNAIDS. Uganda HIV/AIDS Country Progress Report 2014. Geneva: WHO; 2015.Google Scholar
- Uganda AIDS Commission (UAC). National HIV and AIDS strategic plan 2015/2016- 2019/2020. Kampala: UAC; 2015.Google Scholar
- UNAIDS. Access to Antiretroviral therapy in Africa: Status report on progress towards 2015 targets. Geneva: UNAIDS; 2014.Google Scholar
- Zheng Y-L, Zhang L-P, Zhang X-L, Wang K, Zheng Y-J. Forecast model analysis for the morbidity of tuberculosis in Xinjiang, China. PLoS ONE. 2015;10(3):e0116832.View ArticlePubMedPubMed CentralGoogle Scholar
- Wu W, Guan P, Guo J, Zhou B. Comparison of GM(1,1) gray model and ARIMA model in fore- casting the incidence of hemorrhagic fever with renal syndrome. Journal of China Medical University. 2008;37:52–5.Google Scholar
- Siriwan W, Mullica J, Krisanadej J. Development of temporal modeling for prediction of dengue infection in Northeastern Thailand. Asian Pac J Trop Med. 2012;249:25–6.Google Scholar
- Liu Q, Liu X, Jiang B, Yang W. Forecasting incidence of hemorrhagic fever with renal syndrome in China using ARIMA model. BMC Infect Dis. 2011;11:218.View ArticlePubMedPubMed CentralGoogle Scholar
- Earnest A, Chen M, Ng D, Sin L. Using autoregressive integrated moving average (ARIMA) models to predict and monitor the number of beds occupied during a SARS outbreak in a tertiary hospital in Singapore. BMC Health Serv Res. 2005;5:316.View ArticleGoogle Scholar
- Luz P, Mendes B, Codeco C, Struchiner C, Galvani A. Time series analysis of dengue incidence in Rio de Janeiro. Brazil Am J Trop Med Hyg. 2008;79:933–9.PubMedGoogle Scholar
- Nsoesie E, Mekaru S, Ramakrishnan N, Marathe M, Brownstein J. Modeling to predict cases of Hantavirus pulmonary syndrome in Chile. PLoS Negl Trop Dis. 2014;8:e2779.View ArticlePubMedPubMed CentralGoogle Scholar
- Gomez-Elipe A, Otero A, van Herp M, Aguirre-Jaime A. Forecasting malaria incidence based on monthly case reports and environmental factors in Karuzi, Burundi, 1997–2003. Malar J. 2007;6(1):129.View ArticlePubMedPubMed CentralGoogle Scholar
- Stata corp. Time series analysis: introduction to autoregressive intergrated moving average (ARIMA) models. College station: Stata corp LP; 2012.Google Scholar
- Chatfield C. Time series forecasting. Boca Raton: Chapman & Hall/CRC; 2001.Google Scholar
- Yu H-K, Kim N-Y, Kim SS, Chu C, Kee M-K. Forecasting the number of human immunodeficiency virus infections in the Korean population using the autoregressive integrated moving average model. Osong Public Health Res Perspect. 2013;4(6):358e.View ArticleGoogle Scholar
- Rostami M, Jalilian A, Hamzeh B, Laghaei Z. Modeling and forecasting of the under-five mortality rate in Kermanshah province in Iran: a time series analysis. Epidemiol Health. 2015;37:e2015003.View ArticlePubMedPubMed CentralGoogle Scholar
- Calvarese M, Bame S, Bakanume B. Historical analysis of AIDS patients in Uganda using innovative community clinical service: the AIDS Support Organisation (TASO). J Humanit Soc Sci. 2007;1:1–13.Google Scholar
- TASO. TASO annual report. Kampala: TASO; 2013.Google Scholar
- TASO Services and Programmes. [https://www.tasouganda.org/]. Accessed 3 Feb 2015.
- WHO. Antiretroviral therapy for HIV infection in adults and adolescents; recommendations for a public health approach. Geneva: WHO; 2006.Google Scholar
- Uganda Ministry of Health. Antiretroviral treatment policy for Uganda. Kampala: Uganda Ministry of Health; 2003.Google Scholar
- WHO. Antiretroviral therapy for HIV infection in adults and adolescents: recommendations for a public health approach 2010 revision. Geneva: WHO; 2010.Google Scholar
- Uganda Ministry of Health. Antiretroviral treatment policy. Kampala: Uganda Ministry of Health; 2009.Google Scholar
- Box, Jenkins, Reinsel, editors. Time series analysis, forecasting and control. 3rd ed. Englewood Clifs: Prentice Hall; 1994.Google Scholar
- Box G, Jenkins G. Time series analysis, forecasting and control. San Francisco: Holden-Day; 1976.Google Scholar
- Chatfield C. The analysis of time series: an introduction. Washington: Chapman & Hall/CRC; 2004.Google Scholar
- Brooks J, Kaplan J, Holmes K, Benson C, Pau A, Masur H. HIV-associated opportunistic infections--going, going, but not gone: the continued need for prevention and treatment guidelines. Clin Infect Dis. 2009;48(5):609–11.View ArticlePubMedGoogle Scholar
- Bonnet F, Lewden C, May T, Heripret C, Jougla E, Bevilacqua S, Costagliola D, Salmon D, Chene G, Morlat P. Opportunistic infections as cause of death in HIV infected patients in the HAART era in France. Scand J Infect Dis. 2005;37:482–7.View ArticlePubMedGoogle Scholar
- Corey DM, Kim HW, Salazar R, Illescas R, Villena J, Gutierrez L, Sanchez J, Tabet SR. Brief report: effectiveness of combination antiretroviral therapy on survival and opportunistic infections in a developing world setting: an observational cohort study. J Acquir Immune Defic Syndr. 2007;44:451–5.View ArticlePubMedGoogle Scholar
- Sun HY, Chen MY, Hsieh SM, Sheng WH, Chang SY, Hsiao CF, Hung CC, Chang SC. Changes in the clinical spectrum of opportunistic illnesses in persons with HIV infection in the era of highly active antiretroviral therapy. Jpn J Infect Dis. 2006;59:311–6.PubMedGoogle Scholar
- Manosuthi W, Chaovavanich A, Tansphaswadikul S, Prasithsirikul W, Inthong Y, Chottanapund S, Sittibusaya C, Moolasart V, Termvises P, Sungkanuparph S. Incidence and risk factors of major opportunistic infections after initiation of antiretroviral therapy among advanced HIV infected patients in a resource-limited setting. J Infect. 2007;55:464–9.View ArticlePubMedGoogle Scholar
- Paton NI, Sangeetha S, Earnest A, Bellamy R. The impact of malnutrition on survival and the CD4 count response in HIV-infected patients starting antiretroviral therapy. HIV Med. 2006;7:323–30.View ArticlePubMedGoogle Scholar
- Moges N, Kassa M. Prevalence of opportunistic infections and associated factors among HIV positive patients taking anti-retroviral therapy in DebreMarkos Referral Hospital, Northwest Ethiopia. J AIDS Clin Res. 2014;5:301.View ArticleGoogle Scholar
- Wongkoon S, Jaroensutasinee M, Jaroensutasinee K. Development of temporal modeling for prediction of dengue infection in Northeastern Thailand. Asian Pac J Trop Med. 2012;2012:249–52.View ArticleGoogle Scholar