Determining the effects of preseasonal climate factors toward dengue early warning system in Bangladesh ======================================================================================================= * M. Pear Hossain * Wen Zhou * Chao Ren * John Marshall * Hsiang-Yu Yuan ## Abstract The incidence of dengue has increased rapidly in Bangladesh since 2010 with an outbreak in 2018 reaching a historically high number of cases, 10,148. Bangladesh is located to the northeast of India and spans both tropical and subtropical regions. A better understanding of the effects of preseasonal climate variability on the increasing incidence of dengue in Bangladesh may provide insight into whether dengue has been expanding from tropical to temperate regions through the subtropics in Southeast Asia and can enable early warning for preventative measures. We developed a generalized linear model to predict the number of annual dengue cases based on monthly temperature, rainfall and sunshine prior to dengue season. Variable selection and leave-one-out cross-validation were performed to identify the best prediction model and to evaluate the model’s performance. Our model successfully predicted the largest outbreak in 2018, with 10,077 cases (95% CI: [9,912–10,276]), in addition to smaller outbreaks in five different years (2003, 2006, 2010, 2012 and 2014) and successfully identified the increasing trend in cases between 2010 and 2018. We found that the temperature has an apparent dual effect on the annual incidence of dengue, with a positive effect seen prior to April and negative effect seen between April and June. Our results suggest that the optimal minimum temperature for mosquito population expansion is 21–23°*C*. This study has implications for understanding how climate variability has affected recent dengue expansion in northern India and Southeast Asia. **Author summary** Bangladesh has experienced an unprecedented outbreak of dengue in recent years. Climate factors are believed to have played a significant role in the increased mosquito population, the primary vector of dengue. We have observed an upward trend in the number of dengue cases since 2010, which has continued until recent years. However, in previous studies of dengue–climate association, this link was not investigated. Therefore, we developed a model that uses this information to assess the association of annual dengue incidence with climate factors such as temperature, rainfall and sunshine duration. Our findings suggest that warmer springs, with minimum monthly temperatures of 21–23°*C*, are more favorable for mosquito population expansion and subsequent dengue transmission. Shorter duration of sunshine increase the risk of dengue, whereas abundant rainfall reduces the incidence of dengue. This research will help increase understanding of the effect of climate variability on dengue expansion, not only in Bangladesh but also in northern India and Southeast Asia. Keywords * climate change * disease forecasting * dengue outbreak * dengue prediction * warm winter ## Background Dengue fever, one of the most prevalent vector-borne diseases, has led to significant socio-economic costs in many parts of the world [1]. Three-quarters of the global dengue cases occur in Southeast Asian and western Pacific countries, due to the associated favorable weather conditions for mosquito population expansion [2, 3]. Outbreaks of dengue fever can significantly reduce life expectancy due to the possibility of developing severe dengue following secondary infections from different dengue serotypes [4]. Therefore, it is critically important to understand the impacts of the climate on the spread of dengue in these regions, as this can serve as early warning system and enable early preventative measures to be put in place before outbreaks become established. Expansion of dengue in the regions surrounding northern India may have occurred in recent years due to climate change. The prolonged rainy seasons and increasing temperatures in subtropical regions of Southeast Asia may provide favorable conditions for expansion of *Aedes* mosquito populations, the dengue vector [5–8]. In addition, increased incidence of dengue has recently been observed in more temperate regions [9], such as Nepal [10], indicating a possible expansion of the disease from the subtropics to cooler climates, posing a threat to northern India, Pakistan and their neighbors. Bangladesh is located to the northeast of India and to the south of Nepal, and lies along the Tropic of Cancer. Understanding the patterns of recent dengue outbreaks in Bangladesh may provide greater insight into whether dengue has expanded into the region surrounding northern India, a region with more than 140 million inhabitants. Dengue fever was first identified in Bangladesh in 1964 [11] and was not initially considered to be a severe threat to public health. However, in 2000, an outbreak occurred, leading to a total of 5,551 reported cases and 93 confirmed deaths [12, 13]. The average annual number of dengue cases decreased between 2000 and 2010. However, since then the number of annual dengue cases in Bangladesh has been increasing rapidly. A recent outbreak in 2019 was the largest ever experienced by the country, whereas the second largest outbreak was seen only a year prior, in 2018. Whether or how climate variability may have driven this unprecedented rise in outbreak size in 2018 and 2019 is still largely unknown [14]. Several studies have been carried out to estimate dengue incidence in Bangladeshusing climate data prior to 2010 [15–17]; however, the driving factors responsible for the increasing disease burden since 2010 remain to be investigated. In these previous studies, temperature and rainfall were found to be significant contributing factors [18–21]. Previous studies also assumed that the effects of climate variables are independent of the time of year. However, the effects of climate variables can also be time-dependent. Several studies have demonstrated that the effects of rainfall on dengue incidence can vary throughout the year [23, 24]. The abundant rainfall that occurs during monsoon season is likely to have negative effects on mosquito population size, as the rain can disrupt potential mosquito habitats. In contrast, rainfall in winter months may result in stagnant bodies of water suitable for mosquito breeding. Latest studies have mentioned the dengue incidence and mosquito abundance can be affected by weather conditions up to 5 months before the season starts [23, 24, 37] using data in or near subtropical areas. However, most of the studies focus on climate factors during dengue season. A better understanding of the effects of preseasonal climate variability on the increasing incidence of dengue can provide insight into whether dengue has been expanding in a region and allow early warning system to be built. This study aimed to estimate the effects of preseasonal time-dependent climate factors on annual dengue incidence in Bangladesh using historical data from 2000 to 2018. We developed a generalized linear model with a Poisson link function to predict annual dengue cases based on monthly temperature, rainfall and sunshine. We demonstrate that temperature and rainfall have variable effects on dengue incidence depending on the time of year and suggest an ideal temperature range for mosquito population growth based on our findings. ## Methods ### Study location Bangladesh is a Southeast Asian country, as defined by the World Health Organization (WHO) [25]. India surrounds it to the east, west and north, and Myanmar borders it to the southeast (Fig 1). The Bay of Bengal is located in the south. Bangladesh is located at 20°59′ *N* to 26°63′ *N* and 88°03′ *E* to 92°67′ *E*. The Tropic of Cancer line is located at 23°26′ *N* and 88°47′ *E*, where it crosses Bangladesh from east to west [26, 27]. Bangladesh is in a tropical climate region. The seasons in Bangladesh can be broadly characterized as summer (March–June; mostly hot and humid), monsoon (June–October; warm and rainy) and winter (October–March; cold and dry). However, March can also be described as the spring, and the duration between mid-October and mid-November can be called the autumn. The maximum temperature ranges from 30°*C* to 40°*C* during summer, whereas in winter the average temperature reaches as low as 10°*C* in most areas of the country. The average annual rainfall ranges between 1,500 mm and 3,000 mm. Approximately 70–80% of the annual rainfall occurs during monsoon season [28]. The shortest period of sunshine, 5.4–5.8 hours per day, also occurs during this season. In contrast, winter and summer have the longest sunshine duration, 8.5–9.1 hours per day [29]. ![Fig 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/22/2020.09.21.20199190/F1.medium.gif) [Fig 1.](http://medrxiv.org/content/early/2020/09/22/2020.09.21.20199190/F1) Fig 1. Location of Bangladesh on the world map. India surrounds Bangladesh on three sides (east, west and north) and the Bay of Bengal is to the south. The country shares a border with Myanmar in the southeast. The Tropic of Cancer line crosses the middle of the country. ### Dengue data Dengue cases observed in health facilities across the country are generally reported to the Directorate General of Health Services (DGHS), and are classified into suspected, probable and confirmed cases. Individuals having acute febrile illness with or without non-specific signs and symptoms are classified as suspected cases and those having acute febrile illness with serological diagnosis are considered probable cases. The confirmed cases should have an acute febrile illness with positive dengue NS1 antigen or PCR test. Details of dengue case definitions and management are available from the DGHS [30]. The communicable disease control (CDC) unit of the DGHS compiles the reported dengue cases on a daily basis for further circulation. We accessed monthly dengue cases between January 2000 and December 2018 from the DGHS by collaborating with the Institute of Epidemiology, Disease Control and Research. ### Climate data Accumulated weather information was monitored and managed by the Bangladesh Meteorological Department (BMD) at 35 distinct weather stations across the country. We collected these weather records from the BMD including daily mean, minimum and maximum temperature, total and maximum daily rainfall and daily sunshine duration. Temperature and rainfall are measured in degree Celsius (°*C*) and millimeter (mm), respectively, whereas sunshine duration is recorded in hours. Daily information was averaged for each month to obtain monthly information for each station. The national averages for monthly temperature, sunshine duration and rainfall were obtained by averaging the values of all 35 weather stations. ### Model formulation To predict annual dengue incidence, cases reported to the DGHS were considered the outcome variable. A generalized linear model with a Poisson link function was utilized to predict dengue incidence. We assumed that dengue cases reported in January, February and March were belonging to the previous year’s dengue outbreak. Hence, annual dengue incidence is defined as the sum of the number of dengue cases from April to December of a given year and from January to March of the following year. Let *y**j* be the annual dengue cases in the *j**th* year (*j* = 1, 2, …, 19) such that *y**j* ∼ *Poisson*(*λ**j*), where *λ**j* represents the expected number of dengue cases in the *j**th* year, i.e. *E*(*y**j*) = *λ**j*. Therefore, the Poisson regression model can be expressed as ![Formula][1] where *β**i*, *η**i* and *γ**i* represent coefficients of temperature (*T*), rainfall (*R*) and sunshine duration (*S*) in the *i**th* month (*i* = 1, 2, …, 6). Therefore, 15 predictor variables exist in the model as shown in Eq.(1). Monthly minimum temperature and total rainfall from January to June and the sunshine duration from April to June were chosen as potential predictor variables. These predictors were preselected during these months as we aimed to predict dengue outbreaks, which often begin during the early summer. ### Model selection Different models were compared for predicting the annual dengue outbreaks with various combinations of climate variables. A corrected version of the Akaike information criteria (*AIC**c*) [31] was used to extract potential predictor variables. The best prediction model was determined using a two-stage selection approach (Fig 2). In the first stage, variable selection was performed using stepwise *AIC**c* selection for each of the models. In the final step of the stepwise selection, each model provides a set of variables with minimum *AIC**c*. For these models, the model with the lowest *AIC**c* was selected as the lowest *AIC**c* model (*M**L*). In contrast, the models satisfying the condition ![Fig 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/22/2020.09.21.20199190/F2.medium.gif) [Fig 2.](http://medrxiv.org/content/early/2020/09/22/2020.09.21.20199190/F2) Fig 2. Two-stage selection of the best prediction model. *Min*(*AIC**c*) and *Min*(*MSE**V a*) refer to the lowest value of *AIC**c* and lowest value of mean squared error of *LOOCV* validation. ![Graphic][2] represents the average value of *AIC**c* for all candidate models. ![Formula][3] were defined as low *AIC**c* models (*M**l*), where *m*— is the average *AIC**c* of all candidate models. The reason we included both lowest and low *AIC**c* models is that we aim to identify the best prediction model using *LOOCV* among all the models with low *AIC**c*. In the second stage, the *M**L* and *M**l* models were compared using the mean squared errors obtained from leave-one-out cross-validation (*LOOCV*). The best model was the model with the minimum *MSE**Va* in *LOOCV*. ### Model validation Model validation was conducted using *LOOCV*. To perform *LOOCV*, the data for a specific test year were removed and the model was fitted based on the remaining data, which served as a training set. The fitted model was then used to predict the annual dengue cases for the test year. We repeated this procedure for all years from 2000 to 2018. Mean square errors for the validation set *MSE**V a* and for the training set *MSE**T r* were obtained by calculating the difference between predicted and observed numbers of annual dengue cases in the testing set and training set. Next, we checked the mean squared error ratio ![Graphic][4] the ratio of the mean squared errors of the validation set and the training set, to see whether it was less than the factor *c*, ![Graphic][5]. We computed bootstrap confidence intervals for the predicted annual dengue cases in each year. To do this, we simulated 1,000 random samples from a Poisson distribution by considering *LOOCV* -estimated annual cases ![Graphic][6] as the parameter of the distribution. The random numbers were used to refit the model 1,000 times, giving the distribution of the estimated parameters. The lower and upper bounds of the 95% confidence intervals were calculated based on the 2.5% and 97.5% quantiles of the parameter distributions. ### Assessment of the effects of climate factors Interpretation of the estimated model coefficients for a generalized linear model is not as straightforward as it is for an ordinary linear regression model, as the dependent variable *y* is associated with a link function, such as a Poisson link [32]. Therefore, we calculated the marginal effect at the mean (MEM) to understand the effect of each of the predictor variables separately using the *R*-package *ggeffect* [33]. ## Results ### Dengue cases in Bangladesh Dengue cases exhibited a decreasing trend since the outbreak in 2000 until 2010 (Fig 3A). After 2010, dengue cases began to increase rapidly until 2018, with a slight drop in cases in 2014. Over 5,000 infections were reported in 2000, 2002 and 2016. In contrast, in 2018, over 10,000 cases were reported. Dengue fever occurs primarily between July and November each year (Fig 3B). Therefore, monthly climate predictors were selected prior to July to predict the overall annual incidence. ![Fig 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/22/2020.09.21.20199190/F3.medium.gif) [Fig 3.](http://medrxiv.org/content/early/2020/09/22/2020.09.21.20199190/F3) Fig 3. Dengue incidence in Bangladesh for different years and months. (**A**) Annual dengue case trends between 2000 and 2018. A LOESS smoothing function is used to obtain a smooth line to represent the trend over the years. The shaded region shows the pointwise 95% confidence interval. (**B**) Month-wise average dengue cases from January 2000 to December 2018. The shaded area represents the 95% confidence interval. ### Climate variabilities Different combinations of temperature (monthly average, maximum or minimum), rainfall (monthly total or maximum) and sunshine duration were used to predict annual dengue incidence. The monthly minimum temperature in Bangladesh increased after January and continued to increase until June/July, with the highest temperatur*e* of *26*.4°*C* measured in 2010 and 2014 (Fig 4). However, large variations in temperature were observed between March and June. Rainfall increased between April and October (supplementary S2 Fig). The highest amount of rainfall usually occurs between May and August, with low levels of rainfall recorded in January, February and December. Sunshine duration in January to April and in December exhibits a decreasing trend over the years, an indication of warmer winter (Supplementary S1 Fig). Longer sunlight duration is mostly observed in the period from April to June. ![Fig 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/22/2020.09.21.20199190/F4.medium.gif) [Fig 4.](http://medrxiv.org/content/early/2020/09/22/2020.09.21.20199190/F4) Fig 4. Minimum temperatures (in °*C*) in Bangladesh between January 2000 and December 2018. A LOESS smoothing function is used to obtain a smooth line to represent the trend over the years. The shaded region represents the 95% confidence interval. Dots represent the average minimum temperature of a given month for a particular year. ### Model selection and annual dengue prediction To obtain an appropriate set of predictors for annual dengue prediction, we compared six models with different combinations of climate variables (Table 1). The best prediction model was determined using a two-stage model selection approach. In the first stage, Model 3 and Model 6, which belonged to either the low or the lowest *AIC**c* models, were chosen based on the criteria defined in Eq. (2). For details of the stepwise *AIC**c* results, please refer to supplementary S1 Table − S6 Table. In the second stage, *LOOCV* was conducted for the selected models. After *LOOCV* was performed, Model 3 was identified as the best prediction model, with the lowest mean squared error for the validation set, compared with Model 6 (0.29 vs 0.31; see Table 2). We used the best prediction model for further analysis of the impact of climate on dengue incidence. Note that the best fitting model (Model 6) also identified similar climate variables to the best prediction model. The best prediction model can be expressed as, View this table: [Table 1.](http://medrxiv.org/content/early/2020/09/22/2020.09.21.20199190/T1) Table 1. Comparison of candidate models based on different evaluation metrics. View this table: [Table 2.](http://medrxiv.org/content/early/2020/09/22/2020.09.21.20199190/T2) Table 2. The second stage of model selection. ![Formula][7] In Eq.(3) the values inside the brackets represent the standard error of the estimates. We evaluated the performance of the model using an *LOOCV* technique. The best prediction model successfully predicted the largest outbreak in 2018 as well as smaller outbreaks in five different years (2003, 2006, 2010, 2012 and 2014) (Fig 5). The estimated number of annual dengue cases for other years fell within a narrow range of the 95% bootstrap confidence interval (S9 Table). ![Fig 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/22/2020.09.21.20199190/F5.medium.gif) [Fig 5.](http://medrxiv.org/content/early/2020/09/22/2020.09.21.20199190/F5) Fig 5. Comparison of observed and predicted annual dengue cases. The 95% confidence intervals were estimated using a bootstrap estimation technique and the *LOOCV* estimates. ## Marginal effect of climate predictors To check the impact of each climate variable individually on annual dengue incidence, we further assessed the marginal effects of climate predictors. The optimal minimum temperature for mosquito population expansion is around 21–23°*C* (Fig 6). There was an upward trend of temperature from January until June. During this six-month period, the marginal effects of mean minimum temperature from January to March were positive. In contrast, the effects were negative from April to June. Thus, the results indicate that a turning point of marginal effects was located between 21 and 23°*C*. Starting from 2,500 predicted cases with a temperature of 23°*C*, the number of predicted cases gradually declined to below 1000 predicted cases in April if the temperature was increased by two degrees. The similar patterns were also evident in May and June. ![Fig 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/22/2020.09.21.20199190/F6.medium.gif) [Fig 6.](http://medrxiv.org/content/early/2020/09/22/2020.09.21.20199190/F6) Fig 6. Marginal effect of minimum temperature from January to June on annual dengue cases. *T*1 represents minimum temperature in January, *T*2 represents minimum temperature in February, etc. The shaded area denotes the 95% confidence interval of annual dengue cases at different values of minimum temperature, whereas the dotted vertical line represents mean minimum temperature for a month. The marginal effect here represents marginal effects at the mean (MEMs). Rainfall also had different effects depending on the time. In February, April and June, rainfall had a positive relationship with dengue incidence (Fig 7). In contrast, rainfall in a cooler winter period (January) had a negative association with dengue incidence. ![Fig 7.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/22/2020.09.21.20199190/F7.medium.gif) [Fig 7.](http://medrxiv.org/content/early/2020/09/22/2020.09.21.20199190/F7) Fig 7. Marginal effect of total rainfall on annual dengue incidence in January, February, April and June. *R*1 represents total rainfall in January, whereas *R*6 represents total rainfall in June, etc. The shaded area denotes the 95% confidence interval of the expected number of annual dengue cases for different values of rainfall whereas the dotted vertical line represents mean rainfall for a month. The marginal effect here represents marginal effects at the mean (MEMs). The mean sunshine duration in April and May were 7.4 and 6.5 hours, respectively and were negatively associated with annual dengue incidence (Fig 8). Comparing the magnitude of all climate predictors, temperature in June (T6), sunshine duration in May (S5) and rainfall in February (R2) had the strongest effects on annual dengue incidence. ![Fig 8.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/22/2020.09.21.20199190/F8.medium.gif) [Fig 8.](http://medrxiv.org/content/early/2020/09/22/2020.09.21.20199190/F8) Fig 8. Marginal effect of sunshine duration on annual dengue incidence in April and May. *S*4 represents sunshine duration in April whereas *S*5 represents sunshine duration in May. The shaded area denotes the 95% confidence interval of the expected number of annual dengue cases for different values of sunshine duration, whereas the dotted vertical line represents mean sunshine duration for a month. The marginal effect here represents marginal effects at the mean (MEMs). ## Discussion Climate change poses a great threat to global health, particularly for subtropical and tropical climate regions due to the expansion of dengue fever. Since 2010, Bangladesh has had an increasing total number of dengue cases during each seasonal epidemic. In 2018, the recorded number of confirmed and suspected cases was more than 10,000, including 26 confirmed deaths [3, 35]. Although many studies have attempted to estimate dengue incidence in Bangladesh using climate data, most of these studies focused on data collected prior to 2010 [15–17]. We developed a model to estimate the impact of various climate factors in the lead-up to dengue seasons, which typically occur during the monsoon season in Bangladesh. This is the first study to demonstrate that climate variability before dengue season can explain dengue expansions in Bangladesh in the past 20 years, suggesting that an early warning system can be built for this area. We found that some climate factors exerted opposing effects on the annual number of dengue cases, depending on the time of year. For example, minimum temperatures from January to March were positively associated with dengue cases, whereas a negative association was seen in the subsequent months from April to June, closer to the start of dengue season. This may be due to a complex dependency between the population dynamics of the dengue vector and the changing environment, such as the seasonal transition from winter to summer and the associated increasing temperatures. One study claimed that temperatures of 21.3 − 34.0°*C* are optimal for expansion of *Aedes aegypti* populations [21]. As the average daily minimum temperature is lower than 21°*C* before April and higher than 23°*C* during and after April, our results suggest that the optimal daily minimum temperature falls in the range of 21–23°*C*. Total rainfall in a winter month (January) was found to have a negative relationship with dengue cases, whereas a positive relationship was seen for later months mainly in summer such as April and June. Rainfall is thought to have both beneficial and harmful effects on mosquito population growth. Rainfall can provide standing water for mosquito breeding. However, an excessive amount of rainfall (e.g., heavy rainfall during monsoons or cyclones) has been commonly thought to be able to disrupt potential mosquito habitats [36]. Our study has identified a negative relationship between rainfall and dengue incidence in January. This may be because during this period the number of adult mosquitoes is low, meaning that rainfall has a larger negative impact by flushing away mosquito eggs than a positive impact due to creating habitats required by adult mosquitoes. A similar pattern of a negative association of pre-dengue-season rainfall with dengue cases has been seen in recent studies [23, 24, 37]. Additionally, sunshine duration was found to be closely linked to mosquito-related activities, such as frequency of mosquito bites [22]. However, sunshine duration has not been included in any prediction models thus far. Therefore, an evaluation of the increasing dengue incidence since 2010, with respect to these climate variables, is warranted. A shorter duration of sunshine is more favorable for dengue transmission. In general, mosquitoes are more active in darker environments, and there is a greater chance of dengue being transmitted during periods of less sunshine due to the increasing frequency of mosquito bites [39]. The marginal effect of sunshine duration in Fig 8 reveals that the shorter the duration of sunlight, the higher the number of dengue cases, supporting the biological characteristics of mosquito activity as described by [38]. A hour reduction in sunshine duration in April and May was predicted to result in a 3-fold increase in annual cases. These estimates are consistent with a previous study that found a negative association between sunshine duration and dengue incidence [39]. While this study has an implication, some limitations exist in this study. Meteorological data for 2019 were not disclosed at the time of this study; hence, 2019 data were excluded from the models. Secondly, the small number of outcome values (19) and a relatively large number of predictor variables might lead to over-fitting the results. Moreover, the epidemiological data includes both laboratory confirmed cases and probable case, hence, actual estimate might be affected due to under/over reporting biased. Because we aim to estimate the effects of preseasonal climate variability on annual incidence rate in this study. We did not plan to include monthly predictors within dengue season. Given the sample size is not huge, too many predictors may lead to over-fitting. In conclusion, our research offers a potential alert system by modeling annual dengue outbreaks before the season begins using climate variables. As an early warning system is required to improve public health and safety, the model we developed may improve disease control systems in Bangladesh. This research will aid our understanding of the effects of climate variability on dengue expansion, not only in Bangladesh but also in northern India and other Southeast Asian countries with similar climates and social-economic conditions. ## Data Availability All data used in this manuscript is available at the Directorate General of Health Services, Bangladesh and Bangladesh Meteorological Department upon formal request. [https://dghs.gov.bd/index.php/en/](https://dghs.gov.bd/index.php/en/) [http://live3.bmd.gov.bd/](http://live3.bmd.gov.bd/) ## Supplementary figures and table Mosquitoes are usually active in darker environments. With an increase in sunshine duration, their activity decreases, and the reverse is true for shorter sunshine duration [2]. Therefore, sunshine duration is expected to has an significant impact on mosquito-borne viral diseases like dengue. Sunshine duration measured in hours was found to have a decreasing trend from January to April and in December for different years between 2000 and 2018. However, it was almost stable from June to November. In addition, sunshine duration was shorter during the monsoon compared to the winter months from December to March (S1 Fig). Therefore, mosquitoes were more active during the monsoon, causing faster dengue transmission, which is reflected by the number of dengue cases in Fig 3 in the main text. ![Fig. S1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/22/2020.09.21.20199190/F9.medium.gif) [Fig. S1.](http://medrxiv.org/content/early/2020/09/22/2020.09.21.20199190/F9) Fig. S1. Sunshine duration (in hours) in Bangladesh between January 2000 and December 2018. A loess smoothing function is used to obtain a smooth line to represent trend over the years. The shaded region shows the confidence interval. Dots represent the average sunshine duration of a given month for a particular year. ![Fig. S2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/22/2020.09.21.20199190/F10.medium.gif) [Fig. S2.](http://medrxiv.org/content/early/2020/09/22/2020.09.21.20199190/F10) Fig. S2. Monthly total rainfall (in mm) in Bangladesh between January 2000 and December 2018. A loess smoothing function is used to obtain a smooth line to represent trend over years. The shaded region shows the confidence interval. Dots represent the average rainfall of a given month for a particular year. Rainfall is a critical meteorological factor that causes standing water in small storage areas, which are considered potential breeding sites for Aedes mosquitoes [1]. Rainfall was higher from May to August compared to other months of the year (S2 Fig). On the other hand, from December to March, the total amount of rain was minimal. However, the chance of rain water persisting in potential mosquito habitat during this time was higher due to light rain. Therefore, we found these three climate variables-minimum temperature, sunshine duration, and total rainfall-to be a good combination of climate variables capable of predicting dengue incidence optimally. We assumed that climate factors in the months prior to the dengue season had a significant impact on dengue transmission. Consequently, the minimum daily temperature and total daily rainfall from January to June and the daily sunshine duration from April to June were chosen as possible predictor variables for dengue prediction. View this table: [Table S1.](http://medrxiv.org/content/early/2020/09/22/2020.09.21.20199190/T3) Table S1. (**Model 1**) Step-by-step forward selection results of the generalized Poisson regression model in each step based on *AIC**c*.*ave*.*T**i*, *S**i* and *tot*.*R**i* represent mean temperature, sunshine duration and total rainfall in the *i**th* month. For each of the variable included in the model, the corresponding *AIC**c*, the leave-one-out mean squared error for the validation set (*MSE**V a*), the leave-one-out mean squared error for the training set (*MSE**T r*), the mean square ratio ![Graphic][8] and the factor *c*,![Graphic][9] were calculated. View this table: [Table S2.](http://medrxiv.org/content/early/2020/09/22/2020.09.21.20199190/T4) Table S2. (**Model 2**) Step-by-step forward selection results of the generalized Poisson regression model in each step based on *AIC**c*. *max*.*T**i*, *S**i* and *tot*.*R**i* represent maximum temperature, sunshine duration and total rainfall in the *i**th* month. For each of the variable included in the model, the corresponding *AIC**c*, the leave-one-out mean squared error for the validation set (*MSE**V a*), the leave-one-out mean squared error for the training set (*MSE**T r*), the mean square ratio ![Graphic][10] and the factor *c*,![Graphic][11] were calculated. View this table: [Table S3.](http://medrxiv.org/content/early/2020/09/22/2020.09.21.20199190/T5) Table S3. (**Model 3**) Step-by-step forward selection results of the generalized Poisson regression model in each step based on *AIC**c*. *T**i*, *S**i* and *R**i* represent minimum temperature, sunshine duration and total rainfall in the *i**th* month. For each of the variable included in the model, the corresponding *AIC**c*, the leave-one-out mean squared error for the validation set (*MSE**V a*), the leave-one-out mean squared error for the training set (*MSE**T r*), the mean square ratio ![Graphic][12] and the factor *c*,![Graphic][13] were calculated. View this table: [Table S4.](http://medrxiv.org/content/early/2020/09/22/2020.09.21.20199190/T6) Table S4. (**Model 4**) Step-by-step forward selection results of the generalized Poisson regression model in each step based on *AIC**c*. *ave*.*T**i*, *S**i* and *max*.*R**i* represent mean temperature, sunshine duration and maximum rainfall in the *i**th* month. For each of the variable included in the model, the corresponding *AIC**c*, the leave-one-out mean squared error for the validation set (*MSE**V a*), the leave-one-out mean squared error for the training set (*MSE**T r*), the mean square ratio ![Graphic][14] and the factor *c*,![Graphic][15]were calculated. View this table: [Table S5.](http://medrxiv.org/content/early/2020/09/22/2020.09.21.20199190/T7) Table S5. (**Model 5**) Step-by-step forward selection results of the generalized Poisson regression model in each step based on *AIC**c*. *max*.*T**i*, *S**i* and *max*.*R**i* represent maximum temperature, sunshine duration and maximum rainfall in the *i**th* month. For each of the variable included in the model, the corresponding *AIC**c*, the leave-one-out mean squared error for the validation set (*MSE**V a*), the leave-one-out mean squared error for the training set (*MSE**T r*), the mean square ratio ![Graphic][16] and the factor *c*, ![Graphic][17] were calculated. View this table: [Table S6.](http://medrxiv.org/content/early/2020/09/22/2020.09.21.20199190/T8) Table S6. (**Model 6**) Step-by-step forward selection results of the generalized Poisson regression model in each step based on *AIC**c*. *min*.*T**i*, *S**i* and *max*.*R**i* represent minimum temperature, sunshine duration and maximum rainfall in the *i**th* month. For each of the variable included in the model, the corresponding *AIC**c*, the leave-one-out mean squared error for the validation set (*MSE**V a*), the leave-one-out mean squared error for the training set (*MSE**T r*), the mean square ratio ![Graphic][18] and the factor *c*, ![Graphic][19] were calculated. View this table: [Table S7.](http://medrxiv.org/content/early/2020/09/22/2020.09.21.20199190/T9) Table S7. Leave-one-out cross-validation (loocv) results for **Model 6** by omitting the *j**th* year in the *j**th* iteration, where *j* = 1, …, 19, and *j* = 1 indicates the year 2000, *j* = 2 indicates 2001 and so on. Italic font represents the predicted annual dengue cases when the *j**th* year is removed. View this table: [Table S8.](http://medrxiv.org/content/early/2020/09/22/2020.09.21.20199190/T10) Table S8. Leave-one-out cross-validation (loocv) results for **Model 3** by omitting the *j**th* year in the *j**th* iteration, where *j* = 1, …, 19, and *j* = 1 indicates the year 2000, *j* = 2 indicates 2001 and so on. Italic font represents the predicted annual dengue cases when the *j**th* year is removed. View this table: [Table S9.](http://medrxiv.org/content/early/2020/09/22/2020.09.21.20199190/T11) Table S9. Comparison between observed and predicted annual dengue cases in Bangladesh between 2000 and 2018. The lower bound and the upper bound represents the lower and upper limit of 95% bootstrap confidence interval respectively for the predicted value. ## Acknowledgments We thank Dr. Iqbal Ansary Khan, Principal Scientific Officer, and Head of Medical Social Science, Institute of Epidemiology, Disease Control and Research, Directorate General of Health Services, Dhaka, for generously providing dengue surveillance data. The authors are indebted to City University of Hong Kong for providing excellent research facilities. The authors also acknowledge support from grants funded by City University of Hong Kong [#7200573 and #9610416]. * Received September 21, 2020. * Revision received September 21, 2020. * Accepted September 22, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NoDerivs 4.0 International), CC BY-ND 4.0, as described at [http://creativecommons.org/licenses/by-nd/4.0/](http://creativecommons.org/licenses/by-nd/4.0/) ## References 1. 1.Guo C, Zhou Z, Wen Z, Liu Y, Zeng C, Xiao D, et al. Global Epidemiology of Dengue Outbreaks in 1990^a€”2015: A Systematic Review and Meta-Analysis. Frontiers in Cellular and Infection Microbiology. 2017;7:317. doi:10.3389/fcimb.2017.00317. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fcimb.2017.00317&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28748176&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.21.20199190.atom) 2. 2.Ferreira GLC. Global Dengue Epidemiology Trend; 2012. 3. 3.(WHO) WHO. Dengue and severe dengue; 2020. [https://www.who.int/news-room/fact-sheets/detail/dengue-and-severe-dengue](https://www.who.int/news-room/fact-sheets/detail/dengue-and-severe-dengue). 4. 4.Stanaway JD, Shepard DS, Undurraga EA, Halasa YA, Coffeng LE, Brady OJ, et al. The global burden of dengue: an analysis from the Global Burden of Disease Study 2013. The Lancet Infectious Diseases. 2016;16(6):712–723. doi:10.1016/S1473-3099(16)00026-8. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(16)00026-8&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26874619&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.21.20199190.atom) 5. 5.Pasin CI, Elizabeth Halloran M, Gilbert PB, Langevin E, Leon Ochiai R, Pitisuttithum P, et al. Periods of high dengue transmission defined by rainfall do not impact efficacy of dengue vaccine in regions of endemic disease. Plos One. 2018;doi:10.1371/journal.pone.0207878. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0207878&link_type=DOI) 6. 6.Loo YY, Billa L, Singh A. Effect of climate change on seasonal monsoon in Asia and its impact on the variability of monsoon rainfall in Southeast Asia. Geoscience Frontiers. 2015;6(6):817–823. doi:10.1016/j.gsf.2014.02.009. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.gsf.2014.02.009&link_type=DOI) 7. 7.Louis VR, Montenegro Quin∼onez CA, Kusumawathie P, Palihawadana P, Janaki S, Tozan Y, et al. Characteristics of and factors associated with dengue vector breeding sites in the City of Colombo, Sri Lanka. Pathogens and Global Health. 2016;110(2):79–86. doi:10.1080/20477724.2016.1175158. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/20477724.2016.1175158&link_type=DOI) 8. 8.Getachew D, Tekie H, Gebre-Michael T, Balkew M, Mesfin A. Breeding sites of Aedes aegypti: Potential dengue vectors in Dawa, east Ethiopia. Interdisciplinary Perspectives on Infectious Diseases. 2015;2015. doi:10.1155/2015/706276. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1155/2015/706276&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26435712&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.21.20199190.atom) 9. 9.Evelyn N, Murray A, Quam MB, Wilder-Smith A. Epidemiology of dengue: past, present and future prospects. Clinical Epidemiology. 2013;doi:10.2147/CLEP.S34440. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2147/CLEP.S34440&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23990732&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.21.20199190.atom) 10. 10.World Health Organization (WHO). Ten threats to global health in 2019; 2020. [https://www.who.int/vietnam/news/feature-stories/detail/ten-threats-to-global-health-in-2019](https://www.who.int/vietnam/news/feature-stories/detail/ten-threats-to-global-health-in-2019). 11. 11.Russell PK, Buescher EL, McCown JM, Ord∼onez J. Recovery of dengue viruses from patients during epidemics in Puerto Rico and East Pakistan. The American Journal of Tropical Medicine and Hygiene. 1966;15(4):573–579. doi:10.4269/ajtmh.1966.15.573. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoidHJvcG1lZCI7czo1OiJyZXNpZCI7czo4OiIxNS80LzU3MyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA5LzIyLzIwMjAuMDkuMjEuMjAxOTkxOTAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 12. 12.Rahman M, Rahman K, Siddque AK, Shoma S, Kamal AHM, Ali KS, et al. First Outbreak of Dengue Hemorrhagic Fever, Bangladesh. Emerging Infectious Diseases. 2002;8(7). 13. 13.Institute of Epidemiology DC, (IEDCR) R. Web-based Dengue Surveillance; 2020. [https://www.iedcr.gov.bd/index.php/surveillance](https://www.iedcr.gov.bd/index.php/surveillance). 14. 14.Hsan K, Hossain MM, Sarwar MS, Wilder-Smith A, Gozal D. Unprecedented rise in dengue outbreaks in Bangladesh; 2019. 15. 15.Sharmin S, Glass K, Viennet E, Harley D. Geostatistical mapping of the seasonal spread of under-reported dengue cases in Bangladesh. PLoS Neglected Tropical Diseases. 2018;12(11):1–13. doi:10.1371/journal.pntd.0006947. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pntd.0006947&link_type=DOI) 16. 16.Banu S, Hu W, Guo Y, Hurst C, Tong S. Projecting the impact of climate change on dengue transmission in Dhaka, Bangladesh. Environment International. 2014;63:137–142. doi:10.1016/j.envint.2013.11.002. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.envint.2013.11.002&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24291765&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.21.20199190.atom) 17. 17.Hashizume M, Dewan AM, Sunahara T, Rahman MZ, Yamamoto T. Hydroclimatological variability and dengue transmission in Dhaka, Bangladesh: A time-series study. BMC Infectious Diseases. 2012;12. doi:10.1186/1471-2334-12-98. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2334-12-98&link_type=DOI) 18. 18.Lai YH. The climatic factors affecting dengue fever outbreaks in southern Taiwan: An application of symbolic data analysis. BioMedical Engineering Online. 2018;17(2):1–14. doi:10.1186/s12938-018-0575-4. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12938-018-0575-4&link_type=DOI) 19. 19.Lu L, Lin H, Tian L, Yang W, Sun J, Liu Q. Time series analysis of dengue fever and weather in Guangzhou, China. BMC Public Health. 2009;9:1–5. doi:10.1186/1471-2458-9-395. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2458-9-1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19121216&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.21.20199190.atom) 20. 20.Gu H, Leung RKK, Jing Q, Zhang W, Yang Z, Lu J, et al. Meteorological factors for dengue fever control and prevention in South China. International Journal of Environmental Research and Public Health. 2016;13(9):1–12. doi:10.3390/ijerph13090867. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/ijerph13090867&link_type=DOI) 21. 21.Ryan Id SJ, Carlson CJ, Mordecai EA, Johnson LR. Global expansion and redistribution of Aedes-borne virus transmission risk with climate change. PloS Neglected Tropical Diseases. 2019;doi:[https://doi.org/10.1371/journal.pntd.0007213](https://doi.org/10.1371/journal.pntd.0007213). 22. 22.Vu HH, Okumura J, Hashizume M, Tran DN, Yamamoto T. Regional differences in the growing incidence of dengue fever in Vietnam explained by weather variability. Tropical Medicine and Health. 2014;42(1):25–33. doi:10.2149/tmh.2013-24. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2149/tmh.2013-24&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24808744&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.21.20199190.atom) 23. 23.Yuan HY, Wen TH, Kung YH, Tsou HH, Chen CH, Chen LW, et al. Prediction of annual dengue incidence by hydro-climatic extremes for southern Taiwan. International Journal of Biometeorology. 2019;63(2):259–268. doi:10.1007/s00484-018-01659-w. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00484-018-01659-w&link_type=DOI) 24. 24.Yuan HY, Liang J, Lin PS, Sucipto K, Tsegaye MM, Wen TH, et al. The effects of seasonal climate variability on dengue annual incidence in Hong Kong: A modelling study. Scientific Reports. 2020;10(1):1–10. doi:10.1038/s41598-020-60309-7. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-020-60309-7&link_type=DOI) 25. 25.(WHO) WHO. South-East Asia; 2020. [https://www.who.int/southeastasia](https://www.who.int/southeastasia). 26. 26.Worldatlas. Where Is Bangladesh?; 2020. [https://www.worldatlas.com/as/bd/where-is-bangladesh.html](https://www.worldatlas.com/as/bd/where-is-bangladesh.html). 27. 27.Wikipedia. Tropic of Cancer; 020. [https://en.wikipedia.org/wiki/Tropic\_of\_Cancer](https://en.wikipedia.org/wiki/Tropic_of_Cancer). 28. 28.Banglapedia. Rainfall; 2020. [http://en.banglapedia.org/index.php?title=Rainfall](http://en.banglapedia.org/index.php?title=Rainfall). 29. 29.Banglapedia. Bangladesh Geography; 2020. [http://en.banglapedia.org/index.php?title=Bangladesh\_Geography](http://en.banglapedia.org/index.php?title=Bangladesh_Geography). 30. 30.Directorate General of Health Services. National Guideline for Clinical Management of Dengue Syndrome. 4th ed. Dhaka: National Malaria Elimination & Aedes Transmitted Diseases Control Program Disease Control Unit; 2020. 31. 31.Hurvich CM, Tsai CL. Regression and Time Series Model Selection in Small Samples. Biometrika. 1989;76(2):297. doi:10.2307/2336663. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/biomet/76.2.297&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1989AE13400012&link_type=ISI) 32. 32.Fernihough A. mfx: Marginal Effects, Odds Ratios and Incidence Rate Ratios for GLMs; 2019. = [https://CRAN.R-project.org/package=mfx](https://CRAN.R-project.org/package=mfx). 33. 33.LA∼1/4decke D. ggeffects: Tidy Data Frames of Marginal Effects from Regression Models. Journal of Open Source Software. 2018;3(26):772. doi:10.21105/joss.00772. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.21105/joss.00772&link_type=DOI) 34. 34.Zou KH, O’Malley AJ, Mauri L. Receiver-operating characteristic analysis for evaluating diagnostic tests and predictive models. Circulation. 2007;115(5):654–657. doi:10.1161/CIRCULATIONAHA.105.594929. [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MTQ6ImNpcmN1bGF0aW9uYWhhIjtzOjU6InJlc2lkIjtzOjk6IjExNS81LzY1NCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA5LzIyLzIwMjAuMDkuMjEuMjAxOTkxOTAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 35. 35.DGHS. Daily Dengue Status Report; 2020. [https://dghs.gov.bd/index.php/bd/home/5200-daily-dengue-status-report](https://dghs.gov.bd/index.php/bd/home/5200-daily-dengue-status-report). 36. 36.Benedum CM, Seidahmed OME, Eltahir EAB, Markuzon N. Statistical modeling of the effect of rainfall flushing on dengue transmission in Singapore. PLOS Neglected Tropical Diseases. 2018;12(12). doi:10.1371/journal.pntd.0006935. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pntd.0006935&link_type=DOI) 37. 37.Lowe R, Gasparrini A, Van Meerbeeck CJ, Lippi CA, Mahon R, Trotman AR, et al. Nonlinear and delayed impacts of climate on dengue risk in Barbados: A modelling study. PLOS Medicine. 2018;15(7). doi:10.1371/journal.pmed.1002613. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pmed.1002613&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30016319&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.21.20199190.atom) 38. 38.Kim YM, Park JW, Cheong HK. Estimated effect of climatic variables on the transmission of Plasmodium vivax malaria in the Republic of Korea. Environmental Health Perspectives. 2012;120(9):1314–1319. doi:10.1289/ehp.1104577. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1289/ehp.1104577&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22711788&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.21.20199190.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000308786900031&link_type=ISI) 39. 39.Pham HV, Doan HTM, Phan TTT, Tran Minh NN. Ecological factors associated with dengue fever in a Central Highlands province, Vietnam. BMC Infectious Diseases. 2011;11(1):172. doi:10.1186/1471-2334-11-172. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2334-11-172&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21679398&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.21.20199190.atom) ## References 1. 1.Getachew D, Tekie H, Gebre-Michael T, Balkew M, Mesfin A. Breeding sites of aedes aegypti: Potential dengue vectors in dire Dawa, east Ethiopia. Interdisciplinary Perspectives on Infectious Diseases. 2015;2015. doi:10.1155/2015/706276. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1155/2015/706276&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26435712&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.21.20199190.atom) 2. 2.Kim YM, Park JW, Cheong HK. Estimated effect of climatic variables on the transmission of plasmodium vivax malaria in the republic of Korea. Environmental Health Perspectives. 2012;120(9):1314–1319. doi:10.1289/ehp.1104577. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1289/ehp.1104577&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22711788&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.21.20199190.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000308786900031&link_type=ISI) [1]: /embed/graphic-2.gif [2]: F2/embed/inline-graphic-1.gif [3]: /embed/graphic-4.gif [4]: /embed/inline-graphic-2.gif [5]: /embed/inline-graphic-3.gif [6]: /embed/inline-graphic-4.gif [7]: /embed/graphic-9.gif [8]: T3/embed/inline-graphic-5.gif [9]: T3/embed/inline-graphic-6.gif [10]: T4/embed/inline-graphic-7.gif [11]: T4/embed/inline-graphic-8.gif [12]: T5/embed/inline-graphic-9.gif [13]: T5/embed/inline-graphic-10.gif [14]: T6/embed/inline-graphic-11.gif [15]: T6/embed/inline-graphic-12.gif [16]: T7/embed/inline-graphic-13.gif [17]: T7/embed/inline-graphic-14.gif [18]: T8/embed/inline-graphic-15.gif [19]: T8/embed/inline-graphic-16.gif