Forecast analysis of any opportunistic infection among HIV positive individuals on antiretroviral therapy in Uganda
BMC Public Health volume 16, Article number: 766 (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.
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
Four TASO HIV clinics were purposively selected basing on volume, quality of data and geographical representation. The HIV clinics selected were TASO Mulago HIV clinic in central Uganda, TASO Mbarara HIV clinic in south-western Uganda, TASO Tororo HIV clinic in Eastern Uganda and TASO Gulu HIV clinic in Northern Uganda (Fig. 1). All HIV positive adults (15 years and above) who attended at least once at the selected HIV clinics in the period from 1st January 2004 to 31stDecember 2013 were included in the study.
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.
Between 2004 and 2013, a total of 36,133 HIV patients were enrolled on HAART of which two thirds (66 %) were female with a median age of 33 year (IQR, 27–40) (Table 1). In the planning data (2004–2010), it was observed that mean annual any OI prevalence reduced from 56.62 % in 2004 to 36.61 % in 2010. While in the validation data (2011–2013) mean annual any OI prevalence reduced from 35.89 % in 2011 to 27.53 % in 2013 (Table 2).
A time plot of the monthly OI prevalence trends shows several minor peaks along the series (Fig. 2). No seasonal or periodic components were clearly seen in the plot but the smoothed series generally depict a decreasing trend (Z statistic = −10.23, p <0.0001, nptrend) (Fig. 2). The Augmented Dickey-Fuller (ADF) test shows that the original monthly series had a unit root (z (t) = −2.353, p = 0.1556, lags = 20) suggesting that it’s not stationary. However after the first differencing the ADF test rejected the null hypothesis of a unit root in the series (z (t) = −16.7, p <0.001, lags = 20) implying that the differenced monthly series were stationary (Fig. 2). All further statistical procedures were performed on the stationary series.
Model identification started with autocorrelation analysis. Plots of autocorrelation function (ACF) and partial autocorrelation function (PACF) (Fig. 3) showed only the first lag of the ACF was significant (i.e. laying outside the grey 95 % CI band). It was also observed that the first few lags of PACF were decaying with time. Based on the autocorrelation structure, several potential models were identified.
Using the Akaike information criteria (AIC) and Bayesian information criteria (BIC), several first order differenced models were evaluated. ARIMA (1, 1, 1), ARIMA (1, 1, 2) and ARIMA (2,1,1) were initially selected basing on their low AIC and BIC values (Table 3). Parameter estimation for ARIMA (1, 1, 1), showed that AR (1) coefficient of 0.23 and MA (1) coefficient of −0.86 were significant (p <0.05) (Table 4).
Plots of autocorrelation function (ACF) and partial autocorrelation (PACF) of the model residuals showed that ACFs and PACFs were not significantly different from zero (white noise) (Fig. 4). The portmanteau Q-test (Q15 = 8.20, p = 0.9153) was also in favour of the null hypothesis of no autocorrelation in the residuals. A plot of the residual histogram showed no volatility clustering and so assumed the residuals were homoscedastic (Fig. 5). The skewness-kurtosis test (X2 = 3.42, p = 0.154) was consistent with the normality assumption of the model residuals (Fig. 4.13).
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).
A plot of the observed and predicted values of any OI matched reasonably well (Fig. 6). A t-test of the observed and predicted mean annual any OI prevalence showed no significant difference (p >0.05) (Table 5). The predicted linear trend from the model and the fitted linear trend to the observed data had the same slope = −0.0023 (95 % CI: −0.0028 to −0.0019) (Fig. 7).
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 (Table 6). The prediction error estimated in terms of mean absolute percentage error (MAPE) was 8.2 % (SD = 6.0).
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
WHO. HIV/AIDS facts sheet. Geneva: WHO; 2015.
UNAIDS. 90-90-90 : an ambitious treatment target to help end the AIDS epidemic. Geneva: WHO; 2014.
UNAIDS. Fast-Track: ending the AIDS epidemic by 2030. Geneva: UNAIDS; 2015.
UNAIDS. Uganda HIV/AIDS Country Progress Report 2014. Geneva: WHO; 2015.
Uganda AIDS Commission (UAC). National HIV and AIDS strategic plan 2015/2016- 2019/2020. Kampala: UAC; 2015.
UNAIDS. Access to Antiretroviral therapy in Africa: Status report on progress towards 2015 targets. Geneva: UNAIDS; 2014.
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.
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.
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.
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.
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.
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.
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.
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.
Stata corp. Time series analysis: introduction to autoregressive intergrated moving average (ARIMA) models. College station: Stata corp LP; 2012.
Chatfield C. Time series forecasting. Boca Raton: Chapman & Hall/CRC; 2001.
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.
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.
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.
TASO. TASO annual report. Kampala: TASO; 2013.
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.
Uganda Ministry of Health. Antiretroviral treatment policy for Uganda. Kampala: Uganda Ministry of Health; 2003.
WHO. Antiretroviral therapy for HIV infection in adults and adolescents: recommendations for a public health approach 2010 revision. Geneva: WHO; 2010.
Uganda Ministry of Health. Antiretroviral treatment policy. Kampala: Uganda Ministry of Health; 2009.
Box, Jenkins, Reinsel, editors. Time series analysis, forecasting and control. 3rd ed. Englewood Clifs: Prentice Hall; 1994.
Box G, Jenkins G. Time series analysis, forecasting and control. San Francisco: Holden-Day; 1976.
Chatfield C. The analysis of time series: an introduction. Washington: Chapman & Hall/CRC; 2004.
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.
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.
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.
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.
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.
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.
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.
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.
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.
About this article
Cite this article
Rubaihayo, J., Tumwesigye, N.M., Konde-Lule, J. et al. Forecast analysis of any opportunistic infection among HIV positive individuals on antiretroviral therapy in Uganda. BMC Public Health 16, 766 (2016). https://doi.org/10.1186/s12889-016-3455-5
- Opportunistic infections
- Antiretroviral therapy