 Research
 Open Access
 Published:
Describing interaction effect between lagged rainfalls on malaria: an epidemiological study in south–west China
Malaria Journalvolume 16, Article number: 53 (2017)
Abstract
Background
When discussing the relationship between meteorological factors and malaria, previous studies mainly focus on the interaction between different climatic factors, while the possible interaction within one particular climatic predictor at different lag periods has been largely neglected. In this study, this issue was investigated by exploring the interaction of lagged rainfalls and its impact on malaria epidemics, which is a typical example of those meteorological variables.
Methods
The weekly data of malaria cases and three climatic variables of 30 counties in southwest China from 2004 to 2009 were analysed with the varying coefficientdistributed lag nonlinear model. The correlation patterns of the 6th, 9th and 12th week lags would vary over different rainfall levels at the 4thweek lag.
Results
The nonlinear patterns for rainfall at different rainfall levels are distinct from each other. In the low rainfall level at the 4th week lag, the increasing rainfall may promote the transmission of malaria. However, for the high rainfall level at the 4th week lag, evidence shows that the excessive rainfall decreases the risk of malaria.
Conclusion
This study reports for the first time that the interaction effect between lagged rainfalls on malaria exists, and highlights the importance of integrating the interaction between lagged predictors in relevant studies, which could help to better understand and predict malaria transmission.
Background
Malaria is one of the most severe mosquitoborne infectious diseases that threaten human health around the globe [1]. Since the primary vector for malaria is Anopheles, the dynamics and distribution of malaria are closely correlated with meteorological conditions [2]. From a biological perspective, climatic factors have a profound effect on malaria incidence as they affect both the development of mosquitoes and malaria parasites within mosquitoes [3, 4]. For instance, plentiful rainfall and increased humidity offer mosquitoes suitable sites to breed, resulting in increasing number of mosquitoes. Moreover, appropriate temperatures promote malaria prevalence by enhancing the mosquito’s growth and biting rate, as well as increasing the viability and development rate of parasites within vectors [5].
When exploring the meteorological effect on malaria incidence, special attention should be paid to two particular issues, which are the lag and nonlinear patterns [6, 7]. In some previous studies, lag time was assumed to be single fixed [8–11]. This is unreasonable, especially when describing the relationship between climatic variables and malaria risk at the large population level. Biologically speaking, there are at least two stages that should be considered for the lag effect, such as the development of mosquitoes and the incubation of parasites within the mosquito. The lag time of each stage may vary indefinitely based on climatic conditions, thus generating a smoothly varied lag distribution at population level. The nonlinear correlation between rainfall and malaria incidence has been acknowledged and validated both experimentally and epidemiologically in a series of existing studies [12, 13]. It is proposed that similar nonlinear effect may also exist for temperature [14–16]. Therefore, both lag and nonlinear patterns should be considered in the model. For this purpose, the distributed lag nonlinear model proves a valuable and effective method [17].
When discussing the relationship between meteorological factors and malaria incidence, however, in most previous epidemiological studies, the effect of betweenlag interaction has long been overlooked. Unlike the interaction between different exposure variables where the various climatic factors are of the same timeperiod and simultaneously affect malaria incidence, the betweenlag interaction is defined as the interaction between one covariate at different lag time, such as the interaction effect of the rainfalls four and five weeks previously on malaria risk for the current week. The concept of lagged interaction was first developed by Heaton in 2014, as he reported the corresponding statistical methods to examine the relationship between heat exposure and mortality [18]. When considering the betweenlag interaction, the total effect of the climatic variables is not simply an accumulation of lagged effects, as is conventionally supposed. This is explained by the fact that the number of mosquitoes at week t are substantially dependent on the rainfall at week t − 1. More specifically, adequate rainfall at week t − 1 provides mosquitoes with abundant breeding sites and promotes their development, which leads to an increased number of mosquitoes at week t, and consequently a greater risk of malaria transmission.
This study aims to investigate the interaction effect between climatic factors at different lag periods on malaria risk, exemplified by rainfall. So far, it has not been directly reported in previous studies. Specifically, with weekly data during the period of 2004–2009 among 30 counties in southwest China, a varying coefficient distributed lag nonlinear model was applied to investigate the association between rainfall and malaria cases. The correlation pattern between rainfall and malaria incidence was set to change depending on the three levels of rainfall at the fourth week lag. The results can help researchers to better understand the complex relationship between climatic factors and malaria transmission, and develop potential better prediction models.
Methods
Study sites
The southwestern region of China has been severely threatened by malaria in the last century [19]. After decades of continuous effort, malaria in these provinces has been brought under control. However, due to the effect of global warming, as well as the complex meteorological condition in this region (most counties are of subtropical climate; a few southern counties in the tropical region) [20], malaria poses a potential threat to the health of its populace.
The southwest region of China encompasses a large area from 21°14′N to 34°31′N and 97°35′E to 110°19′E. It covers primarily four provinces of Sichuan, Chongqing, Yunnan, and Guizhou with 483 counties in total (countylevel cities and districts). This region has a population of 189,977,077 (results from the sixth national census in 2010) and occupies an area of 1,137,570 sq km. Malaria data covered 483 counties, while only 131 counties had daily meteorological records. Chosen from those counties where malaria and climatic data were available, 30 counties with the highest average annual incidence were chosen for this study. The relevant principles of sample selection have been described in a previous study [21]. Additional file 1 shows the 483 counties and the selected 30 counties in southwest China.
Data description
The weekly meteorological data from July 2003 to December 2009 was collected from the Chinese Meteorological Data Sharing Service System [22], in which the mean temperature (°C), rainfall (mm) and relative humidity (%) had been recorded. There are 131 meteorological monitoring stations among the 483 counties in the southwest region, and those that are relevant to the counties with high malaria incidence were used.
As daily malaria records would bring many zero counts and jeopardize the stability of the parameter estimation, weekly case reports were used in this study and collected through the Chinese Information System for Infectious Diseases Control and Prevention (CISIDCP) from 2004 to 2009 among the 30 counties in southwest China [23]. Although both malaria subtypes (Plasmodium vivax and Plasmodium falciparum) could be the potential cause of the cases reported, most records did not include the type of parasites. The population data associated with the selected counties were obtained from the National Bureau of Statistics of China, from 2004 to 2009.
Basic distributed lag nonlinear model (DLNM)
The methodology of distributed lag nonlinear model (DLNM) is used to describe the dependencies that are both nonlinear and delayed [17]. Because malaria cases are ordinary count data, the association between the expected number of malaria cases E(Y _{ it }) at week t in county i and climatic variables in the previous weeks was modelled by the Poisson regression,
where d _{ it } denotes the population in county i at week t; β _{ i0} denotes the intercept effect for county i. x _{ it,r }, x _{ it,h } and \(\, x_{{it,T_{m} }}\) are the weekly meteorological variables for county i at week t, representing the rainfall, relative humidity and mean temperature, respectively.
The lag ranges for meteorological factors were determined based on relevant biological knowledge. Considering the lag effect, it would be logical to assume that cases in a specific week will be affected by climatic factors several weeks before. Consequently, Model (1) could be used to estimate the cumulative contributions across the entire lag range, rather than on a single fixed lag time. The lag ranges were determined by referring to those biological factors and empirical results from previous studies [24]. The period from the fourth to the 15th week was used as the lag range for weekly rainfall and mean relative humidity, while for weekly mean temperature it was from the third to the tenth week [24, 25].
There are two basis functions included in Model (1) to represent the nonlinear and lag effects. Taking rainfall as an example, the first is f(x _{ i(t − l),r }, β _{ rl }), which describes the nonlinear effect of rainfall that happened l weeks before. It could be interpreted by many functional forms, for instance, the polynomial function. As for the other function, its purpose is to constrain the parameter β _{ rl }, thus to refrain from the high collinearity caused by the significant correlation between rainfall on consecutive weeks. With the introduction of the constraining function, a reduction of the noise in the unconstrained distributed lag model could be achieved with less bias [7]. Subsequently, the secondorder natural cubic spline was applied for the investigation of both basis functions in Model (1) due to the fact that the meteorological variables are unimodal [26], as well as the requirement for parsimony.
Correlations between climatic factors and malaria incidences are not equal from place to place, and that the correlation in one county would be comparatively larger than that within two counties, in most circumstances. The inequality was caused by some unmeasured (or even unmeasurable) countyspecific variables. To deal with the potential confounding, β _{ i0} is modelled as a multilevel random intercept, obeying a normal distribution that β _{ i0} ∼ N(β _{0}, δ _{0} ^{2} ). β _{0} is the average intercept of all counties, while δ _{0} ^{2} represents the variation of countyspecific intercepts around β _{0}.
Varying coefficient distributed lag nonlinear model
The lagged interaction was excluded in Model (1) due to the hypothesis that the rainfall at week t has the same effect at different levels of rainfall at week t − k. However, the effect of rainfall at week t may also depend on the level of rainfall at week t − k. In this regard, a varying coefficient model [27] was used to examine the dependencies of the rainfall effect at week t on the rainfall level at week t − k and to investigate the betweenlag interaction.
The rainfall at the fourth week lag was set as the stratification variable in this study, mainly considering the lag range. As x _{ i(t − 4),r } denotes the rainfall at the fourth week lag, x _{ i(t4),r } should, in a certain degree, influence the lag nonlinear pattern of rainfall at other lag weeks, providing the hypothesis is correct that the lagged interaction exists and the effect of rainfall at other lag weeks is indeed dependent on the rainfall level at the fourth week lag.
To investigate the possible influence of the lagged interaction, all x _{ i(t − 4),r } were divided into three quantile groups (33.3 and 66.6% percentiles). The three groups were denoted as R _{ i(t − 4),r0}, R _{ i(t − 4),r1} and R _{ i(t − 4),r2}, which represents the x _{ i(t − 4),r } at the low, medium and high level of rainfall at the fourth week lag, respectively. Model (1) is then adjusted to embody these changes, and the modifications are shown below:
Compared to Model (1), Model (2) is different in two major aspects. The first is that β _{ rl } has now become β _{ rl }(R _{ i(t − 4),rg }), reflecting the fact that the coefficient β _{ rl } is now varying over R _{ i(t4),rg }. Consequently, the effect of rainfall at other lag weeks is now dependent on the relevant rainfall level at the fourth week lag. Besides the lag nonlinear pattern, the modified model can also be used to interrogate the lagged interaction that the rainfall effect of other lag weeks changes over different levels of rainfall at the fourth week lag. If the lagged interaction indeed exists, the lag nonlinear patterns for rainfall at the three levels should differ from each other.
The second difference is that the coefficient R _{ i(t − 4),rg } is now included in the intercept. Like the ordinary categorical predictor, R _{ i(t − 4),r0} at the lowest level was used as the reference group, and α _{1} and α _{2} represent the differential effects for R _{ i(t − 4),r1} and R _{ i(t − 4),r2}, respectively. An assumption was made in Model (2) that the mean temperature and relative humidity do not interact with the rainfall. The reference values of all climatic factors are all set at zero.
The analysis may be sensitive to the choice of lag ranges, and therefore the situation where the fourth to 14th weeks were selected as lag range for rainfall instead of the fourth to 15th weeks was also investigated, and the results were robust and showed no significant change. All the analysis was performed with R, which is an open statistical software [28]. Specifically, the package lme4 [29] was loaded to estimate the parameters. The data above have been used in previous studies by the same research group to explore other relationships between meteorological factors and malaria incidences [21, 30].
Results
Descriptive analysis
Among the selected 30 counties in southwest China, 21,944 malaria cases were reported in total over the 6 years from 2004 to 2009, and the descriptive analysis of the collected data could be found in Additional file 2.
The three rainfall ranges at the fourth week lag are (0.0, 1.6 mm), (1.6, 20.0 mm), (20.0, 310.9 mm), and their sample sizes are 3155, 3106 and 3129, respectively. The slight difference between the sample sizes is ascribed to the fact that most weekly rainfall values are concentrated at the 33.3 and 66.6 percentile cutpoints. Nonetheless, these values are eventually grouped into the corresponsive levels and cause the inequality of sample numbers. In Fig. 1, the meteorological variables between different rainfall levels have been compared, and pronounced trends were found: the mean temperature, rainfall and relative humidity were all positively correlated with the rainfall levels at the fourth week lag. Conclusively, the annualized average incidences are 9.73, 14.18 and 13.36 per 100,000 for low, medium and high rainfall levels, respectively.
Varying coefficient distributed lag nonlinear model
Figure 2 demonstrates the estimated lagged nonlinear relationships between rainfall and malaria incidence in the exposure dimension. The results are presented with nine representative panels, as there are three lag conditions and each condition was analysed under three levels of rainfall at the fourth week lag. The Yaxis represents the logarithm value of the relative risk ratio in comparison with the reference value at rainfall 0 mm. In general, the rainfall levels at the fourth week lag affect the nonlinear patterns between rainfall and malaria incidence under each lag time. In the low rainfall level at the fourth week lag, rainfall promotes malaria incidence. On the contrary, in the high rainfall level at the fourth week lag, the risk of malaria incidence would decrease as the rainfall increases. (See figure on previous page).
For each lag time, the nonlinear patterns between rainfall and malaria incidence are distinctively different across the three rainfall levels. In the high rainfall level at the fourth week lag, the magnitude and direction of the association would change according to the rainfall values. More specifically, in the high level, the logarithmic value of relative risk (logRR) increases slightly at first and reaches the maximum at approximate 75 mm, but then it starts to decline sharply. In the range 0–150 mm, rainfall is positively associated with malaria incidence. Afterwards, the correlation between rainfall and malaria risk becomes negative. The same trend is observed for the 12th week lag in the medium level (panel h of Fig. 2). Secondly, in the medium rainfall level at the fourth week lag (except panel h), the curve shows a slowly downward trend. However, this correlation is not statistically significant, as the value of 0 is included in the 95% confidence interval of logRR.
In the low rainfall level at the fourth week lag, rainfall is positively associated with malaria incidence, and the logRR would increase sharply when the rainfall increases. The ninth week lag has the greatest impact on malaria risk at the low and high rainfall levels. The effect of the sixth week lag ranks second, which is greater than that of the 12th week lag.
Figure 3 shows the estimated lagged nonlinear relationships between rainfall and malaria incidence in the lag dimension. The results are presented with three rainfall values at the 25, 50 and 75% percentiles of weekly rainfall from 2004 to 2009, which are 0.2, 15.5 and 30.8 mm, respectively. As the three representative rainfall values were within the three rainfall level ranges (0.0, 1.6 mm), (1.6, 20.0 mm) and (20.0, 310.9 mm), respectively, Fig. 3 also demonstrates the overall relationship between rainfall and malaria incidence for each rainfall level at the fourth week lag. Specifically, Fig. 3a shows the correlation at the low rainfall level, Fig. 3b the medium rainfall level and Fig. 3c the high rainfall level. Yaxis represents the logarithm value of the relative risk ratio in comparison with the reference value at rainfall 0 mm. The three panels present distinct nonlinear patterns between rainfall and malaria incidence at different rainfall values. When the rainfall is at 30.8 mm, the corresponding distributed lag curve shows an apparent inversed Ushape, which goes up first until peaking at the ninth week, then starts going downwards and ends up with nonsignificant correlation from the 13th week. When the rainfall is at 15.5 and 0.2 mm, both curves show similar trends: the correlations are nonsignificant at first and then start to increase substantially during the 11–15th weeks. Finally, by comparing the three panels, it is observed that the logRR increases along with the rainfall values, and the increasing rate is much larger at the lower rainfall. Specifically, the increasing rate at 15.5 mm is nearly 100 times of that at 0.2 mm, while the difference becomes minimal when comparing that under 15.5 and 30.8 mm.
Discussion
The interaction between meteorological factors is important in malaria transmission, as they are closely associated with vector abundance and survival as well as parasite maturation [24]. Exploring the interaction between climatic variables on malaria incidence can help to better understand the relationship between meteorological factors and malaria incidence [31]. The interaction between exposure predictors is common in existing studies, but so far there is no report dedicated to the lagged interaction between climatic factors in the process of malaria transmission, which may also play a crucial role in malaria epidemics. Particularly, the lagged interaction effect on malaria incidence examined in this study was the interaction between rainfall at the 4th week lag and that at the 6th, 9th and 12th week lags.
The results indicate that the rainfall at the 4th week lag affects the correlation between malaria incidence and rainfall at the other lag weeks, implying the interaction effect between lagged rainfalls on malaria. When in the low rainfall level at the fourth week lag, the malaria risk increases along with the increase of rainfall, suggesting that the increasing rainfall promotes malaria transmission when rainfall is low at the 4th week lag. In contrast, excessive rainfall decreases the risk of malaria when rainfall is high at the 4th week lag, which can be observed in panel c, f, i of Fig. 2. These results can be explained by malaria dynamics [19]: rainfall elevates the environmental humidity and brings about many temporary puddles, simultaneously increasing the number of mosquito breeding sites and enhancing mosquito survival. However, excessive rainfall and accumulation of surface water in complicated terrain would potentially destroy mosquitobreeding sites, thus reducing the mosquito density. Abundant rainfall may also prevent people from working outdoors, resulting in lower chances of people being bitten by mosquitoes and consequently, decreasing malaria incidence. Specifically, when the rainfall level at the 4th week lag is low, the greater rainfall at week t would relieve the effect of insufficient rainfall so that rainfall offers more breeding habitats to mosquitoes and increases their number, resulting in the increased risk of malaria incidence. In contrast, when the rainfall level at the fourth week lag is high, abundant rainfall at week t would exacerbate the effect of excessive rainfall, resulting in mosquito breeding sites being destroyed and people reducing their outdoor activities. The effect of high rainfall at week t would be attenuated or even become negative.
It is also observed that the lagged effect of rainfall on malaria incidence was greatest at the ninth lag week, compared to the 6th and 12th weeks. This is biologically acceptable as the effect of rainfall occurring in the current week or too long before is negligible on malaria incidence.
Greater rainfall brings a higher relative risk but a shorter lag for malaria cases, which can be observed from Fig. 3. Specifically, in the low and medium levels of rainfall at the 4th week lag, rainfall starts to be significantly associated with malaria incidence at the 11th week. In high rainfall level at the fourth week lag, the distributed lag curve shows a significant correlation from the 4th to 13th week. Compared with the high level of rainfall at the fourth week lag, rainfall in low levels is associated with delayed malaria risk. The results are consistent with a previous study [32]. These may be a result from the previously mentioned malaria dynamics that rainfall could provide fitted habitats for mosquitoes to breed, thus shortening their life cycles and accelerating the spreading of malaria [33]. Despite the fact that the relative risk for malaria cases is positively correlated with rainfall, the increase in the relative risk is more drastic when rainfall is low, while it becomes minimal when rainfall is high. This phenomenon may be explained by the saturation effect, where the contribution of increasing rainfall to the development of mosquito and parasite becomes negligible or even counterproductive.
Rainfall is selected as an example to reveal the lagged interaction on malaria incidence mainly for biological and epidemiological considerations. From an entomological perspective, rainfall affects most of the stages of the mosquito’s life cycle. For example, plentiful rainfall provides mosquitoes with aquatic breeding sites for their growth and reproduction [34]. While from an epidemiological view, the reported relationship between rainfall and malaria vary in the literature [24]. This is especially true in China, as several studies showed that rainfall was closely correlated with malaria incidence [16, 35, 36], while other studies denied the existence of such correlations [37]. Understanding the interaction effect between rainfall at different lag time on malaria incidence may help to better explore the relationship between rainfall and malaria incidence.
α _{1} and α _{2} are introduced to describe the main effects of the rainfall levels at the fourth week lag. As demonstrated in Fig. 1, the three groups of different rainfall levels at the fourth week lag do not have identical baseline distribution of climatic factors. Even under the same rainfall condition at the fourth week lag, the average effect of rainfall among the three groups should be distinctively different; α _{1} and α _{2} are consequently added as the average deviation, ensuring that the logRR values of all groups would be zero at the reference rainfall, and allowing the comparison of variations for rainfall.
The countyspecific random intercept model allows fitting a regression model to meteorological factors with the systematic unexplained variation among the 30 counties. As with other epidemiological literature on the relationship between meteorological variables and malaria incidence, the final results would be potentially interfered by some confounding factors. For instance, there may be different preventive measures with different enforcement strength that are deployed by an individual county to fight malaria, as well as some behavioural patterns, such as the utilization of nets of different types. The variance δ _{0} ^{2} of the countyspecific random intercept β _{ i0} represents the variation between counties that are not caused by the climatic predictor. The random intercept model has proven efficient in handling the potential bias [38].
There are a few limitations that should be acknowledged. First, the data quality may change over the 6 years. This primarily varied with time, and the best data quality was found in 2009. Second, only 30 counties with malaria prevalence were used in this study. However, this should not introduce significant inherent bias into this study. The range of malaria incidence varied from low to high values. Specifically, the annualized average incidence ranged from 348.2/100,000 to 1.1/100,000. In particular, it is evident (see Additional file 2) that the 30 counties included many lowincidence counties, such as Eshan with just 11 malaria cases in 6 years. Therefore, the selection method in this study should not undermine the credibility of this study. Third, like several existing studies [39, 40], the characteristics of P. vivax and P. falciparum were not analysed separately due to the lack of associated information in this study. As a result, vivax relapses may be mistaken as new infections caused by meteorological factors. The lag nonlinear patterns of the two malaria subtypes may be slightly different from each other. By investigating the potential bias, future studies might provide more details to elucidate the association between climatic factors and malaria incidence in southwest China.
Conclusions
Using weekly data of malaria cases and climatic variables during the period of 2004–2009 among 30 counties in southwest China, the interaction effect between rainfall at different lag time on malaria incidence was examined. As previous studies rarely accounted for the interaction between lagged climatic factors, this work highlights the importance of including the lagged interaction effect in the investigation of malaria incidence, which can provide supplementary evidence to understand and predict malaria transmission.
Abbreviations
 DLNM:

distributed lag nonlinear model
 logRR:

logarithmic value of relative risk
References
 1.
WHO. World malaria report 2011. Geneva: World Health Organization; 2015.
 2.
Zhou G, Minakawa N, Githeko AK, Yan G. Association between climate variability and malaria epidemics in the East African highlands. Proc Natl Acad Sci USA. 2004;101:2375–80.
 3.
Xiao D, Long Y, Wang S, Fang L, Xu D, Wang G, et al. Spatiotemporal distribution of malaria and the association between its epidemic and climate factors in Hainan, China. Malar J. 2010;9:185.
 4.
Parham PE, Michael E. Modeling the effects of weather and climate change on malaria transmission. Environ Health Perspect. 2010;118:620–6.
 5.
Mordecai EA, Paaijmans KP, Johnson LR, Balzer C, BenHorin T, de Moor E, et al. Optimal temperature for malaria transmission is dramatically lower than previously predicted. Ecol Lett. 2013;16:22–30.
 6.
Teklehaimanot HD, Schwartz J, Teklehaimanot A, Lipsitch M. Weatherbased prediction of Plasmodium falciparum malaria in epidemicprone regions of Ethiopia II. Weatherbased prediction systems perform comparably to early detection systems in identifying times for interventions. Malar J. 2004;3:44.
 7.
Schwartz J. The distributed lag between air pollution and daily deaths. Epidemiology. 2000;11:320–6.
 8.
Kilian AH, Langi P, Talisuna A, Kabagambe G. Rainfall pattern, El Nino and malaria in Uganda. Trans R Soc Trop Med Hyg. 1999;93:22–3.
 9.
Abeku TA, De Vlas SJ, Borsboom GJ, Tadege A, Gebreyesus Y, Gebreyohannes H, et al. Effects of meteorological factors on epidemic malaria in Ethiopia: a statistical modelling approach based on theoretical reasoning. Parasitology. 2004;128:585–93.
 10.
Ceccato P, Ghebremeskel T, Jaiteh M, Graves PM, Levy M, Ghebreselassie S, et al. Malaria stratification, climate, and epidemic early warning in Eritrea. Am J Trop Med Hyg. 2007;77:61–8.
 11.
Midekisa A, Senay G, Henebry GM, Semuniguse P, Wimberly MC. Remote sensingbased time series models for malaria early warning in the highlands of Ethiopia. Malar J. 2012;11:165.
 12.
Thomson MC, DoblasReyes FJ, Mason SJ, Hagedorn R, Connor SJ, Phindela T, et al. Malaria early warnings based on seasonal climate forecasts from multimodel ensembles. Nature. 2006;439:576–9.
 13.
Briet OJ, Vounatsou P, Gunawardena DM, Galappaththy GN, Amerasinghe PH. Models for short term malaria prediction in Sri Lanka. Malar J. 2008;7:76.
 14.
Pascual M, Ahumada JA, Chaves LF, Rodo X, Bouma M. Malaria resurgence in the East African highlands: temperature trends revisited. Proc Natl Acad Sci USA. 2006;103:5829–34.
 15.
Gosoniu L, Veta AM, Vounatsou P. Bayesian geostatistical modeling of Malaria Indicator Survey data in Angola. PLoS ONE. 2010;5:e9322.
 16.
Zhou SS, Huang F, Wang JJ, Zhang SS, Su YP, Tang LH. Geographical, meteorological and vectorial factors related to malaria reemergence in Huang–Huai River of central China. Malar J. 2010;9:337.
 17.
Gasparrini A. Modeling exposurelagresponse associations with distributed lag nonlinear models. Stat Med. 2014;33:881–99.
 18.
Heaton MJ, Peng RD. Extending distributed lag models to higher degrees. Biostatistics. 2014;15:398–412.
 19.
Hui FM, Xu B, Chen ZW, Cheng X, Liang L, Huang HB, et al. Spatiotemporal distribution of malaria in Yunnan Province, China. Am J Trop Med Hyg. 2009;81:503–9.
 20.
Ma ZF, Liu J, Yang SQ. Climate change in southwest China during 1961–2010: impacts and adaptation. Adv Clim Change Res. 2013;4:223–9.
 21.
Zhao X, Chen F, Feng Z, Li X, Zhou XH. The temporal lagged association between meteorological factors and malaria in 30 counties in south–west China: a multilevel distributed lag nonlinear analysis. Malar J. 2014;13:57.
 22.
Chinese meteorological data sharing service system. Beijing: Chinses National Meteorological Information Center; 2016. http://data.cma.cn/data/index/0b9164954813c573.html. Accessed 21 Sept 2016.
 23.
Wang L, Wang Y, Jin S, Wu Z, Chin DP, Koplan JP, et al. Emergence and control of infectious diseases in China. Lancet. 2008;372:1598–605.
 24.
Teklehaimanot HD, Lipsitch M, Teklehaimanot A, Schwartz J. Weatherbased prediction of Plasmodium falciparum malaria in epidemicprone regions of Ethiopia. I. Patterns of lagged weather effects reflect biological mechanisms. Malar J. 2004;3:41.
 25.
Craig MH, Kleinschmidt I, Le Sueur D, Sharp BL. Exploring 30 years of malaria case data in KwaZuluNatal, South Africa: part II. The impact of nonclimatic factors. Trop Med Int Health. 2004;9:1258–66.
 26.
Dell AI, Pawar S, Savage VM. Systematic variation in the temperature dependence of physiological and ecological traits. Proc Natl Acad Sci USA. 2011;108:10591–6.
 27.
Hastie T, Tibshirani R. Varyingcoefficient models. J R Stat Soc B. 1996;55:757–96.
 28.
R Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2013. http://www.Rproject.org. Accessed 21 Sept 2016.
 29.
Bates D, Maechler M, Bolker B, Walker S. Lme4: linear mixedeffects models using Eigen and S4. R Package Version. 2013;1:1–5.
 30.
Zhao X, Chen F, Feng Z, Li X, Zhou XH. Characterizing the effect of temperature fluctuation on the incidence of malaria: an epidemiological study in south–west China using the varying coefficient distributed lag nonlinear model. Malar J. 2014;13:192.
 31.
Erhart A, Ngo DT, Phan VK, Ta TT, Van Overmeir C, Speybroeck N, et al. Epidemiology of forest malaria in central Vietnam: a large scale crosssectional survey. Malar J. 2005;4:58.
 32.
Luo Y, Zhang YH, Pei FQ, Liu T, Zeng WL, Xiao JP, et al. Timeseries analysis on the malaria morbidity affected by meteorological factors in Guangdong province. Chin J Prev Med. 2012;46:892–7.
 33.
Tian L, Bi Y, Ho SC, Liu W, Liang S, Goggins WB, et al. Oneyear delayed effect of fog on malaria transmission: a timeseries analysis in the rain forest area of Mengla County, south–west China. Malar J. 2008;7:110.
 34.
le Sueur D, Sharp BL. Temperaturedependent variation in Anopheles merus larval head capsule width and adult wing length: implications for anopheline taxonomy. Med Vet Entomol. 1991;5:55–62.
 35.
Huang F, Zhou S, Zhang S, Zhang H, Li W. Meteorological factorsbased spatiotemporal mapping and predicting malaria in central China. Am J Trop Med Hyg. 2011;85:560–7.
 36.
Huang F, Zhou S, Zhang S, Wang H, Tang L. Temporal correlation analysis between malaria and meteorological factors in Motuo County, Tibet. Malar J. 2011;10:54.
 37.
Zhang Y, Bi P, Hiller JE. Meteorological variables and malaria in a Chinese temperate city: a twentyyear timeseries data analysis. Environ Int. 2010;36:439–45.
 38.
Gelman A. Multilevel (Hierarchical) modeling: what it can and cannot do. Technometrics. 2006;48:432–5.
 39.
Bi Y, Hu W, Yang H, Zhou XN, Yu W, Guo Y, et al. Spatial patterns of malaria reported deaths in Yunnan Province, China. Am J Trop Med Hyg. 2013;88:526–35.
 40.
Kigozi R, Zinszer K, Mpimbaza A, Sserwanga A, Kigozi SP, Kamya M. Assessing temporal associations between environmental factors and malaria morbidity at varying transmission settings in Uganda. Malar J. 2016;15:511.
Authors’ contributions
YYW performed the statistical analysis and drafted the manuscript. ZJQ cleared the data. ZJF provided the original data. NW helped to revise the manuscript. XZ conceived of the project concept. XZ, XSL and HJY gave technical assists. All authors read and approved the final manuscript.
Acknowledgements
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Availability of data and materials
The meteorological data analysed during the current study are available in the Chinese Meteorological Data Sharing Service System, [http://data.cma.cn/data/index/0b9164954813c573.html].
The malaria cases data that support the findings of this study are available from the Chinese Centre for Disease Control and Prevention (CCDC) through the Chinese Information System for Infectious Diseases Control and Prevention (CISIDCP) but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available. Data are however available from the authors upon reasonable request and with permission of the Chinese Centre for Disease Control and Prevention.
The data were analysed by R. R is a free software programming language and a software environment for statistical computing and graphics [28]. Specifically, the addon package lme4 [29] was used for the parameter estimation.
Funding
This study was supported by the Natural Science Foundation of China (Nos. 81502894, 61572332), Sichuan Applying Basic Research Project (No. 2016JY0073), and National Science Fund for Distinguished Young Scholars (81525023). All the funding was used in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Author information
Additional files
12936_2017_1706_MOESM1_ESM.png
12936_2017_1706_MOESM2_ESM.docx
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Received
Accepted
Published
DOI
Keywords
 Malaria
 Rainfall
 Lag
 Nonlinear
 Interaction
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate. Please note that comments may be removed without notice if they are flagged by another user or do not comply with our community guidelines.