Daily forecast of dengue fever incidents for urban villages in a city
© Chan et al.; licensee BioMed Central. 2015
Received: 15 December 2014
Accepted: 26 January 2015
Published: 31 January 2015
Instead of traditional statistical models for large spatial areas and weekly or monthly temporal units, what public health workers urgently need is a timely risk prediction method for small areas. This risk prediction would provide information for early warning, target surveillance and intervention.
Daily dengue cases in the 457 urban villages of Kaohsiung City, Taiwan from 2009 to 2012 were used for model development and evaluation. There were in total 2,997 confirmed dengue cases during this period. A logistic regression model was fitted to the daily incidents occurring in the villages for the past 30 days. The fitted model was then used to predict the incidence probabilities of dengue outbreak for the villages the next day. A percentile of the 457*30 fitted incidence probabilities was chosen to determine a cut-point for issuing the alerts. The covariates included three different levels of spatial effect, and with four lag time periods. The population density and the meteorological conditions were also included for the prediction.
The performance of the prediction models was evaluated on 122 consecutive days from September 1 to December 31, 2012. With the 80th percentile threshold, the median sensitivity was 83% and the median false positive rate was 23%. We found that most of the coefficients of the predictors of having cases at the same village in the previous 14 days were positive and significant for the 122 daily updated models. The estimated coefficients of population density were significant during the peak of the epidemic in 2012.
The proposed method can provide near real-time dengue risk prediction for a small area. This can serve as a useful decision making tool for front-line public health workers to control dengue epidemics. The precision of the spatial and temporal units can be easily adjusted to different settings for different cities.
KeywordsLogistic regression Dengue Dynamic threshold Decision support
The increasing economic and disease burden of dengue not only is an important public health issue in tropical and sub-tropical countries , but it has also become an increasing threat to intemperate countries in Europe and North America due to the effect of rising temperatures [2, 3]. The World Health Organization (WHO) estimated that there are 50–100 million dengue infections globally every year, including 500,000 people with severe conditions that require hospitalization, and about 2.5 % of those severe cases die . One study in Southeast Asia found that the annual economic burden was US$950 million (95% certainty level: US$610 million- US$1,384 million) and disability-adjusted life years (DALYs) per million inhabitants was 372 (210–520) . Dengue is a vector-borne viral disease transmitted by vectors such as Aedes aegypti and Aedes albopictus to humans. So far, environmental cleaning or peridomestic space spraying of insecticides have been the two major approaches for controlling and preventing dengue epidemics in communities . Although insecticide resistance is a critical problem for controlling epidemics or outbreaks , the combined strategy of proper chemical spraying  and removal of vector-breeding sites  can indeed minimize risks of dengue infection. In addition to such control strategies, public health workers urgently need to know where the next high-risk areas are in order for them to monitor intensively, identify the cases early, and clean out the environment as soon as possible.
Therefore, the dengue surveillance involves routine collection on both the clinical infection cases and mosquito density, such as the House Index (HI), Container Index (CI) and Breteau Index (BI) . However, because the mosquito surveillance is not systematic and does not collect data in real time, this indicator is not very sensitive or reliable for risk prediction . In the empirical findings, environmental factors such as temperature and precipitation are correlated with the dengue virus’ activities and vector’s life cycle [12–14]. The other challenge for dengue prevention is the unpredictability of human movements. Unlike people, mosquitoes tend to have short movement distances (i.e. flying mostly shorter than 150 meters) . Moreover, the infected persons usually find it difficult to recall where they got their mosquito bites. Therefore, we can only use the location of confirmed dengue cases coupled with the incubation time lag to differentiate the possible previous local or neighboring effect on the future risk of a dengue outbreak. In addition, the underlying transportation network between the sites of infection (i.e. where the person was bitten) and their residence might be estimated from the cumulative spatio-temporal trajectory of the confirmed cases. House-to-house human interaction is an important factor in dengue virus transmission . Population density may also play a role in increasing the risk of dengue infection . It would be beneficial to consider all these factors to enhance the accuracy of risk prediction.
Instead of traditional statistical models for large spatial areas and weekly or monthly temporal units, what public health workers need is a timely and local risk prediction method. The alert threshold should evolve with time, and the risk level should be adjusted throughout the different stages of an epidemic. This method should balance accuracy and false positives, and can help with early warning, the allocation of resources, risk communication and disease control. This study proposes a simple and timely method to predict incidence probability of dengue outbreak at a local level.
This study was approved by the institutional review board (IRB) of Academia Sinica (IRB#: AS-IRB01-13060). The databases we used were all stripped of identifying information and thus informed consent was not needed.
We denoted the daily reported dengue cases at the ith urban village by y it for the tth day. Using the daily reported cases at the n villages in the city for m consecutive days before the current day, say d - 1, we constructed a logistic regression model for predicting the probabilities Pr(Y id > 0) and determining a threshold to classify the villages as having cases or not on the dth day. In addition to limited covariates such as population density and weather variables, we tried to construct various predictors using the reported cases in the affected villages. Basically, we considered three categories of autoregressive factors of cases reported in the four periods, which are the previous day, 2–7 days, 8–14 days and 15–30 days before the current day among the n villages. We decided these time intervals based on three criteria. First, we considered the extrinsic and intrinsic incubation period of dengue infection. The mean intrinsic incubation period (IIP) has been determined to be 5.9 days (95% CI: 3–10 days), with an extrinsic incubation period (EIP) at 25°C of 5 days to 33 days . In general, 75% of infected patients have been found to develop symptoms by 7.1 days (95% CI: 6.7–7.6) . That’s why we used 7 days for the first cut-points, and considered the EIP for the other 14 days. Second, we considered both the maximum lifetime of the Aedes aegypti, 22 days and IIP . Thus, we decided on 30 days as our maximum observation window. Third, from a model selection perspective, we have also tried modeling with different combinations of time intervals, and found the current interval had good performance and is easy to remember in practical usage.
On day d, the logistic regression model was fitted to the dengue cases y it for i = 1, 2, …, n and t = d - 1, …, d - m. We have considered four weather covariates in the regression model. If all the m values of a weather covariate are zero, the covariate is excluded from the model. We applied the glm function in the R package to estimate the model coefficients and calculated the m × n estimates of the probabilities Pr(Y it > 0). We denoted c p as the pth percentile of those estimates of Pr(Y it > 0) with y it = 0. The value c p will be used as a threshold for classifying predicted probabilities of the villages on the dth day as having cases or not. The predicted probabilities can be obtained simply from the logistic regression model with the estimated coefficients and new predictors A qid , B qid , C qid , Z jid and X jd . If the predicted probability at a village is larger than c p , we conclude that cases will occur at this village and none otherwise.
Summary of sensitivities and false positive rates under different threshold cut-points from the 122 predictive models during September 1 to December 31, 2012
(B) False positive rates
In this study, we propose a small area-based dengue prediction method with timeliness improved by the daily prediction, like a weather forecast. In addition to the labor-intensive  and costly  vector surveillance, this sensitive, comprehensive and near real-time disease surveillance is also crucial to enhancing the accuracy of the prediction and improving risk management of dengue epidemics and outbreaks. One previous study found that implementing a rapid measure for controlling dengue outbreaks in the early stage can significantly reduce the epidemic size . Therefore, this early warning statistical method will give front-line public health workers an advantage in identifying high-risk areas for intensive surveillance and early intervention.
Good quality of dengue surveillance is the basic infrastructure for any prediction model. In Taiwan, dengue is a category 2 notifiable infectious disease. The physicians have the responsibility to report the suspected dengue cases to the local health department within 24 hours of clinical diagnosis . Thus, the quality of surveillance was ensured by the Communicable Disease Prevention Act, and the efficiency was enhanced by the internet reporting system in recent years. According to a previous study, the case fatality rate of dengue hemorrhagic fever (DHF) ranged from 9.1% to 50% during 2002–2007 . The ratio of symptomatic and asymptomatic cases is 1.78 for adult dengue in Taiwan, which means about two-thirds of the dengue-infected adults have clinical symptoms . Even when patients had clinical symptoms, a few cases would still be missed by physicians in hospitals . Continuous education for physicians will be needed to reduce the under-reporting of dengue cases.
In the real situation, the dengue infected cases will be either asymptomatic or might be under-reported by the clinicians. Thus, when the surveillance system detects one confirmed case, it may mean some infected cases also exist in the same village or other villages. The current policy will initiate an epidemiological investigation after receiving a report of suspected or confirmed dengue cases, and the investigation is required to be completed within 24 hours . At the same time, the environmental cleaning and mosquitos’ surveillance also need to be completed within 48 hours. The effect of peridomestic space spraying does not persist very long. The major tasks suggested by Taiwan’s CDC are source reduction [26, 27]. In relation to our prediction model, we can help the decision maker in the local health department to decide where the first priority villages to do intensive source reduction and environmental cleaning are. With these approaches, instead of chemical spraying, economic loss can be minimized and the participation and acceptance rate in the village will be elevated.
Unlike other dengue surveillance models with larger spatial and temporal units [28, 29], a small area or village (mean size: 0.36 km2) was our major targeted and surveillance unit. Therefore, we quantify the influence of local villages using different time lags. We found that the local influence was very strong from the previous day to 14 days. This was quite consistent with the incubation period of dengue infection . We also found that when the epidemic became more severe, the estimated coefficient of the mean number of cases at the local village also increased positively and significantly. This finding meant that the local clustering and transmission was really present [30, 31]. The neighboring and other influential villages were not at all significant. However, including these two factors into the model will enhance the sensitivity and incorporate the risk estimation when the other villages have dengue infection. It might somewhat capture the underlying human movements. On the other hand, the specificity was also very important in minimizing false positives. In this study, we considered the fact that some villages had no confirmed cases of dengue in the previous 30 days, helping the model to identify and adjust for those low risk areas.
The tradeoff between sensitivity and false positive rate is always a difficult issue for the model prediction. In this study, the key point for determining these two values was by the thresholds or cut-points. The presented median sensitivity and false positive rate were obtained based on the threshold corresponding to an expected 20% false alarm rate. The sensitivity we presented here was 83% and specificity was 77% (1-false positive rate). In other dengue prediction models [28, 32], the sensitivities were around 60% and the specificities were around 97%. If we choose the threshold corresponding to an expected 10% false alarm rate, the sensitivity will decline to 66.7% but our specificity can be elevated to near 90%. Most of the other models were applied to weekly data and large regions or cities. The challenge for the small area and daily data in this study might be even higher than the current methods.
Although climate conditions such as temperature, precipitation, humidity, and sea surface temperature were found to be highly correlated with the cases of dengue infection [33–35], the influence of meteorological conditions was not detected in this study. Our study area was confined to one city with only one weather station. Kaohsiung City is located in the tropical climate zone, so weather conditions do not change very much with time, and there were no weather variations in the small areas we studied. This might be the reason why the weather conditions could not significantly explain dengue epidemics in small areas.
A complete dengue transmission network requires vectors (mosquitoes) and human beings. However, data from the vector surveillance were not systematically collected, and might be misleading if the surveillance was implemented after cases of dengue were reported or after insecticide spraying. Previous studies also found that the vector surveillance data were not good enough for identifying any correlation with dengue epidemics . Therefore, we did not include this information into our final model. The other important factor was population density. One previous study found that around 3,000-7,000 persons/km2 in Vietnam was the condition for the highest risk for dengue epidemics, and the risk did not increase with density beyond that . Some studies found that dengue cases clustered in urban areas . In our studied area, the population densities were higher than the mentioned one in Vietnam. During our model validation period, the population density only showed a significant influence during the peak of the epidemic in 2012. High population density might only facilitate the dengue transmission cycle between human beings and mosquitoes when the number of infected cases is large enough, and the trend will still need to be validated by more data in the future.
Although this proposed method can be easily applied to routine dengue surveillance, there are still many challenges limiting the prediction power of the model. The first one is not knowing where the actual infection of each dengue case took place. The mosquitoes themselves cannot fly long distances, but the infected human may travel longer distances within the city or even beyond, and subsequently spread the infection. In this study, we have tried to use three different levels of influence, namely the local, neighboring, and other influential villages to capture possible impact of human movement. In this approach, we borrowed information mainly from the dengue surveillance data. If the information from epidemiological investigation can be applied to the surveillance model, the dynamic pattern of human movement might increase its prediction power. The second limitation is the lack of completed vector surveillance data. Although the current data were not suitable for the prediction, systematic collection of vector surveillance will be needed and might be beneficial for risk prediction in the future.
Timely risk prediction and early intervention in the local villages are important and urgently needed to control dengue epidemics. With simple logistic regression and a dynamic alert threshold, sensitivity of nearly 80% and a false positive rate of 20% can be achieved in daily dengue risk prediction. The local effect of the same village was the significant predictor of dengue outbreaks within 14 days after the latest dengue infection occurred. In addition, this small area-based risk prediction will be practical for targeted surveillance and intervention. It will serve as a useful and near real-time decision-making tool for front-line public health workers. The precision of the spatial and temporal units can be easily adjusted to different settings for different cities.
This research is supported by a grant from the Ministry of Science and Technology, Taiwan (MOST 103-2621-M-001 -002). We gratefully acknowledge the Data Bank for Atmospheric Research, sponsored by the National Science Council, Taiwan, for providing the daily weather data. In addition, we thank the Department of Health, Kaohsiung City Government for the administrative support.
- Bhatt S, Gething PW, Brady OJ, Messina JP, Farlow AW, Moyes CL, et al.: The global distribution and burden of dengue. Nature 2013,496(7446):504–7. 10.1038/nature12060View ArticlePubMed CentralPubMedGoogle Scholar
- Lourenco J, Recker M: The 2012 Madeira dengue outbreak: epidemiological determinants and future epidemic potential. PLoS Negl Trop Dis 2014,8(8):e3083. 10.1371/journal.pntd.0003083View ArticlePubMed CentralPubMedGoogle Scholar
- Liu-Helmersson J, Stenlund H, Wilder-Smith A, Rocklov J: Vectorial capacity of Aedes aegypti: effects of temperature and implications for global dengue epidemic potential. PLoS One 2014,9(3):e89783. 10.1371/journal.pone.0089783View ArticlePubMed CentralPubMedGoogle Scholar
- Dengue and severe dengue: Dengue and severe dengue. [http://www.who.int/mediacentre/factsheets/fs117/en/]
- Shepard DS, Undurraga EA, Halasa YA: Economic and disease burden of dengue in Southeast Asia. PLoS Negl Trop Dis 2013,7(2):e2055. 10.1371/journal.pntd.0002055View ArticlePubMed CentralPubMedGoogle Scholar
- Esu E, Lenhart A, Smith L, Horstick O: Effectiveness of peridomestic space spraying with insecticide on dengue transmission; systematic review. Trop Med Int Health 2010,15(5):619–31.PubMedGoogle Scholar
- Maciel-de-Freitas R, Avendanho FC, Santos R, Sylvestre G, Araujo SC, Lima JB, et al.: Undesirable consequences of insecticide resistance following Aedes aegypti control activities due to a dengue outbreak. PLoS One 2014,9(3):e92424. 10.1371/journal.pone.0092424View ArticlePubMed CentralPubMedGoogle Scholar
- Oki M, Sunahara T, Hashizume M, Yamamoto T: Optimal timing of insecticide fogging to minimize dengue cases: modeling dengue transmission among various seasonalities and transmission intensities. PLoS Negl Trop Dis 2011,5(10):e1367. 10.1371/journal.pntd.0001367View ArticlePubMed CentralPubMedGoogle Scholar
- Ocampo CB, Mina NJ, Carabali M, Alexander N, Osorio L: Reduction in dengue cases observed during mass control of Aedes (Stegomyia) in street catch basins in an endemic urban area in Colombia. Acta Trop 2014, 132:15–22.View ArticlePubMedGoogle Scholar
- Sarfraz MS, Tripathi NK, Tipdecho T, Thongbu T, Kerdthong P, Souris M: Analyzing the spatio-temporal relationship between dengue vector larval density and land-use using factor analysis and spatial ring mapping. BMC Public Health 2012., 12: Google Scholar
- Bowman LR, Runge-Ranzinger S, McCall PJ: Assessing the relationship between vector indices and dengue transmission: a systematic review of the evidence. PLoS Negl Trop Dis 2014,8(5):e2848. 10.1371/journal.pntd.0002848View ArticlePubMed CentralPubMedGoogle Scholar
- Carrington LB, Armijos MV, Lambrechts L, Barker CM, Scott TW: Effects of Fluctuating Daily Temperatures at Critical Thermal Extremes on Aedes aegypti Life-History Traits. PLoS One 2013.,8(3): Google Scholar
- Costa EAPD, Santos EMD, Correia JC, de Albuquerque CMR: Impact of small variations in temperature and humidity on the reproductive activity and survival of Aedes aegypti (Diptera, Culicidae). Rev Bras Entomol 2010,54(3):488–93. 10.1590/S0085-56262010000300021View ArticleGoogle Scholar
- Wee LK, Weng SN, Raduan N, Wah SK, Ming WH, Shi CH, et al.: Relationship between rainfall and Aedes larval population at two insular sites in Pulau Ketam, Selangor, Malaysia. SE Asian J Trop Med 2013,44(2):157–66.Google Scholar
- Harrington LC, Scott TW, Lerdthusnee K, Coleman RC, Costero A, Clark GG, et al.: Dispersal of the dengue vector Aedes aegypti within and between rural communities. Am J Trop Med Hyg 2005,72(2):209–20.PubMedGoogle Scholar
- Stoddard ST, Forshey BM, Morrison AC, Paz-Soldan VA, Vazquez-Prokopec GM, Astete H, et al.: House-to-house human movement drives dengue virus transmission. Proc Natl Acad Sci U S A 2013,110(3):994–9. 10.1073/pnas.1213349110View ArticlePubMed CentralPubMedGoogle Scholar
- Schmidt WP, Suzuki M, Thiem VD, White RG, Tsuzuki A, Yoshida LM, et al.: Population density, water supply, and the risk of dengue fever in Vietnam: cohort study and spatial analysis. PLoS Med 2011,8(8):e1001082. 10.1371/journal.pmed.1001082View ArticlePubMed CentralPubMedGoogle Scholar
- Chan M, Johansson MA: The incubation periods of Dengue viruses. PLoS One 2012,7(11):e50972. 10.1371/journal.pone.0050972View ArticlePubMed CentralPubMedGoogle Scholar
- Rudolph KE, Lessler J, Moloney RM, Kmush B, Cummings DA: Incubation periods of mosquito-borne viral infections: a systematic review. Am J Trop Med Hyg 2014,90(5):882–91. 10.4269/ajtmh.13-0403View ArticlePubMedGoogle Scholar
- Brady OJ, Golding N, Pigott DM, Kraemer MU, Messina JP, Reiner RC Jr, et al.: Global temperature constraints on Aedes aegypti and Ae. albopictus persistence and competence for dengue virus transmission. Parasites and Vectors 2014, 7:338. 10.1186/1756-3305-7-338View ArticlePubMed CentralPubMedGoogle Scholar
- Sivagnaname N, Gunasekaran K: Need for an efficient adult trap for the surveillance of dengue vectors. Indian J Med Res 2012,136(5):739–49.PubMed CentralPubMedGoogle Scholar
- Pepin KM, Marques-Toledo C, Scherer L, Morais MM, Ellis B, Eiras AE: Cost-effectiveness of novel system of mosquito surveillance and control, Brazil. Emerg Infect Dis 2013,19(4):542–50. 10.3201/eid1904.120117View ArticlePubMed CentralPubMedGoogle Scholar
- Chowell G, Fuentes R, Olea A, Aguilera X, Nesse H, Hyman JM: The basic reproduction number R0 and effectiveness of reactive interventions during dengue epidemics: the 2002 dengue outbreak in Easter Island, Chile. Math Biosci Eng 2013,10(5–6):1455–74.PubMedGoogle Scholar
- Lin CC, Huang YH, Shu PY, Wu HS, Lin YS, Yeh TM, et al.: Characteristic of dengue disease in Taiwan: 2002–2007. Am J Trop Med Hyg 2010,82(4):731–9. 10.4269/ajtmh.2010.09-0549View ArticlePubMed CentralPubMedGoogle Scholar
- Lee MS, Hwang KP, Chen TC, Lu PL, Chen TP: Clinical characteristics of dengue and dengue hemorrhagic fever in a medical center of southern Taiwan during the 2002 epidemic. J Microbiol Immunol Infect 2006,39(2):121–9.PubMedGoogle Scholar
- Taiwan CDC: Guidelines for Dengue Control. Edited by: Taiwan CDC. Taipei, Taiwan: Department of Health; 2012.Google Scholar
- Tuan Yen-Chang WY-Z, Chi-Chuan H, Chin-Hsien W, Chao-Ching C, Li-Jen L, Tsuey-Feng L, et al.: An analysis of the Implementation of the Source Reduction Program in Response to 2010 Dengue Epidemic. Taiwan Epidemiol Bull 2011 2011,27(10):164–76.Google Scholar
- Buczak AL, Baugher B, Babin SM, Ramac-Thomas LC, Guven E, Elbert Y, et al.: Prediction of high incidence of dengue in the Philippines. PLoS Negl Trop Dis 2014,8(4):e2771. 10.1371/journal.pntd.0002771View ArticlePubMed CentralPubMedGoogle Scholar
- Restrepo AC, Baker P, Clements AC: National spatial and temporal patterns of notified dengue cases, Colombia 2007–2010. Trop Med Int Health 2014,19(7):863–71. 10.1111/tmi.12325View ArticlePubMedGoogle Scholar
- Wen TH, Lin NH, Chao DY, Hwang KP, Kan CC, Lin KC, et al.: Spatial-temporal patterns of dengue in areas at risk of dengue hemorrhagic fever in Kaohsiung, Taiwan, 2002. Int J Infect Dis 2010,14(4):e334–343. 10.1016/j.ijid.2009.06.006View ArticlePubMedGoogle Scholar
- Vazquez-Prokopec GM, Kitron U, Montgomery B, Horne P, Ritchie SA: Quantifying the spatial dimension of dengue virus epidemic spread within a tropical urban environment. PLoS Negl Trop Dis 2010,4(12):e920. 10.1371/journal.pntd.0000920View ArticlePubMed CentralPubMedGoogle Scholar
- Buczak AL, Koshute PT, Babin SM, Feighner BH, Lewis SH: A data-driven epidemiological prediction method for dengue outbreaks using local and remote sensing data. BMC Med Inform Decis 2012., 12: Google Scholar
- Xu HY, Fu X, Lee LK, Ma S, Goh KT, Wong J, et al.: Statistical modeling reveals the effect of absolute humidity on dengue in Singapore. PLoS Negl Trop Dis 2014,8(5):e2805. 10.1371/journal.pntd.0002805View ArticlePubMed CentralPubMedGoogle Scholar
- Eastin MD, Delmelle E, Casas I, Wexler J, Self C: Intra- and interseasonal autoregressive prediction of dengue outbreaks using local weather and regional climate for a tropical environment in Colombia. Am J Trop Med Hyg 2014,91(3):598–610. 10.4269/ajtmh.13-0303View ArticlePubMed CentralPubMedGoogle Scholar
- Shang CS, Fang CT, Liu CM, Wen TH, Tsai KH, King CC: The role of imported cases and favorable meteorological conditions in the onset of dengue epidemics. PLoS Negl Trop Dis 2010,4(8):e775. 10.1371/journal.pntd.0000775View ArticlePubMed CentralPubMedGoogle Scholar
- Jeefoo P, Tripathi NK, Souris M: Spatio-temporal diffusion pattern and hotspot detection of dengue in Chachoengsao province, Thailand. Int J Environ Res Public Health 2011,8(1):51–74.View ArticlePubMed CentralPubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. 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.