<title>Background</title> <p>Qualitative changes in the dynamics of populations, or regime shifts, are common phenomena across all living organisms [##UREF##0##1##,##UREF##1##2##]. Originally defined in fisheries science [##UREF##2##3##], the concept that at some time (termed a "breakpoint") there are disturbances that push a biological system beyond its normal dynamical pattern and can qualitatively change its behavior. Recently, this has become a major concern for vector-borne diseases in the context of global climatic change [##REF##16571662##4##, ####REF##17989780##5##, ##REF##18256410##6####18256410##6##]. Such "breakpoints" derive from ecological analysis that has come to be known as Schmalhausen's law [##UREF##1##2##], which states that systems at the border of their limits of tolerance to one factor become more sensitive to small changes along any other dimension of existence [##UREF##1##2##]. Schmalhausen's law implies that if a system is pushed away from a state of exacerbation, its mean value and variability should decrease. This principle is strongly connected with the idea of resilience [##UREF##3##7##], the robustness of an ecological system before changing to a qualitatively different state, which in principle should be less susceptible to the effects of climatic variability as populations become less vulnerable to infection [##REF##10379454##8##].</p> <p>Malaria in the archipelago of Vanuatu has historically been a major public health problem as shown by the early entomological surveys of Buxton and Hopkins [##UREF##4##9##], followed by the extensive work of Bastien [##REF##2403624##10##], where an increase in the burden of the disease in the early 1980s was reported [##REF##3542257##11##], as well as its possible association to the evolution of quinine resistant parasites [##REF##3310255##12##,##REF##3046769##13##], numerous studies have shown this disease to be a major burden for Vanuatu inhabitants. Although occasionally hyperendemic, like in some areas of sub-Saharan Africa, malaria patterns are very different from this region in several aspects. In Vanuatu, the frequency of fatal cases is greatly diminished [##UREF##5##14##,##REF##9463668##15##], the number of inapparent infections changes seasonally, disease depends on <italic>Plasmodium </italic>species [##REF##9015495##16##], the diversity of parasites is reduced [##REF##10811274##17##], and the genetic make-up of the native populations presents signatures of evolutionary changes driven by malaria. The latter is expressed in an increased frequency of α-thalassaemia associated with mild cases of malaria [##REF##8849722##18##], and an increased frequency of G6PDH enzyme deficiency [##REF##9777715##19##], which is different from sickle cell anaemia, the most common one seen in Africa [##REF##8849722##18##,##REF##9777715##19##].</p> <p>Malaria control efforts also are important to analysis of this time pattern. In 1988, a major control intervention was launched, with a massive distribution of insecticide-treated nets (ITNs), following the abandon of indoor residual spraying for controlling malaria [##UREF##6##20##]. Although focused studies have demonstrated the use of ITNs to be very effective on small islands of this archipelago, as demonstrated by the elimination of the disease in Aneytium [##REF##11075770##21##], another study analysing the effects of this policy at the level of the whole country has not been undertaken. In the present study, the dynamics of malaria before and after the introduction of ITNs into the archipelago are evaluated in an attempt to determine whether there were breakpoints where dynamics shifted transmission patterns, and quantified the effects of climate on these patterns before and after this intervention took effect.</p>
<title>Methods</title> <title>Malaria data and monitored population at risk</title> <p>Monthly records of malaria were obtained from health centers of people who presented with fever or a recent history of fever, and whose standard blood slide analysis indicated infection with either <italic>Plasmodium vivax </italic>or <italic>Plasmodium falciparum</italic>, from January 1983 to December 1999. Malaria cases detected by this passive surveillance were the basis of the analysis. During this period total population increased (Figures ##FIG##0##1##, ##FIG##1##2##). Data on distributed ITNs with permethrin and re-impregnations were available for the same period (Figure ##FIG##2##3A##). Data collection was done under the guidance of the World Health Organization, and controlled by two of the authors (AK, GT) who maintained quality controls on the reporting system and diagnosis reliability during the studied period. All data were obtained from the Malaria and other Vector Borne Diseases Control Unit, Ministry of Health, Port Vila, Vanuatu.</p> <p>This passive case detection system changed in January 1991, as slide examination in small rural health posts was discouraged by the central government of Vanuatu [##REF##9777715##19##]. This policy change reduced the number of people being monitored, however, it remained representative of the whole population [##REF##9777715##19##]. To account for the possible effects of this policy change, changes in the rate of slide examination were measured before and after the breakpoint obtained for the rate of slide examination, and assumed it to be linearly correlated to changes in the population monitored (Figure ##FIG##0##1F##). That is, the population at risk (corresponding to the population in districts where malaria was present) was multiplied by the fraction obtained by dividing the average rate of examination before and after the breakpoint to evaluate this possible source of error. A 50% reduction in average rate of slide examination was found during 1990–1991 (Figure ##FIG##0##1F##), as described in [##REF##9777715##19##].</p> <title>Environmental data</title> <p>Weather data included Sea Surface Temperature (SST) indexes: 1+2, 3, 3.4 and 4 (also known as the Niño 1+2, 3, 3.4 and 4 [##UREF##7##22##]; Additional file ##SUPPL##0##1##), and precipitation and temperature data from the climate database for political areas [##UREF##8##23##,##UREF##9##24##]. These data were used as predictors in models to assess changes in the magnitude of forcing by climatic variables in the dynamics of malaria incidence.</p> <title>Statistical analysis</title> <title>Breakpoints and regime shifts</title> <p>Tests of structural changes in time series can be undertaken using at least three different strategies: F tests that compare the null hypothesis of no regime shift to the presence of a regime shift, generalized fluctuation tests that do not assume any particular pattern of deviation from the absence of regime shifts [##UREF##10##25##,##UREF##11##26##] and adaptive filtering of signals [##UREF##12##27##]. A total of three approaches were used in the present study to assess the robustness of the findings. The F statistic is obtained by comparing the residuals (<italic>i</italic>) of a segmented regression at time i with the residuals  from an unsegmented regression using the following expression:</p> <p></p> <p>Where n is the time series length and k the number of parameters. The null hypothesis is rejected when the supremum of the statistic is larger than the value of a distribution SupF derived by Hansen [##UREF##13##28##,##UREF##14##29##]. This approach has been generalized for <italic>l </italic>breaks, with arbitrary but fixed <italic>l </italic>[##UREF##15##30##,##UREF##16##31##]; where the number of breaks can be selected using conventional tools for model selection like the Akaike Information Criterion (AIC) [##UREF##17##32##].</p> <p>The other two approaches, the generalized fluctuation test and the adaptive filtering, include formal significance tests, yet reveal regime shifts graphically instead of assuming specific types of departure in advance. For the generalized fluctuation test a parametric model is fitted to the data and an empirical process (EFP) is derived that captures the fluctuation either in residuals or parameter estimates [##UREF##10##25##,##UREF##11##26##]. Under the null hypothesis the fluctuations are governed by central limit theorems while under the alternative (regime shifts) the fluctuation is increased [##UREF##11##26##]. In the present analysis, the ordinary least squares (OLS) based CUSUM tests introduced in [##UREF##18##33##] was used. This test is based in cumulative sums of residuals from a linear regression:</p> <p></p> <p>where a regime shift is evidenced by a single peak around the breakpoint, provided that the limiting process for  is the standard Brownian bridge <italic>W</italic><sub>0</sub>(<italic>t</italic>) = <italic>W</italic>(<italic>t</italic>) - <italic>tW</italic>(l), where <italic>W</italic>(.) denotes Brownian motion. Significance for the CUSUM was tested using the derivations presented in [##UREF##10##25##,##UREF##11##26##]. For equation (1) and (2) the residuals  came from a linear seasonal autoregressive [##UREF##17##32##] model fitted by OLS:</p> <p></p> <p>The third approach is totally non-parametric, and is based on recovering a signal and its breaks. The Kolmogorov-Zurbenko adaptive filter (KZAF) [##UREF##12##27##] is based on filtering the time series <italic>y </italic>using:</p> <p></p> <p>Where</p> <p></p> <p>And q is half-length of a k iterative moving average (<italic>x</italic><sub><italic>t</italic></sub>) applied to the original time series <italic>y</italic><sub><italic>t</italic></sub>. The term <italic>f</italic>(<italic>D</italic>(<italic>t</italic>)) is defined by:</p> <p></p> <p>And D(t) is the absolute difference defined by:</p> <p></p> <p>And D'(t) as:</p> <p></p> <p>Once <italic>z</italic><sub><italic>t </italic></sub>is obtained quantitative estimates of discontinuity can be based on an analysis of the sample variances of <italic>z</italic><sub><italic>t</italic></sub>, defined by:</p> <p></p> <p>When there are no breaks, maxima in the estimated variance of (9) are approximately independent and exponentially distributed with a expected number of peaks of about n/(2qk<sup>0.5</sup>), allowing to consider a breakpoint when the  value exceeds the 95% upper tail of the exponential distribution with such parameter.</p> <p>Regime shift analyses were carried out on: (i) the monthly rate of slide examination (No. Slides examined*1,000/Total population at risk); (ii) the monthly rate of the two malaria parasites (No. slides examined*1,000/Monitored population at risk) and (iii) weather variables (rainfall and temperature).</p> <title>Threshold for ITN coverage</title> <p>Time series for total number of bed nets distributed per month were accumulated and divided by the total population at risk estimated from the annual population data. It was assumed that the annual data corresponded to December, and interpolated the rest of the months using a smoothing splines regression as explained in [##REF##17211661##34##]. The probability density [##UREF##17##32##] of the percentage of people locally covered with the distributed ITNs was also studied.</p> <title>Seasonality</title> <p>The seasonality of vivax and falciparum malaria rates (cases/population size) were assessed by using box diagrams before and after the regime shift [##UREF##17##32##].</p> <title>Non-stationary patterns of association</title> <p>The wavelet transform can be used to study the patterns of association between two nonstationary time series [##UREF##19##35##,##REF##16903778##36##]. Specifically, the wavelet coherency analysis can determine whether the presence of a particular frequency at a given time in the disease corresponds to the presence of that same frequency at the same time in a covariate (e.g., rainfall and temperature). The cross-wavelet phase analysis can determine the time lag separating these two series as well.</p> <title>Changes in the effects of climate on the dynamics</title> <p>Once breakpoints for the regime shift were identified in the falciparum and vivax malaria rate series, the splitted series around the breakpoints were studied using seasonal auto-regressive (SAR) models [##UREF##17##32##]. The procedure for model building was similar to the one described in [##REF##16903778##36##]: (i) a null model was fitted to the rate of the falciparum and vivax malaria (ii) temperature and rainfall were filtered with the coefficients of the null model, and (iii) cross-correlation functions were computed using the residuals of the null model and those of the filtered climatic variables.</p> <p>The full model for <italic>P. falciparum </italic>considered precipitation (P) with lags of 2 and 29 months, and temperature (T) with lags of 3 and 12 months, as follows:</p> <p></p> <p>For <italic>P. vivax </italic>the full model considered precipitation (P) a lag 9 months, and temperature (T) with a lag 10 months, as follows:</p> <p></p> <p>In both cases, the error was assumed as independent and normally distributed: <italic>ε</italic><sub><italic>t</italic></sub>~<italic>N</italic>(0,). After the initial fitting, models were simplified using a process of backward elimination[##REF##16903778##36##]: (i) taking out one predictor at a time, (ii) finding the minimum AIC for models with similar complexity, i.e., number of parameters, (iii) comparing the likelihood of the best model (minimum AIC) for each level of complexity with the full model, and simplifying while differences were not statistically significant. For the analyses the climatic covariates were demeaned in order to not affect the intercept value [##UREF##17##32##].</p>
<title>Results</title> <p>Temporal patterns of malaria in Vanuatu present a clear shift in the incidence rate by the end of 1993 and beginning of 1994, for both parasite species (Figures ##FIG##1##2A## and ##FIG##1##2D##).</p> <p>Breakpoints were confirmed by all three different methods (Additional files ##SUPPL##1##2## and ##SUPPL##2##3##). For the incidence rate in both malaria species, breakpoints were statistically significant according to the F statistic and the variance of the KZAF. Even though the EFP estimates were not significant, peaks were detectable in both cases in January 1992 (Additional file ##SUPPL##2##3##). During that same time period no significant changes were found for climatic time series (Additional file ##SUPPL##3##4##). By the time changes were detected, bed net coverage (Figure ##FIG##2##3B##) was as low as 6% (EFP estimate) or slightly above 20% of the population at risk (KZAF). At a more local scale, villages where bed nets were distributed mostly had ~80% of the population covered (Figure ##FIG##2##3C##).</p> <p><italic>Plasmodium falciparum </italic>seasonality was qualitatively very similar before and after the breakpoint (Figures ##FIG##1##2B## and ##FIG##1##2C##), showing maximum incidence during the first quarter of the year (January-March), and minimum incidence during the third quarter of the year (July-September). For <italic>P. vivax </italic>(Figures ##FIG##1##2E## and ##FIG##1##2F##) a similar change was observed, although the patterns were not so clear as for <italic>P. falciparum</italic>, due to greater seasonal variability. With the exception of a brief period during 1992–1996 where cases due to both parasites were synchronous (i.e., with peaks at the same time), the dynamics of the infections were mainly asynchronous and not coherent (i.e., not associated in the frequency domain) at the seasonal scale. However, both diseases were significantly cross-coherent at an interannual scale, with the dynamics of <italic>P. falciparum </italic>cases being mostly synchronous with that of <italic>P. vivax </italic>(Additional file ##SUPPL##4##5##).</p> <p>Regarding the effects of climate during the studied period, the cross-coherence wavelet analysis showed malaria to be correlated with temperature at the seasonal scale; both <italic>P. falciparum </italic>and <italic>P. vivax </italic>incidence rates were led by temperature (Figure ##FIG##3##4##). A similar pattern was seen between the two parasites and rainfall at the seasonal scale, despite the presence of some gaps. A significant coherence with rainfall at interannual scales was also found. For <italic>P. vivax</italic>, coherence was statistically significant for periods between two and four years, during 1992–1996. No evidence that El Niño indices were leading the dynamics of the disease was identified.</p> <p>Finally, Table ##TAB##0##1## presents the parameter estimates for rate models of <italic>P. falciparum </italic>and <italic>P. vivax</italic>, including exogenous forcing by temperature, before and after the breakpoint. Model selection by backward elimination showed that rainfall was not a significant covariate (detailed values in Additional file ##SUPPL##5##6##). Following the qualitative change in the dynamics,<italic>P. falciparum </italic>had a proportional (~66%) and absolute (7.7 cases/1000 population) decline in incidence that was greater than that for <italic>P. vivax </italic>(~52% and 2.6 cases/1000 population, respectively). The importance of temperature in driving the dynamics also declined after the breakpoint for both species, between 31% and 49% for <italic>P. falciparum </italic>and 80% for <italic>P. vivax </italic>(the coefficient after the breakpoint became statistically non-significant). The importance of temperature in driving the dynamics also declined after the breakpoint for both species, between 31% and 49% for <italic>P. falciparum </italic>and, although not statistically significant, 80% for <italic>P. vivax</italic>. This suggests that the average effect of 1°C increase in temperature will increase incidence in a reduced amount when compared with its effect before the breakpoint. For example, preceding the shift each degree Celsius above the three-month lagged mean temperature value used to increase the rate by 1.43 cases/1000 people at risk, for <italic>Plasmodium falciparum</italic>. In contrast, after the breakpoint this change only increased the rate by half of its previous magnitude, i.e., 0.72/1000 people at risk (Table ##TAB##0##1##). A similar phenomenon was also seen for the variability that was not explained by the models, which also was reduced by 54% and 33% in <italic>P. vivax </italic>and <italic>P. falciparum</italic>, respectively, as shown by the decrease in the error variance of the models after the breakpoint (Table ##TAB##0##1##).</p>
<title>Discussion</title> <p>Following a disturbance, biological systems can either return to their normal state of variability or can move far away from such a state [##UREF##0##1##,##UREF##1##2##,##REF##11595939##37##,##UREF##20##38##]. Transients, i.e., the anomalous behavior between regimes or basins [##UREF##21##39##,##REF##16701224##40##], can obscure the qualitative changes of a system, because jumps from one state to another are not always instantaneous, complicating our ability to identify regime shifts [##UREF##21##39##,##REF##17479460##41##]. This is likely one of the main differences between the dynamics of <italic>P. falciparum </italic>and <italic>P. vivax</italic>, since a consistent estimate for the breakpoint was easy to find for the former, while the estimates for the latter differed significantly. This was especially true for KZAF, which identified a later breakpoint. Assumptions underlying the employed techniques [##UREF##10##25##, ####UREF##11##26##, ##UREF##12##27##, ##UREF##13##28##, ##UREF##14##29##, ##UREF##15##30##, ##UREF##16##31##, ##UREF##17##32##, ##UREF##18##33####18##33##] might favour the estimate from KZAF, since the F statistic is quite sensitive to the stationarity (i.e., constant mean) of the time series, while the CUSUM EFP may be too sensitive given the quality of the data examined, identifying the change of policy in slide examination. By contrast, the KZAF is an adaptive technique that allows control of the time scale at which changes may be occurring [##UREF##12##27##]. This is a very useful characteristic for addressing one of the major recurrent problems in the study of ecological systems, i.e. finding the appropriate temporal scale of a natural phenomenon [##UREF##22##42##]. In this study, the adaptive ability of KZAF allowed for breaks to be distinguished from natural cycles associated with exogenous factors (i.e., climate). The fact that the basin (or regime) shift in the time series can be attributed to the effects of bed net use appears robust. During the study period no other major changes in control strategies, landscape cover, medication or drug resistance were reported [##REF##2403624##10##,##REF##3542257##11##,##REF##9777715##19##] after controlling for the policy change in data collection [##REF##9777715##19##].</p> <p>The analysis identified a major difference between <italic>P. falciparum </italic>and <italic>P. vivax</italic>, namely the earlier breakpoint for <italic>P. falciparum</italic>. This pattern would not be expected under conditions of cross or heterologous immunity [##REF##4201093##43##], and its evaluation with cross-infection studies is limited because quality data that are necessary to make such inferences [##REF##6969036##44##] are lacking [##REF##9777715##19##,##REF##11075770##21##]. However, this pattern should be studied further, because it might reflect the dynamics of immunity in the population, where a generalized density-dependent immunity may be triggered by the within-host density of each parasite species [##UREF##23##45##]. Alternatively, if <italic>P. falciparum </italic>was the first species to be cleared, as shown in the classical co-infection neuro-syphilis malariotherapy experiments of Boyd and Kitchen [##REF##10657296##46##], temporal patterns can only be appreciated when studying the dynamics of the within-host parasitic infection [##REF##17690394##47##]. In addition, the pattern simply could arise by the ability of <italic>P. vivax </italic>to relapse [##REF##9777715##19##,##REF##11075770##21##], possibly in conjunction with the immunity dynamics described above.</p> <p>Although regime shifts tend to be thought of in terms of increased variability as the best diagnostic condition [##REF##16958897##48##], they can occur in the opposite direction, with systems becoming more stable. For both <italic>P. falciparum </italic>and <italic>P. vivax </italic>not only did the mean value of incidence decrease, but also the variance of the models decreased, which is a more robust measure of stability [##UREF##20##38##] than just looking at mean values [##UREF##0##1##] in dynamical systems. The patterns seen for the two species differed: falciparum malaria declined more abruptly, in total and relative terms, than in vivax malaria. Perhaps there are differences in the life history strategies of the parasites under different scenarios for transmission, with the most virulent parasite (<italic>P. falciparum</italic>) being more successful in environments with high transmission rates and the least virulent (<italic>P. vivax</italic>) being less sensitive to the intensity of transmission.</p> <p>A surprising result was that the breakpoint occurred after just 20% of the population was covered with bed nets, which is half that predicted for <italic>Anopheles gambiae </italic>transmission by Killeen <italic>et al </italic>[##REF##17608562##49##]. Perhaps <italic>Anopheles farauti</italic>, the main vector in Vanuatu [##UREF##4##9##,##REF##7668922##50##] is less efficient. Regardless, the fact that such ITN coverage could explain the decrease has a robust theoretical explanation as presented in the groundbreaking work of Becker and Dietz [##REF##7795319##51##], later confirmed using field data as the 80/20 rule for several infectious diseases [##REF##8990210##52##,##REF##17686448##53##] where the control, which targets 20% of the population, could benefit the other 80% of people.</p> <p>Interestingly, this rule has been derived by looking at local populations, but the pattern seen in Vanuatu is more likely to arise from the subdivided nature of the population in villages, or patches if seen from the perspective of metapopulations [##UREF##24##54##]. The coverage per patch was high enough (80% with a very low dispersion around this value) to guarantee the local interruption of transmission according to mechanistic models of bed net action in settings with a higher entomological inoculation rate [##REF##17608562##49##,##REF##17254339##55##] than that observed in Vanuatu [##REF##9015495##16##,##REF##7668922##50##].</p> <p>As a control strategy, ITNs outperform similar strategies aimed at reducing vectorial capacity, such as the indoor residual spraying, mainly because of its cost-effectiveness, as well as for its ease of implementation and distribution [##UREF##25##56##,##REF##15106149##57##]. Several studies have shown that bed nets reduce total infant mortality in endemic areas [##REF##9623929##58##,##REF##11157527##59##], are a sustainable option for control in terms of the reduction of relative risk of malaria death in the medium- to long-term time scales [##REF##15173148##60##], and are successful across several cultural settings [##REF##15106149##57##,##REF##15462999##61##, ####REF##16172493##62##, ##REF##17313511##63##, ##REF##17988171##64####17988171##64##]. The advantages of bed nets also go beyond the immediate effects, since so far there is no evidence for selection of insecticide-resistant mosquitoes [##REF##11579872##65##], and they are protective even in areas where mosquito resistance to the insecticides used for bed net impregnation has been reported [##REF##16472385##66##]. This result also has been theoretically reinforced by models that consider the use of bed nets in conjunction with other control strategies, such as zooprophylaxis [##REF##17631372##67##], provided that both measures in conjunction are likely to counteract any selective pressure for the development of insecticide resistance, since mosquito fitness would not be under a selective pressure, and may even be under selection for feeding preferences in non-human hosts [##UREF##26##68##,##REF##12460393##69##]. However, urban settings pose a major challenge since effective zooprophylaxis might be diminished because of higher human densities. Behavioral changes in mosquitoes and decreased bed net effectiveness have been documented in urban areas [##REF##17880679##70##]. From a wider perspective, bed nets are also a more ecologically-sound strategy since they reduce impacts on natural enemies of vectors via positive feedbacks loops that can be generated by large scale insecticide spraying [##UREF##26##68##,##REF##17373942##71##,##UREF##27##72##]. A large body of literature supports that idea that in relatively undisturbed environments mosquito abundance is regulated by interactions with other animals, e.g., tadpoles, fish and other insects [e.g., [##REF##17373942##71##, ####UREF##27##72##, ##UREF##28##73##, ##REF##16978142##74##, ##REF##15707287##75##, ##REF##17883002##76####17883002##76##]], however such natural control is diminished by anthropogenic disturbances of food webs.</p>
<title>Conclusion</title> <p>The success of the Vanuatu malaria control programme also stems from the strategy of bed net distribution, where large fractions of the population were locally covered at the village level, ensuring the reduction in transmission, even leading to local elimination in some islands [##REF##11075770##21##]. As stressed by Killeen <italic>et a</italic>l [##REF##17608562##49##] and Ilboudo-Sanogo <italic>et al </italic>[##REF##11579872##65##], an efficient bed net programme needs to cover a large proportion of the population in order to ensure that both sources (e.g., asymptomatic people) and sinks (e.g., pregnant women and young children) of infection are effectively covered. The erroneous targeting of transmission groups for control can exacerbate the conditions for transmission [##REF##17524806##77##]. Additionally, as suggested by Mathanga <italic>et al </italic>[##REF##16930258##78##], for ethical and humanitarian reasons the goal should be to cover as much of the population present in the endemic setting as possible, retaining traditional practices (e.g., voluntary work) for the exchange of goods when mainstream means of commercialization are not enough to achieve such a goal. In Vanuatu, special care was taken to address these factors by implementing a strategy where children under five years of age, their mothers and pregnant women received free nets. Cost was half price for school children and other adults were charged the full price, ensuring an equitable coverage of the population [##REF##11075770##21##] and an equitable distribution of this valued resource.</p> <p>A factor that deserves further study is the role that concomitant knowledge transfer associated to the distribution of bed nets have on the awareness of the population about the risk leading to malaria transmission. Unlike insecticide residual spraying whose effectiveness depends mostly on being applied correctly, the effective use of bed nets requires knowledge for its proper use. In Vanuatu, parents' awareness was likely to play a role in diminishing incidence among young children (&lt;5 years), because of the free distribution to this age group and training to parents about the benefits of using the nets [##REF##11075770##21##]. But, the positive effects of knowledge transfer are likely to be more comprehensive. For example, Mathanga <italic>et al </italic>[##REF##16930258##78##] showed that even though children didn't regularly use bed nets, those in communities where malaria transmission plummeted after the introduction of widespread bed net use were aware of the benefits. Similar knowledge transfers are known to be present among some Native American tribes whose mythology has associated malaria risk with the blossoming of water-retaining flowers where vector larvae develop [##UREF##29##79##]. Changes in collective behavior in villages that were stricken by malaria have been seen before community-based educational campaigns were implemented [##REF##16760516##80##, ####REF##17556618##81##, ##REF##16457724##82####16457724##82##] and more generally, traditional knowledge has been shown to be a robust strategy to handle issues of pest management by native populations in Meso-America [##UREF##30##83##].</p> <p>The association between climatic forces and malaria dynamics in Vanuatu presents features that make it unusual when compared to other settings, where the climate and ecological dynamics have been studied [e.g., [##REF##16903778##36##], reviewed in [##REF##14561270##84##]]. None of the ENSO indices led the dynamics of malaria, yet clear signals of association at interannual time scales were found with local climatic variables. This may be a result of the relationship of ENSO with the local climate in the area [##REF##16292302##85##] which influences rainfall during a season, October to January [##UREF##31##86##, ####UREF##32##87##, ##UREF##33##88####33##88##], that probably is not relevant for the biology of mosquitoes in regards to transmission. The unusual pattern is less likely because of a demographic effect of small insular population size as suggested in [##REF##9924003##89##]. Mechanisms for the action of rainfall across a wide range of landscapes have been very well described, it increases the rate of a disease when new mosquito habitats are created by increased precipitation [##UREF##34##90##], and the additional weakening of inter-specific interactions regulating mosquito populations [##REF##9388155##91##]. However, ecological studies of vectors are needed to understand their local population dynamics in Vanuatu. Similarly it may be understood why hotter temperatures can increase the transmission of vector-borne diseases, because of known effects of temperature on the rate of insect and parasite development [##REF##16292302##85##,##REF##16595623##92##]. However, increased resilience to the effects of climate in an infectious disease as a result of control measures, in our knowledge, has not been reported before. The fact that such a measure also decreases the incidence of malaria under changing climatic conditions is a remarkable fact strengthening the usefulness of this strategy.</p> <p>Finally, a precautionary note on bed nets should be posed. Even though they are a very robust strategy to control malaria from evolutionary, ecological, conservation and cost-effectiveness perspectives [##UREF##25##56##,##REF##15106149##57##,##REF##11579872##65##,##REF##16930258##78##], the use of bed nets should not be viewed as a exhaustive solution if the long-term goal of population health is to be pursued. As shown in [##REF##10499079##93##] a fraction of the death toll that was avoided by controlling malaria through the use of insecticide treated curtains in areas of Burkina Faso was shifted to meningococcal meningitis. Evidence also suggests that in urban settings, for a series of factors that go from the absence of alternative hosts to behavioural shifts in humans, insecticide treated nets are not going to be a sufficient strategy to keep malaria under control [##REF##17880679##70##]. To achieve this goal, a wide research agenda, fully integrated with policies beyond disease control is a path that needs to be taken [##REF##17211661##34##,##REF##17096840##94##, ####REF##14522266##95##, ##REF##3294623##96##, ##UREF##35##97####35##97##], where ultimate goals are aimed at pushing out the stressful contextual conditions that make human populations vulnerable to infectious diseases [##UREF##1##2##], especially malaria.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>Malaria is an important public-health problem in the archipelago of Vanuatu and climate has been hypothesized as important influence on transmission risk. Beginning in 1988, a major intervention using insecticide-treated bed nets (ITNs) was implemented in the country in an attempt to reduce <italic>Plasmodium </italic>transmission. To date, no study has addressed the impact of ITN intervention in Vanuatu, how it may have modified the burden of disease, and whether there were any changes in malaria incidence that might be related to climatic drivers.</p> <title>Methods and findings</title> <p>Monthly time series (January 1983 through December 1999) of confirmed <italic>Plasmodium falciparum </italic>and <italic>Plasmodium vivax </italic>infections in the archipelago were analysed. During this 17 year period, malaria dynamics underwent a major regime shift around May 1991, following the introduction of bed nets as a control strategy in the country. By February of 1994 disease incidence from both parasites was reduced by at least 50%, when at most 20% of the population at risk was covered by ITNs. Seasonal cycles, as expected, were strongly correlated with temperature patterns, while inter-annual cycles were associated with changes in precipitation. Following the bed net intervention, the influence of environmental drivers of malaria dynamics was reduced by 30–80% for climatic forces, and 33–54% for other factors. A time lag of about five months was observed for the qualitative change ("regime shift") between the two parasites, the change occurring first for <italic>P. falciparum</italic>. The latter might be explained by interspecific interactions between the two parasites within the human hosts and their distinct biology, since <italic>P. vivax </italic>can relapse after a primary infection.</p> <title>Conclusion</title> <p>The Vanuatu ITN programme represents an excellent example of implementing an infectious disease control programme. The distribution was undertaken to cover a large, local proportion (~80%) of people in villages where malaria was present. The successful coverage was possible because of the strategy for distribution of ITNs by prioritizing the free distribution to groups with restricted means for their acquisition, making the access to this resource equitable across the population. These results emphasize the need to implement infectious disease control programmes focusing on the most vulnerable populations.</p>
<title>Competing interests</title> <p>The authors declare that they have no competing interests.</p> <title>Authors' contributions</title> <p>LFC conceived the study and carried out the analyses. AK and GT collected the data and performed the field studies. AK and GT provided input on methods. AK, MP, MLW provided input on results interpretation. LFC drafted the manuscript to which AK, GT, MP, and MLW made contributions.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>We thank the people of Vanuatu for their generosity and help throughout years of study. We are also grateful to Dr. Kazuyo Ichimori at WHO for coordinating and encouraging the efforts for data collection, Mr. Takayuki Kurita for compiling the database, Dr. Patrick Bastien at Université de Montpellier, France, who provided us with reprints of his papers, Dr. Christian Lengeler at the Swiss Tropical Institute, for his reprints, advice and insights on the role of bed nets for controlling malaria, Dr. Thomas Smith for his hospitality during a visit by LFC to the Swiss Tropical Institute, Basel, Switzerland, Dr. Benjamin Cash at COLA and an anonymous referee for comments on climatic events in the area of Vanuatu, Dr. Edward L. Ionides for his comments and suggestions on statistical techniques, and the Tropical Biology study group at University of Michigan, especially Mr. Brian Sedio for his comments on the manuscript. LFC was supported by Fundación Polar (Caracas, Venezuela) and The University of Michigan through: The International Institute, The Rackham Graduate School Graduate Student Research Grant, a Summer Fellowship and Block Grant from the Department of Ecology and Evolutionary Biology. AK was supported by an Institutional Grant of the Swedish Foundation for International Cooperation in Research and Higher Education (STINT) and a Project Grant of the Swedish Research Council (VR, 2005–6836). GT was supported by the Government of the Republic of Vanuatu. MP was supported by The Graham Institute for Environmental Sustainability at University of Michigan, the National Oceanic and Atmospheric Administration (Oceans and Health Program NA 040 AR 460019) and NSF-NIH (Ecology of Infectious Diseases Grant EF 0430 120). MLW was supported by the Global Health Program, School of Public Health, The University of Michigan. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p>Time Series: <bold>A </bold><italic>Plasmodium falciparum </italic>malaria cases, <bold>B </bold><italic>P. vivax </italic>malaria cases, <bold>C </bold>Temperature (°C), <bold>D </bold>Precipitation (mm), <bold>E </bold>Population at risk (solid), <bold>F </bold>Monthly slide examination rate (slides*1000/population at risk), the dashed line corresponds to the breakpoint, August 1990, estimated using the F statistic, and the solid lines at the bottom of the graph to the confidence intervals (February 1990, May 1992) the thick-black solid line is the Kolmogorov-Zurbenko adaptive filter implemented with a half window size, q, of 36 months, the breakpoint is December 1991, the blue line corresponds to the breakpoint obtained using the CUSUM (march 1990). The mean rate (± S.D.) of slide examination before the breakpoint (August 1990) was (51.22 ± 11.40) being reduced to (25.92 ± 11.56) after it. Statistical tests of significance can be seen in Additional file ##SUPPL##2##3##.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p>Regime Shift for falciparum and vivax malaria: <bold>A </bold>falciparum malaria rate, the dashed line corresponds to the breakpoint, January 1992, estimated using the F statistic, and the solid lines at the bottom of the graph to the confidence intervals (June 1989, June 1994), the thick-black solid line is the Kolmogorov-Zurbenko adaptive filter implemented with a half window size, q, of 36 months, the breakpoint corresponds to August 1993, the blue line corresponds to the breakpoint obtained using the CUSUM (January 1992).<bold>B </bold>&amp;<bold>C </bold>seasonal falciparum malaria rate before and after breakpoint (January 1992) <bold>D </bold>vivax malaria rate, the dashed line corresponds to the breakpoint, May 1991, estimated using the F statistic, and the solid lines at the bottom of the graph to the confidence intervals (June 1989, November 1992), the black solid line is the Kolmogorov-Zurbenko adaptive filter implemented with a half window size, q, of 36 months, the breakpoint corresponds to February 1994, the blue line corresponds to the breakpoint obtained using the CUSUM (January 1992) <bold>E </bold>&amp;<bold>F </bold>seasonal vivax malaria rate before and after breakpoint. For the F statistics the 30% percent of the data belonging to the extremes (15% each) was left out. For the Kolmogorov-Zurbenko adaptive filter q was set to 36, in order to avoid the misidentification of cycles shorter than 6 years.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p>Bed nets <bold>A </bold>Monthly number of distributed bed nets(black line) and number Re-impregnated bed nets (green line) <bold>B </bold>Probability density of the percentage of people locally covered with bed nets between 1988 and 1997, bandwidth of 0.027 <bold>C </bold>Percent (%) of population covered by bed nets for the lower and upper time limit for the breakpoints, the green-blue line corresponds to January 1992 (<italic>Plasmodium falciparum </italic>and <italic>P. vivax</italic>), the green line to September 1992 (<italic>P. falciparum</italic>) and the blue line to December 1992 (<italic>P. vivax</italic>).</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p>Cross-wavelet coherency and phase of <italic>Plasmodium falciparum </italic>malaria rate with <bold>A </bold>temperature and <bold>B </bold>rainfall and of <italic>P. vivax </italic>malaria rate with <bold>C </bold>temperature and <bold>D </bold>rainfall. The coherency scale is from zero (blue) to one (red). Red regions in the upper part of the plots indicate frequencies and times for which the two series share variability. The cone of influence (within which results are not influenced by the edges of the data) and the significant (<italic>p </italic>&lt; 0.05) coherent time-frequency regions are indicated by solid lines. The colors in the phase plots correspond to different lags between the variability in the two series for a given time and frequency, measured in angles from -PI to PI. A value of PI corresponds to a lag of 17 mo. The procedures and software are those described in [31,32]. A smoothing window of 15 mo (2w + 1 = 31) was used to compute the cross-wavelet coherence.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Parameter values and % reduction for <italic>Plasmodium falciparum </italic>and <italic>Plasmodium vivax </italic>rate before and after the breakpoint obtained by using the variance of the Kolmogorov Zurbenko adaptive filter.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center"><bold>Species</bold></td><td align="center"><bold>Parameter</bold></td><td align="center"><bold>Before</bold></td><td align="center"><bold>After</bold></td><td align="center"><bold>% Reduction</bold></td><td align="center"><bold>P</bold></td></tr></thead><tbody><tr><td align="center"><bold><italic>P.falciparum</italic></bold></td><td align="center"></td><td align="center">11.56 ± 1.02</td><td align="center">3.90 ± 1.15</td><td align="center">66.26</td><td align="center">B/A</td></tr><tr><td/><td align="center"></td><td align="center">1.43 ± 0.42</td><td align="center">0.72 ± 0.35</td><td align="center">48.59</td><td align="center">B/A</td></tr><tr><td/><td align="center"></td><td align="center">1.36 ± 0.44</td><td align="center">0.94 ± 0.33</td><td align="center">30.88</td><td align="center">B/A</td></tr><tr><td/><td align="center"></td><td align="center">5.76</td><td align="center">3.84</td><td align="center">33.31</td><td align="center">-</td></tr><tr><td colspan="6"><hr/></td></tr><tr><td align="center"><bold><italic>P.vivax</italic></bold></td><td align="center"></td><td align="center">4.83 ± 0.62</td><td align="center">2.33 ± 0.43</td><td align="center">51.76</td><td align="center">B/A</td></tr><tr><td/><td align="center"></td><td align="center">0.55 ± 0.18</td><td align="center">0.11 ± 0.16</td><td align="center">80.00</td><td align="center">B</td></tr><tr><td/><td align="center"></td><td align="center">1.17</td><td align="center">0.53</td><td align="center">54.31</td><td align="center">-</td></tr></tbody></table></table-wrap>
<inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M1" name="1475-2875-7-100-i1" overflow="scroll"><mml:semantics><mml:mover accent="true"><mml:mi>ε</mml:mi><mml:mo>^</mml:mo></mml:mover></mml:semantics></mml:math></inline-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M2" name="1475-2875-7-100-i1" overflow="scroll"><mml:semantics><mml:mover accent="true"><mml:mi>ε</mml:mi><mml:mo>^</mml:mo></mml:mover></mml:semantics></mml:math></inline-formula> <disp-formula id="bmcM1"><label>(1)</label><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M3" name="1475-2875-7-100-i2" overflow="scroll"><mml:semantics><mml:mrow><mml:msub><mml:mi>F</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msup><mml:mover accent="true"><mml:mi>ε</mml:mi><mml:mo>^</mml:mo></mml:mover><mml:mi>T</mml:mi></mml:msup><mml:mover accent="true"><mml:mi>ε</mml:mi><mml:mo>^</mml:mo></mml:mover><mml:mo>−</mml:mo><mml:mover accent="true"><mml:mi>ε</mml:mi><mml:mo>^</mml:mo></mml:mover><mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mi>T</mml:mi></mml:msup><mml:mover accent="true"><mml:mi>ε</mml:mi><mml:mo>^</mml:mo></mml:mover><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mrow><mml:mover accent="true"><mml:mi>ε</mml:mi><mml:mo>^</mml:mo></mml:mover><mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mi>T</mml:mi></mml:msup><mml:mover accent="true"><mml:mi>ε</mml:mi><mml:mo>^</mml:mo></mml:mover><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>/</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mi>n</mml:mi><mml:mo>−</mml:mo><mml:mn>2</mml:mn><mml:mi>k</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mfrac></mml:mrow></mml:semantics></mml:math></disp-formula> <disp-formula id="bmcM2"><label>(2)</label><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M4" name="1475-2875-7-100-i3" overflow="scroll"><mml:semantics><mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:mrow><mml:msubsup><mml:mi>W</mml:mi><mml:mi>n</mml:mi><mml:mn>0</mml:mn></mml:msubsup><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:mover accent="true"><mml:mi>σ</mml:mi><mml:mo>^</mml:mo></mml:mover><mml:msqrt><mml:mi>n</mml:mi></mml:msqrt></mml:mrow></mml:mfrac><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mrow><mml:mo>⌊</mml:mo><mml:mrow><mml:mi>n</mml:mi><mml:mi>t</mml:mi></mml:mrow><mml:mo>⌋</mml:mo></mml:mrow></mml:mrow></mml:munderover><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>ε</mml:mi><mml:mo>^</mml:mo></mml:mover><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mstyle></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mn>0</mml:mn><mml:mo>≤</mml:mo><mml:mtext>t</mml:mtext><mml:mo>≤</mml:mo><mml:mn>1</mml:mn><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:semantics></mml:math></disp-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M5" name="1475-2875-7-100-i4" overflow="scroll"><mml:semantics><mml:mrow><mml:msubsup><mml:mi>W</mml:mi><mml:mi>n</mml:mi><mml:mn>0</mml:mn></mml:msubsup><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:semantics></mml:math></inline-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M6" name="1475-2875-7-100-i1" overflow="scroll"><mml:semantics><mml:mover accent="true"><mml:mi>ε</mml:mi><mml:mo>^</mml:mo></mml:mover></mml:semantics></mml:math></inline-formula> <disp-formula id="bmcM3"><label>(3)</label><italic>y</italic><sub><italic>t </italic></sub>= <italic>μ </italic>+ <italic>φ</italic><sub>1</sub><italic>y</italic><sub><italic>t</italic>-1 </sub>+ <italic>φ</italic><sub>12</sub><italic>y</italic><sub><italic>t</italic>-12 </sub>+ <italic>ε</italic><sub><italic>t</italic></sub></disp-formula> <disp-formula id="bmcM4"><label>(4)</label><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M7" name="1475-2875-7-100-i5" overflow="scroll"><mml:semantics><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mi>t</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:msub><mml:mi>q</mml:mi><mml:mi>H</mml:mi></mml:msub><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>+</mml:mo><mml:msub><mml:mi>q</mml:mi><mml:mi>T</mml:mi></mml:msub><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mfrac><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mo>−</mml:mo><mml:msub><mml:mi>q</mml:mi><mml:mi>T</mml:mi></mml:msub><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mrow><mml:msub><mml:mi>q</mml:mi><mml:mi>H</mml:mi></mml:msub><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:munderover><mml:mrow><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mstyle></mml:mrow></mml:semantics></mml:math></disp-formula> <disp-formula id="bmcM5"><label>(5)</label><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M8" name="1475-2875-7-100-i6" overflow="scroll"><mml:semantics><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mi>q</mml:mi><mml:mi>H</mml:mi></mml:msub><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:mi>q</mml:mi></mml:mtd><mml:mtd><mml:mrow><mml:mtext>if</mml:mtext></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mi>D</mml:mi><mml:mo>'</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>&lt;</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mi>f</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>D</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">)</mml:mo><mml:mi>q</mml:mi></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mtext>if</mml:mtext></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mi>D</mml:mi><mml:mo>'</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>≥</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>q</mml:mi><mml:mi>T</mml:mi></mml:msub><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:mi>q</mml:mi></mml:mtd><mml:mtd><mml:mrow><mml:mtext>if</mml:mtext></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mi>D</mml:mi><mml:mo>'</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>&gt;</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mi>f</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>D</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">)</mml:mo><mml:mi>q</mml:mi></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mtext>if</mml:mtext></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mi>D</mml:mi><mml:mo>'</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>≤</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:semantics></mml:math></disp-formula> <disp-formula id="bmcM6"><label>(6)</label><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M9" name="1475-2875-7-100-i7" overflow="scroll"><mml:semantics><mml:mrow><mml:mi>f</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>D</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mfrac><mml:mrow><mml:mi>D</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mrow><mml:mi>max</mml:mi><mml:mo>⁡</mml:mo><mml:mo stretchy="false">[</mml:mo><mml:mi>D</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">]</mml:mo></mml:mrow></mml:mfrac></mml:mrow></mml:semantics></mml:math></disp-formula> <disp-formula id="bmcM7"><label>(7)</label><italic>D</italic>(<italic>t</italic>) = |<italic>x</italic><sub><italic>t</italic>+<italic>q </italic></sub>- <italic>x</italic><sub><italic>t</italic>-<italic>q</italic></sub>|</disp-formula> <disp-formula id="bmcM8"><label>(8)</label><italic>D</italic>'(<italic>t</italic>) = <italic>D</italic>(<italic>t </italic>+ 1) - <italic>D</italic>(<italic>t</italic>)</disp-formula> <disp-formula id="bmcM9"><label>(9)</label><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M10" name="1475-2875-7-100-i8" overflow="scroll"><mml:semantics><mml:mrow><mml:msubsup><mml:mover accent="true"><mml:mi>σ</mml:mi><mml:mo>^</mml:mo></mml:mover><mml:mi>t</mml:mi><mml:mn>2</mml:mn></mml:msubsup><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mi>q</mml:mi><mml:mi>T</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>q</mml:mi><mml:mi>H</mml:mi></mml:msub></mml:mrow></mml:munderover><mml:mrow><mml:mo>{</mml:mo><mml:msub><mml:mi>z</mml:mi><mml:mi>t</mml:mi></mml:msub><mml:mo>−</mml:mo><mml:mover accent="true"><mml:mi>z</mml:mi><mml:mo>¯</mml:mo></mml:mover><mml:mo>}</mml:mo></mml:mrow></mml:mstyle></mml:mrow><mml:mrow><mml:msub><mml:mi>q</mml:mi><mml:mi>T</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>q</mml:mi><mml:mi>H</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mrow></mml:semantics></mml:math></disp-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M11" name="1475-2875-7-100-i9" overflow="scroll"><mml:semantics><mml:mrow><mml:msubsup><mml:mover accent="true"><mml:mi>σ</mml:mi><mml:mo>^</mml:mo></mml:mover><mml:mi>t</mml:mi><mml:mn>2</mml:mn></mml:msubsup></mml:mrow></mml:semantics></mml:math></inline-formula> <disp-formula id="bmcM10"><label>(10)</label><italic>y</italic><sub><italic>t </italic></sub>= <italic>μ </italic>+ <italic>φ</italic><sub>1</sub>(<italic>y</italic><sub><italic>t</italic>-1 </sub>- <italic>μ</italic>) + <italic>φ</italic><sub>12</sub>(<italic>y</italic><sub><italic>t</italic>-12 </sub>- <italic>μ</italic>) + <italic>β</italic><sub>1</sub><italic>P</italic><sub><italic>t</italic>-2 </sub>+ <italic>β</italic><sub>2</sub><italic>P</italic><sub><italic>t</italic>-29 </sub>+ <italic>α</italic><sub>1</sub><italic>T</italic><sub><italic>t</italic>-3</sub>+ <italic>α</italic><sub>2</sub><italic>T</italic><sub><italic>t</italic>-12 </sub>+ <italic>ε</italic><sub><italic>t</italic></sub></disp-formula> <disp-formula id="bmcM11"><label>(11)</label><italic>y</italic><sub><italic>t </italic></sub>= <italic>μ </italic>+ <italic>φ</italic><sub>1</sub>(<italic>y</italic><sub><italic>t</italic>-1 </sub>- <italic>μ</italic>) + <italic>φ</italic><sub>12</sub>(<italic>y</italic><sub><italic>t</italic>-12 </sub>- <italic>μ</italic>) + <italic>β</italic><sub>1</sub><italic>P</italic><sub><italic>t</italic>-9 </sub>+ <italic>α</italic><sub>1</sub><italic>T</italic><sub><italic>t</italic>-10 </sub>+ <italic>ε</italic><sub><italic>t</italic></sub></disp-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M12" name="1475-2875-7-100-i10" overflow="scroll"><mml:semantics><mml:mrow><mml:msubsup><mml:mi>σ</mml:mi><mml:mi>ε</mml:mi><mml:mn>2</mml:mn></mml:msubsup></mml:mrow></mml:semantics></mml:math></inline-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M13" name="1475-2875-7-100-i11" overflow="scroll"><mml:semantics><mml:mover accent="true"><mml:mi>μ</mml:mi><mml:mo>^</mml:mo></mml:mover></mml:semantics></mml:math></inline-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M14" name="1475-2875-7-100-i12" overflow="scroll"><mml:semantics><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>α</mml:mi><mml:mo>^</mml:mo></mml:mover><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:semantics></mml:math></inline-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M15" name="1475-2875-7-100-i13" overflow="scroll"><mml:semantics><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>α</mml:mi><mml:mo>^</mml:mo></mml:mover><mml:mn>2</mml:mn></mml:msub></mml:mrow></mml:semantics></mml:math></inline-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M16" name="1475-2875-7-100-i14" overflow="scroll"><mml:semantics><mml:mrow><mml:msubsup><mml:mover accent="true"><mml:mi>σ</mml:mi><mml:mo>^</mml:mo></mml:mover><mml:mi>ε</mml:mi><mml:mn>2</mml:mn></mml:msubsup></mml:mrow></mml:semantics></mml:math></inline-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M17" name="1475-2875-7-100-i11" overflow="scroll"><mml:semantics><mml:mover accent="true"><mml:mi>μ</mml:mi><mml:mo>^</mml:mo></mml:mover></mml:semantics></mml:math></inline-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M18" name="1475-2875-7-100-i12" overflow="scroll"><mml:semantics><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>α</mml:mi><mml:mo>^</mml:mo></mml:mover><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:semantics></mml:math></inline-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M19" name="1475-2875-7-100-i14" overflow="scroll"><mml:semantics><mml:mrow><mml:msubsup><mml:mover accent="true"><mml:mi>σ</mml:mi><mml:mo>^</mml:mo></mml:mover><mml:mi>ε</mml:mi><mml:mn>2</mml:mn></mml:msubsup></mml:mrow></mml:semantics></mml:math></inline-formula>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p>Time Series for The El Niño Southern Oscillation: <bold>A </bold>SST 1+2, <bold>B </bold>SST 3, <bold>C </bold>SST 3.4, <bold>D </bold>SST 4.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional file 2</title><p>Breakpoints for the rate of slide examination <bold>A </bold>F statistic for the falciparum malaria rate, the solid line is the 95% upper tail of the distribution for the F Statistic [##UREF##9##24##,##UREF##10##25##]<bold>B </bold>Empirical fluctuation period of the CUSUM test, the maximum value is in may 1992, the redline is the threshold value for breakpoint significance [##UREF##8##23##,##UREF##9##24##]<bold>C </bold>Variance of the Kolmogorov-Zurbenko adaptive filter, the dashed line is the 95% upper tail of the exponential distribution for this statistic, the maximum value corresponds to may 1992.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional file 3</title><p>Breakpoints for the <italic>Plasmodium falciparum </italic>and <italic>P. vivax </italic>rates <bold>A </bold>F statistic for the falciparum malaria rate, the solid line is the 95% upper tail of the distribution for the F Statistic [##UREF##9##24##,##UREF##10##25##]<bold>B </bold>Empirical fluctuation period of the CUSUM test, the maximum value is in may 1992 <bold>C </bold>Variance of the Kolmogorov-Zurbenko adaptive filter, the dashed line is the 95% upper tail of the exponential distribution for this statistic, the maximum value corresponds to may 1992 <bold>D </bold>F statistic for the vivax malaria rate, the solid line is the 95% upper tail of the distribution for the F Statistic [##UREF##9##24##,##UREF##10##25##]<bold>E </bold>Empirical fluctuation period of the CUSUM test for, the maximum value is in January 1992 <bold>F </bold>Variance of the Kolmogorov-Zurbenko adaptive filter. In <bold>A</bold>, <bold>B </bold>and <bold>D, E </bold>the redline is the threshold value for breakpoint significance [##UREF##8##23##,##UREF##9##24##]. In <bold>C </bold>and <bold>F </bold>the blue dashed line is the 95% upper tail of the exponential distribution for this statistic, the maximum value corresponds to December 1992.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S4"><caption><title>Additional file 4</title><p>Breakpoints for <bold>A,B,C </bold>Temperature and <bold>D,E,F </bold>Rainfall in Vanuatu using the F statistics(<bold>A,D</bold>), the empirical fluctuation period of the CUSUM (<bold>B,E</bold>) and the Kolmogorov-Zurbenko Adaptive Filter, KZAF (<bold>C,F</bold>). For the F statistics the 30% percent of the data belonging to the extremes was left out. For the KZAF only the seasonality was filtered (i.e., q = 6)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S5"><caption><title>Additional file 5</title><p>Cross-wavelet coherency and phase of <italic>Plasmodium vivax </italic>malaria rate with <italic>P. falciparum </italic>malaria rate. For technical details see legend of Figure ##FIG##3##4##.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S6"><caption><title>Additional file 6</title><p>Model selection and parameter values for models of <italic>Plasmodium falciparum </italic>and <italic>Plasmodium vivax </italic>rates before and after the breakpoint found using the Kolgomorov Zurbenko Adaptive filter.</p></caption></supplementary-material>
<table-wrap-foot><p>. % Reduction is defined as 1- (parameter value before breakpoint/parameter value after breakpoint). For <italic>P. falciparum </italic>the final model was:</p><p><italic>y</italic><sub><italic>t </italic></sub>= <italic>μ </italic>+ <italic>φ</italic><sub>1</sub>(<italic>y</italic><sub><italic>t</italic>-1 </sub>- <italic>μ</italic>) + <italic>φ</italic><sub>12</sub>(<italic>y</italic><sub><italic>t</italic>-12 </sub>- <italic>μ</italic>) + <italic>α</italic><sub>1</sub><italic>T</italic><sub><italic>t</italic>-3 </sub>+ <italic>α</italic><sub>2</sub><italic>T</italic><sub><italic>t</italic>-12 </sub>+ <italic>ε</italic><sub><italic>t </italic></sub>and for <italic>P.vivax</italic>:</p><p><italic>y</italic><sub><italic>t </italic></sub>= <italic>μ </italic>+ <italic>φ</italic><sub>1</sub>(<italic>y</italic><sub><italic>t</italic>-1 </sub>- <italic>μ</italic>) + <italic>φ</italic><sub>12</sub>(<italic>y</italic><sub><italic>t</italic>-12 </sub>- <italic>μ</italic>) + <italic>α</italic><sub>1</sub><italic>T</italic><sub><italic>t</italic>-10 </sub>+ <italic>ε</italic><sub><italic>t </italic></sub>Model selection process and all parameter values can be seen in Additional file 6. Column P (&lt;0.05) indicates the significance of any parameter B (before breakpoint)/A(after breakpoint)</p></table-wrap-foot>
<graphic xlink:href="1475-2875-7-100-1"/> <graphic xlink:href="1475-2875-7-100-2"/> <graphic xlink:href="1475-2875-7-100-3"/> <graphic xlink:href="1475-2875-7-100-4"/>
<media xlink:href="1475-2875-7-100-S1.eps" mimetype="application" mime-subtype="postscript"><caption><p>Click here for file</p></caption></media> <media xlink:href="1475-2875-7-100-S2.eps" mimetype="application" mime-subtype="postscript"><caption><p>Click here for file</p></caption></media> <media xlink:href="1475-2875-7-100-S3.eps" mimetype="application" mime-subtype="postscript"><caption><p>Click here for file</p></caption></media> <media xlink:href="1475-2875-7-100-S4.eps" mimetype="application" mime-subtype="postscript"><caption><p>Click here for file</p></caption></media> <media xlink:href="1475-2875-7-100-S5.eps" mimetype="application" mime-subtype="postscript"><caption><p>Click here for file</p></caption></media> <media xlink:href="1475-2875-7-100-S6.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media>
[{"surname": ["May"], "given-names": ["RM"], "article-title": ["Thresholds and breakpoints in ecosystems with a multiplicity of stable states"], "source": ["Nature"], "year": ["1977"], "volume": ["269"], "fpage": ["471"], "lpage": ["477"], "pub-id": ["10.1038/269471a0"]}, {"surname": ["Awerbuch", "Kiszewski", "Levins", "Martens P, McMichael AJ"], "given-names": ["T", "AE", "R"], "source": ["Surprise, nonlinearity and Complex behaviour Environmental Change, Climate and Health"], "year": ["2002"], "publisher-name": ["Cambridge University Press"], "fpage": ["96"], "lpage": ["119"]}, {"surname": ["Ricker"], "given-names": ["WE"], "article-title": ["Big effects from small causes: 2 examples from fish population dynamics"], "source": ["J Fish Res Board Can"], "year": ["1963"], "volume": ["20"], "fpage": ["257"], "lpage": ["264"]}, {"surname": ["Holling"], "given-names": ["CS"], "article-title": ["Resilience and stability of ecological systems"], "source": ["Ann Rev Ecol Sys"], "year": ["1973"], "volume": ["4"], "fpage": ["1"], "lpage": ["23"], "pub-id": ["10.1146/annurev.es.04.110173.000245"]}, {"surname": ["Buxton", "Hopkins"], "given-names": ["PA", "GHE"], "source": ["Researches in Polynesia and Melanesia An Account of investigations in Samoa, Tonga, The Ellice group and The New Hebrides, in 1925 Parts I-IV (relating principally to medical entomology)"], "year": ["1924"], "publisher-name": ["London: Memoir of the London School of Tropical Medicine and Hygiene"]}, {"surname": ["Bastien"], "given-names": ["P"], "article-title": ["Particularit\u00e9s \u00e9pid\u00e9miologiques des acc\u00e8s pernicieux \u00e0 "], "italic": ["Plasmodium falciparum "], "source": ["Med Trop"], "year": ["1987"], "volume": ["47"], "fpage": ["125"], "lpage": ["131"]}, {"surname": ["Mouchet"], "given-names": ["J"], "source": ["Le Paludisme au Vanuatu, Rapport de Mission Sante No5"], "year": ["1997"], "publisher-name": ["Port Vila: ORSTOM"]}, {"article-title": ["Enso Indices"]}, {"surname": ["Mitchell", "Hulme", "New"], "given-names": ["TD", "M", "M"], "article-title": ["Climate data for political areas"], "source": ["Area"], "year": ["2002"], "volume": ["34"], "fpage": ["109"], "lpage": ["112"], "pub-id": ["10.1111/1475-4762.00062"]}, {"article-title": ["Tyndall Center"]}, {"surname": ["Zeileis"], "given-names": ["A"], "article-title": ["A unified approach to structural change tests based on ML scores, F statistics and OLS residuals"], "source": ["Econ Rev"], "year": ["2005"], "volume": ["24"], "fpage": ["445"], "lpage": ["466"], "pub-id": ["10.1080/07474930500406053"]}, {"surname": ["Zeileis", "Kleiber", "Kr\u00e4mer", "Hornik"], "given-names": ["A", "Ch", "W", "K"], "article-title": ["Testing and dating structural changes in practice"], "source": ["Comp Stat Dat Anal"], "year": ["2003"], "volume": ["44"], "fpage": ["109"], "lpage": ["123"], "pub-id": ["10.1016/S0167-9473(03)00030-6"]}, {"surname": ["Zurbenko", "Porter", "Rao", "Ku", "Gui", "Eskridge"], "given-names": ["I", "PS", "ST", "JY", "R", "RE"], "article-title": ["Detecting discontinuities in time series of upper-air data: Development and demonstration of an adaptive filter technique"], "source": ["J Clim"], "year": ["1996"], "volume": ["9"], "fpage": ["3548"], "lpage": ["3560"], "pub-id": ["10.1175/1520-0442(1996)009<3548:DDITSO>2.0.CO;2"]}, {"surname": ["Hansen"], "given-names": ["BE"], "article-title": ["Tests for parameter instability in regressions with I(1) processes"], "source": ["J Bus Econ Stat"], "year": ["1992"], "volume": ["10"], "fpage": ["321"], "lpage": ["335"], "pub-id": ["10.2307/1391545"]}, {"surname": ["Hansen"], "given-names": ["BE"], "article-title": ["Approximate asymptotic p values for structural changes tests"], "source": ["J Bus Econ Stat"], "year": ["1997"], "volume": ["15"], "fpage": ["60"], "lpage": ["67"], "pub-id": ["10.2307/1392074"]}, {"surname": ["Bai", "Perron"], "given-names": ["J", "P"], "article-title": ["Estimating and testing linear models with multiple linear models"], "source": ["Econometrica"], "year": ["1998"], "volume": ["66"], "fpage": ["47"], "lpage": ["78"], "pub-id": ["10.2307/2998540"]}, {"surname": ["Bai", "Perron"], "given-names": ["J", "P"], "article-title": ["Computation and analysis of multiple structural change models"], "source": ["J Appl Econometrics"], "year": ["2003"], "volume": ["18"], "fpage": ["1"], "lpage": ["22"], "pub-id": ["10.1002/jae.659"]}, {"surname": ["Shumway", "Stoffer"], "given-names": ["RH", "DS"], "source": ["Time series analysis and its applications"], "year": ["2000"], "publisher-name": ["New York: Springer"], "fpage": ["572"]}, {"surname": ["Ploberger", "Kr\u00e4mer"], "given-names": ["W", "W"], "article-title": ["The CUSUM test with OLS residuals"], "source": ["Econometrica"], "year": ["1992"], "volume": ["60"], "fpage": ["271"], "lpage": ["285"], "pub-id": ["10.2307/2951597"]}, {"surname": ["Maraun", "Kurths"], "given-names": ["D", "J"], "article-title": ["Cross wavelet analysis: Significance testing and pitfalls"], "source": ["Nonlinear Proc Geophys"], "year": ["2004"], "volume": ["11"], "fpage": ["505"], "lpage": ["514"]}, {"surname": ["Levins", "Tchuenche"], "given-names": ["R", "JM"], "article-title": ["Stability in variance: a new criterion for stability"], "source": ["Far East J Dyn Sys"], "year": ["2006"], "volume": ["8"], "fpage": ["27"], "lpage": ["36"]}, {"surname": ["Vandermeer", "Yodzis"], "given-names": ["J", "P"], "article-title": ["Basin boundary collision as a model of discontinuous change in ecosystems"], "source": ["Ecology"], "year": ["1999"], "volume": ["80"], "fpage": ["1817"], "lpage": ["1827"]}, {"surname": ["Levin"], "given-names": ["SA"], "article-title": ["The problem of pattern and scale in ecology"], "source": ["Ecology"], "year": ["1992"], "volume": ["73"], "fpage": ["1943"], "lpage": ["1967"], "pub-id": ["10.2307/1941447"]}, {"surname": ["Boyd", "Kitchen"], "given-names": ["MF", "SF"], "article-title": ["Vernal vivax activity in persons simultaneously inoculated with "], "italic": ["Plasmodium vivax ", "Plasmodium falciparum"], "source": ["Am J Trop Med"], "year": ["1938"], "volume": ["18"], "fpage": ["505"], "lpage": ["514"]}, {"surname": ["Levins"], "given-names": ["R"], "article-title": ["Some demographic and genetic consequences of environmental heterogeneity for biological control"], "source": ["Bull Entomol Soc Amer"], "year": ["1969"], "volume": ["15"], "fpage": ["237"], "lpage": ["240"]}, {"surname": ["Lengeler", "Sharpe"], "given-names": ["C", "B"], "source": ["Indoor residual spraying and insecticide-treated nets"], "year": ["2003"], "publisher-name": ["Washington: Global Health Council"], "fpage": ["17"], "lpage": ["24"]}, {"surname": ["Lehane"], "given-names": ["M"], "source": ["The biology of blood sucking in insects"], "year": ["2005"], "publisher-name": ["Cambridge University Press"], "fpage": ["321"]}, {"surname": ["Chaves", "Cohen", "Pascual", "Wilson"], "given-names": ["LF", "JM", "M", "ML"], "article-title": ["Social exclusion modifies the effects of deforestation and climatic variability on a vector-borne disease"], "source": ["PLoS Neg Trop Dis"], "year": ["2008"], "volume": ["2"], "fpage": ["e176"], "pub-id": ["10.1371/journal.pntd.0000176"]}, {"surname": ["Kumar", "Hwang"], "given-names": ["R", "JS"], "article-title": ["Larvicidal efficiency of aquatic predators: A perspective for mosquito biocontrol"], "source": ["Zool Stud"], "year": ["2006"], "volume": ["45"], "fpage": ["447"], "lpage": ["466"]}, {"surname": ["Molez"], "given-names": ["JF"], "article-title": ["Les mythes repr\u00e9sentant la transmission palustre chez les Indiens d'Amazonie et leurs rapports avec deux modes de transmission rencontr\u00e9s en for\u00eat"], "source": ["Sante"], "year": ["1999"], "volume": ["9"], "fpage": ["159"], "lpage": ["162"]}, {"surname": ["Morales", "Perfecto"], "given-names": ["H", "I"], "article-title": ["Traditional knowledge and pest management in the Guatemalan highlands"], "source": ["Agr Hum Val"], "year": ["2000"], "volume": ["17"], "fpage": ["49"], "lpage": ["63"], "pub-id": ["10.1023/A:1007680726231"]}, {"surname": ["He", "Barnston", "Hilton"], "given-names": ["Y", "AG", "AC"], "source": ["A Precipitation Climatology for Stations in the Tropical Pacific Basin; Effects of ENSO"], "year": ["1997"], "publisher-name": ["NCEP/Climate Prediction Center Atlas No. 5, U. S. Dept. of Commerce"]}, {"surname": ["He", "Barnston"], "given-names": ["YX", "AG"], "article-title": ["Long-lead forecsts of seasonal precipitation in the tropical pacific islands using CCCA"], "source": ["J Climate"], "year": ["1996"], "volume": ["9"], "fpage": ["2020"], "lpage": ["2035"], "pub-id": ["10.1175/1520-0442(1996)009<2020:LLFOSP>2.0.CO;2"]}, {"surname": ["Barnston", "Chelliah", "Goldenberg"], "given-names": ["AG", "M", "SB"], "article-title": ["Documentation of a highly ENSO-related SST region in the equatorial Pacific"], "source": ["Atmosphere-Ocean"], "year": ["1997"], "volume": ["35"], "fpage": ["367"], "lpage": ["383"]}, {"surname": ["Chase", "Knight"], "given-names": ["JM", "TM"], "article-title": ["Drought-induced mosquito outbreaks in wetlands"], "source": ["Ecol Lett"], "year": ["2003"], "volume": ["6"], "fpage": ["1017"], "lpage": ["1024"], "pub-id": ["10.1046/j.1461-0248.2003.00533.x"]}, {"surname": ["Levins"], "given-names": ["R"], "article-title": ["Toward an integrated epidemiology"], "source": ["Trends Ecol Evol"], "year": ["1995"], "volume": ["10"], "fpage": ["304"], "pub-id": ["10.1016/0169-5347(95)90043-8"]}]
Malar J. 2008 Jun 2; 7:100