Development of temporal modelling for forecasting and prediction of malaria infections using time-series and ARIMAX analyses: a case study in endemic districts of Bhutan.
BACKGROUND
Malaria still remains a public health problem in some districts of Bhutan despite marked reduction of cases in last few years. To strengthen the country's prevention and control measures, this study was carried out to develop forecasting and prediction models of malaria incidence in the endemic districts of Bhutan using time series and ARIMAX.


METHODS
This study was carried out retrospectively using the monthly reported malaria cases from the health centres to Vector-borne Disease Control Programme (VDCP) and the meteorological data from Meteorological Unit, Department of Energy, Ministry of Economic Affairs. Time series analysis was performed on monthly malaria cases, from 1994 to 2008, in seven malaria endemic districts. The time series models derived from a multiplicative seasonal autoregressive integrated moving average (ARIMA) was deployed to identify the best model using data from 1994 to 2006. The best-fit model was selected for each individual district and for the overall endemic area was developed and the monthly cases from January to December 2009 and 2010 were forecasted. In developing the prediction model, the monthly reported malaria cases and the meteorological factors from 1996 to 2008 of the seven districts were analysed. The method of ARIMAX modelling was employed to determine predictors of malaria of the subsequent month.


RESULTS
It was found that the ARIMA (p, d, q) (P, D, Q)s model (p and P representing the auto regressive and seasonal autoregressive; d and D representing the non-seasonal differences and seasonal differencing; and q and Q the moving average parameters and seasonal moving average parameters, respectively and s representing the length of the seasonal period) for the overall endemic districts was (2,1,1)(0,1,1)12; the modelling data from each district revealed two most common ARIMA models including (2,1,1)(0,1,1)12 and (1,1,1)(0,1,1)12. The forecasted monthly malaria cases from January to December 2009 and 2010 varied from 15 to 82 cases in 2009 and 67 to 149 cases in 2010, where population in 2009 was 285,375 and the expected population of 2010 to be 289,085. The ARIMAX model of monthly cases and climatic factors showed considerable variations among the different districts. In general, the mean maximum temperature lagged at one month was a strong positive predictor of an increased malaria cases for four districts. The monthly number of cases of the previous month was also a significant predictor in one district, whereas no variable could predict malaria cases for two districts.


CONCLUSIONS
The ARIMA models of time-series analysis were useful in forecasting the number of cases in the endemic areas of Bhutan. There was no consistency in the predictors of malaria cases when using ARIMAX model with selected lag times and climatic predictors. The ARIMA forecasting models could be employed for planning and managing malaria prevention and control programme in Bhutan.
Wangdi
Kinley
K
Department of Tropical Hygiene, Faculty of Tropical Medicine, Mahidol University, Bangkok, Thailand.
Singhasivanon
Pratap
P
Silawan
Tassanee
T
Lawpoolsri
Saranath
S
White
Nicholas J
NJ
Kaewkungwal
Jaranit
J
D001631:Bhutan / Q000453:epidemiology; D005544:Forecasting; D006801:Humans; D015994:Incidence; D008288:Malaria / Q000453:epidemiology* / Q000635:transmission*; D015233:Models, Statistical; D012189:Retrospective Studies; D012621:Seasons; D013997:Time Factors; D014887:Weather
J Vector Borne Dis. 2006 Mar;43(1):21-8
16642782
Lancet. 2005 Apr 23-29;365(9469):1487-98
15850634
Am J Trop Med Hyg. 2001 Jan-Feb;64(1-2 Suppl):1-11
11425172
Malar J. 2008 May 06;7:76
18460204
Science. 2001 Jul 6;293(5527):48-9
11444289
Am J Trop Med Hyg. 2007 Jan;76(1):33-8
17255225
Am J Trop Med Hyg. 2005 Mar;72(3):256-62
15772317
Public Health Rep. 2003 Jan-Feb;118(1):65-71
12604766
Lancet. 2003 May 17;361(9370):1705-6
12767739
Environ Health Perspect. 1995 May;103(5):458-64
7656875
Malar J. 2007 Sep 24;6:129
17892540
Trop Med Int Health. 1998 Oct;3(10):818-27
9809915
Acta Trop. 2009 Feb;109(2):118-23
19013420
Trans R Soc Trop Med Hyg. 2006 Jul;100(7):623-31
16406037
Malar J. 2008 Jun 19;7:110
18565224
BMC Public Health. 2007 Jun 08;7:101
17559638
Lancet. 2006 Mar 11;367(9513):859-69
16530580
Int J Health Geogr. 2006 May 15;5:20
16700905
J Infect. 2006 Feb;52(2):131-9
16442438
Cochrane Database Syst Rev. 2004;(2):CD000363
15106149
Nature. 2005 Mar 10;434(7030):214-7
15759000
Trop Med Int Health. 2008 Dec;13(12):1501-8
19000157
East Mediterr Health J. 2007 Jan-Feb;13(1):17-24
17546901
Malar J. 2008 May 06;7:77
18460205
Int J Health Geogr. 2009 Jan 26;8:5
19171051