When climate variables improve the dengue forecasting: a machine learning approach

Dengue is a viral vector-borne infectious disease that affects many countries worldwide, infecting around 390 million people per year. The main outbreaks occur in subtropical and tropical countries. We study here the influence of climate on dengue in Natal (2016-2019), Brazil, Iquitos (2001-2012), Peru, and Barranquilla (2011-2016), Colombia. For the analysis and simulations, we apply Machine Learning (ML) techniques, especially the Random Forest (RF) algorithm. In addition, regarding a feature in the ML technique, we analyze three possibilities: only dengue cases (D); climate and dengue cases (CD); humidity and dengue cases (HD). Depending on the city, our results show that the climate data can improve or not the forecast. For instance, for Natal, D induces a better forecast. For Iquitos, it is better to use CD. Nonetheless, for Barranquilla, the forecast is better, when we include cases and humidity data. For Natal, when we use more than 64\% and less than 80\% of the time series for training, we obtain results with correlation coefficients ($r$) among 0.917 and 0.949 and mean absolute errors (MAE) among 57.783 and 71.768 for the D case in forecasting. The optimal range for Iquitos is obtained when 79\% up to 88\% of the time series is considered for training. For this case, the best case is CD, having a minimum $r$ equal to 0.850 and maximum 0.887, while values of MAE oscillate among 2.780 and 4.156. For Barranquilla, the optimal range occurs between 72\% until 82\% of length training. In this case, the better approach is HD, where the measures exhibit a minimum $r$ equal to 0.942 and a maximum 0.953, while the minimum and maximum MAE vary between 6.085 and 6.669. We show that the forecast of dengue cases is a challenging problem and climate variables do not always help. However, when we include the mentioned climate variables, the most important one is humidity.


Introduction
The number of dengue cases has increased over the past 2 decades.The world-wide reported cases in 2000 were 500,000 but have an incredible increase to 5.2 million in 2019 [1].Most of these cases were reported by countries in regions of Americas.For several countries belonging to these region, the beginning of 2024 is marked by an exponential increase in the dengue cases [2].In regions of Americas, 4,565,911 cases were reported in 2023, where individual acquires permanent immunity against it and temporary cross-immunity to other ones [8].Secondary infections increase the probability of dengue hemorrhagic fever [9].Besides the consequences for human health, this disease has significant impacts on the society and economy [10].Although the dengue incidence is spatially distributed, higher outbreaks occur in the tropical countries, mainly in South America [11].Some reasons for that are due to the environmental conditions and climate factors [12], such as temperature [13] and precipitation [14].The temperature acts directly in the life cycle and reproductive of the mosquitoes [15,16], while the rain precipitation provides containers for eggs and habitat for the mosquitoes [17,18].Moreover, the relationship between rainfall and temperature are essential for regulating the water habits.Therefore, the climate variables play an important role in the seasonality of dengue epidemics [19,20].
The correlation between climate and dengue cases has been studied by many approaches [20,21,22,23,24,25,26], among them Machine learning (ML) techniques [27].Considering a multi-stage ML approach, Appice et al. investigated the evolution of dengue cases in Mexico [28].To improve forecasting, they analyzed the temperature influence.Guo et al., by means of five ML algorithms, developed a forecast model for dengue data from China [29].Additionally, they took into account climate data, such as mean temperature, relative humidity and rainfall.They showed that the support vector regression algorithm is the best one to forecast the outbreaks in China.An ensemble neural network model was proposed to forecast dengue outbreaks based on rainfall data from San Juan, Iquitos and Ahmedabad [30].The authors used a framework available to forecast long-term cases around 52 weeks.Deep learning techniques were employed by Zhao et al. to study and forecast the dengue cases in Singapore [31].They also considered the average temperature and rainfall time series.The best performance for 2, 3 and 4 weeks in advanced forecasting was 84.61%.Other algorithms have been employed to forecast dengue cases based on climate data [32,33,34,35,36,37].
Roster et al. used ML based on meteorological variables to forecast dengue cases in Brazil [38].They considered Random Forest (RF), gradient boosting regression, multi-layer perceptron and support vector regression methods.The dengue cases and meteorological variables were monthly recorded from 2007 to 2019 and from 2005 to 2019, respectively.Training the algorithms, the best performance was obtained by the RF method.With regard to RF, Ong et al. investigated the dengue transmission in Singapore utilizing data from the spread of dengue, such as population, entomological and environmental data [39].Their results showed that RF has high accuracy to reproduce dengue cases, obtaining a correlation coefficient greater than 0.86.They demonstrated that spatial risk of dengue cases can be modeled by means of RF.For Iquitos (Peru), San Juan (Puerto Rico) and Singapore, Benedum et al. combined weather data with dengue cases to con-struct a forecast model [40].Comparing time series, regression and RF methods, they verified that the RF method has 21% and 33% of error less than regression and time series models, respectively, for near short predictions (4 up to 12 weeks).Mussumeci and Coelho also employed the RF method to forecast dengue incidence in 790 cities in Brazil.Concerning the climate data, the authors considered information about incidence cases in social networks.
In this work, we employ the RF method [41] to forecast dengue incidence in three cities localized in South America: Natal (Brazil), Iquitos (Peru) and Barranquilla (Colombia).In our simulations, we use the data delayed by up to one week to forecast the new cases.We also employ three combinations of features: i) we took into account only dengue cases (D); ii) we combine dengue and climate data (CD); iii) we utilize the data of humidity and dengue cases (HD).As an important new finding, we show that, depending on the city and the training length, the results can be improved with a given combination of features.For instance, for Natal when we consider climate variables the forecasting is not improved.For Iquitos, the strategy ii is better for forecasting, while for Barranquilla the best strategy is iii.Depending on the training length, we find an optimal region for each city where the correlation among real and simulated data increase.
The structure of this work is given by the following order: In Section 2, we describe the data acquisition processing and the RF method.Section 3 is dedicated to exhibit and extract information of each time series.Forecasting results are discussed in Section 4. Finally, our conclusions are drawn in Section 5.

Data acquisition
In this work, we consider weekly dengue cases and average week climate variables from three localities: (i) Natal (Brazil, elevation 30m, latitude −5.81 and longitude −35.25) from the 16th week of 2016 until the 52th week of 2019, with time series totalling a length equal to 193 weeks; (ii) Iquitos (Peru, elevation 106m, latitude −3.74 and longitude −73.25) from the 28th week of 2001 until the 52th week of 2012, whose time series length is equal to 597 weeks; (iii) Barranquilla (Colombia, elevation 18m, latitude 10.96 and longitude −74.79) from the 2nd week of 2011 until the 47th week of 2016, with a time series length of 307 weeks.For Natal, we extract the dengue cases from Sanchez-Gendriz et al. [42] and the climate data (precipitation, relative humidity and air temperature) from National Institute of Meteorology [43].For Iquitos, we obtain both data from the Dengue Forecasting Project Data Repository [44].The climate variables for Iquitos are: minimum and average temperature, relative and absolute humidity and precipitation.Dengue cases from Barranquilla are available on Sivigila Portal [45] and meteorological (maxima temperature, relative humidity and precipitation) in Ref. [46].
The data and codes employ in this research are available on GitHub [47].

Data processing and statistical analysis
When there are many outliers in the data set, a standardization technique helps in decreasing the error in the results.In this work, we use the Robust Scaler method, which works by subtracting the median (med(X)) of the data (X) and scaling in the interval between the 1st (Q 1 ) and the 3rd (Q 3 ) quartiles.The Robust Scaler equation is given by In this procedure, the median and interquartile range (Q 3 − Q 1 ) are stored and used on future data as the transformation applied on the forecasting.One important characteristic of time series is to know if is stationary or non-stationary.There are some techniques that are employed to answer this question.In this work, we consider the augmented Dickey-Fuller test (ADF), which belongs to the unit root test.A unit root test verify if the time series is non-stationary (null hypothesis) or not (alternative hypothesis).The ADF test can be described by ∆y t = α+βt+γy t−1 +δ 1 ∆y t−1 +...+δ p−1 ∆y t−p+1 +ǫ t , (2) where α is a constant, β is the time coefficient, γ is the coefficient that present the root (is the focus of test), p is the lag order of the first differences in autoregressive process, and ǫ is a noise term.The parameter that we focus is γ.If γ = 0, or positive, then we assume the null hypothesis, i.e., a non-stationary time series.On the other hand, if γ is negative we assume the stationarity.In addition to this analyses, we also test the null hypothesis by computing the p-value.If p-value < 0.05 we assume that the time series is stationary.We conducted the ADF test by the URCA package in R [48].

Random Forest
We adopt here the Random Forest (RF) algorithm [52], which is a supervised ML algorithm based on the ensemble learning method.It is formed by many decision trees.The ensemble learning techniques are more accurate to made predictions than individual models.It occurs due to the fact that the ensemble techniques combine many basic ML algorithms and their predictions.Given an initial dataset, the model split it in K random subsets, generating different subsets from the original dataset.These final subsets are namely terminal or leaf nodes and the intermediate subsets are called internal nodes.The prediction of the results in each terminal node is made by the average outcome of the training data.In this way, the RF algorithm generates predictions or rankings from a set of multiple decision trees, as schematically represented in Fig. 1.
By combining the outputs of the trees, the algorithm provides a consolidated and more accurate result than the result from basic ML algorithms.The strength of the algorithm lies in its ability to handle complex data sets and mitigate overfitting, making it a valuable tool for various predictive tasks via ML.
The choice of partitions in each decision tree (Fig. 1) is made by considering the minimum Mean Squared Error (MSE) where y is the average value of the partition and y i is the value of each data point within that partition.Then, it is important to use the prune method as a regularization.This method makes a balance between the lowest MSE value and its depth according to where α is a tuning parameter that is determined by cross validation and T is the number of terminal nodes [49].Our simulations are implemented in Python and we use sklearn libraries for statically analyses [50].

Error analysis
In the forecasting range we compute the error among simulated and real points by considering the absolute error (e), mean absolute error (MAE), and the error (∆E).Given the real time series denoted by y i and the simulated points x i , the absolute error in the i-week is defined as the notation e D , e CD and e HD means the absolute error computed using x i from strategy D, CD, and HD, respectively.Where the index i is omitted by economy of notation.Another measure that we use is where n is the number of weeks.
Another important question to study is whether our model leads to overestimating or underestimating cases.To extract this information, we compute then, when ∆E i < 0 the model overestimating, and ∆E i > 0 the forecasting is underestimating.

Natal (Brazil)
Time series for Natal are displayed in Fig. 2, where the panel (a) shows the dengue cases (×10 2 ), (b) precipitation (mm) and (c) relative humidity (%).The panel (d) exhibits the minimum (blue line), maximum (red line) and average temperature (black line) in o C. The bottom xaxis display the year and its relative epidemiological week, while the upper x-axis shows the number of week for the whole time series.We compute the power spectra density using the package PDSR implemented in R [51] for all the variables.Our results show that the period is approximately 1 year for all variables.Specifically, for the dengue cases we obtain a period equal to 53 ± 2 weeks.Before employing the ML algorithm, we check the stationary of the data by means of the ADF method.The result of the test is displayed in Table 1.In the ADF test present in this paper, we set the p-value equal to zero when the algorithm returns a value less than 10 −6 .As the p-value for the dengue cases is greater than 0.05 and γ is close to zero, we conclude that this time series is non-stationary.Figure 2(a) exhibits the data during the four analysed years.From the 16th week of 2016 until the end of this year, 3422 cases were reported.In 2017, 2018 and 2019, there were 4747, 15178 and 17197 cases reported, respectively.In 2018 and 2019, there were more cases than expected, then in these years outbreaks occurred.Considering the time range of the first outbreak from the 10th until the 37th week of 2018 (where more than 200 cases occur every week), there were 12442 reported cases.Using the same criteria, the second outbreak exists between the 12th until the 38th week of 2019, where 14687 cases were reported.Just these two outbreaks were responsible for 27129 infections.

Iquitos (Peru)
The time series for the data from Iquitos are displayed in Fig. 3, where the panel (a) is for dengue cases (×10), Along this time series, we observe the presence of 9 outbreaks along these 12 years.Two of them call attention due to the high amplitude.The first one occurs between the 22th and the 27th weeks of 2004 and the second between the 25th and the 35th weeks of 2010, occurring a total of infection equal to 265 and 413, respectively.

Barranquilla (Colombia)
Figure 4 displays the time series for Barranquilla.The panels (a), (b), (c) and (d) exhibit the time evolution of the dengue cases (×10 2 ), maximum air temperature o C, precipitation (mm) and relative humidity (%), respectively.The data for Barranquilla start in the first week of 2011 and ends in the 47th week of 2016 (bottom x-axis), totalling 307 weeks (top x-axis).From the spectral analyzes, the period associated with the climate variables is 1 year.The period for dengue cases in Barranquilla correspond to 52 ± 4 weeks.The time series from Barraquilla are stationary, as observed in Table 3. Again, p-values less than 10 −6 are considered equal to zero.If we consider an outbreak when more than 50 cases are reported, we see 3 outbreaks in Fig. 4(a).The first one starts in the 8th week of 2013 and ends in the 52th week of the same year, having 2425 reported cases.The second one starts in the 39th and finishes in the 53th week of 2014 where 1210 cases were reported.The last one occurs from the 43th week of 2015 and finishes in the 3rd week of 2016, with 979 cases.Due to the huge outbreaks in 2013 and 2014, they are the years with more reported cases.For instance, 655, 946, 1364 and 617 cases were reported in 2011, 2012, 2015 and 2016.However, the data show 2748 and 2737 cases in 2013 and 2014, respectively.

ML based forecasting
After analysing the general characteristics of the time series, we train the RF algorithm.We consider three ap- proaches to compare the predictive capacity of the algorithm: climate data and dengue cases (CD) in the i-th week to predict dengue cases in the (i+1)-th week; dengue cases (D) from the i-th week to predict disease cases in the (i + 1)-th week; humidity and dengue cases (HD) from the i-th week to predict the cases in the (i + 1)-th week.The data are separated into the training and testing stages according to the length of the time series.

Natal forecasting
Firstly, we employ the algorithm in Natal with training a length equal to 144 weeks and the 49 remaining weeks.In our algorithm, we consider the number of trees, maximum features and depth equals to 1000, 0.5 and 100, respectively.The criteria to choose partitions is Friedman MSE [49].For CD, we consider previous cases, minimum, maximum and average temperature, humidity and precipitation.For each feature, the algorithm gives the following importance: 0.67, 0.04, 0.03, 0.03, 0.12 and 0.08.The dengue cases and humidity are more often used by the algorithm.Figure 5(a) displays the dengue cases (black points) and our simulated results (colored lines).The training region shows the red (CD approach) and the dotted blue (D approach) line, while the test range (grey background) exhibits the magenta (CD) and green (D) lines.As expected, in the train region, we observe a good accordance among the points generated by the ML algorithm and the real data.In the test region, magnified in the panel (b), the error increases when compared with the training.In the panel (b), the new curves exhibit the absolute error (e) among the real and simulated points.The light magenta colour, denoted by e CD , is the absolute error that is associated with the CD approach, while the light green (e D ) with D approach.Our results show that the error increases in the peak region, in both approaches.The mean absolute error (MAE) for CD in the test range is 97.92 and the correlation coefficient (r) is 0.90, while for D approach these values change to 62.27 and 0.94, respectively.Considering that the peak starts in the 156th week and ends in the 173th week, the MAE before the peak is 49.24 and 42.08 for the CD and D approach.During the peak, the MAE is 166.67 and 88.34 for CD and D, respectively.After the peak the MAE is 35.44 and 52.22 for CD and D. During and before the peak, D performs better than the CD approach.This scenario changes after the peak.During the test region, the maximum value of reported cases is 1113 and the MAE in peak range represents approximately 14.9% and 7.9% of this value.Outside this range, the MAE is around 4%.The result for ∆E i in the testing range, is displayed in Fig. 5(c) by the blue line for CD and the black line for D. Before the peak, we have mostly underestimating values.After the peak, we identify a mix.This overestimation occurs due to the fact that the algorithm learns this huge peak and this information remains in its memory.
Figure 6 displays the comparison between each approach, i.e., CD (blue line), D (dotted black line) and HD (dotted red line), as a function of the training length.The panel (a) exhibits the MAE and the panel (b) shows the r computed in the remaining test region.The main results are: when we train our algorithm for a few weeks, i.e., less than 55% of the time series, the RF produces considerable error in the forecast.Increasing the training length, the RF starts to perform better in the forecast.If we consider more than 90% of our time series for training, the error increases and the correlation decreases.After 90% of training, just a few weeks remain for forecast, then we do not have enough weeks to obtain a reasonable statistic.Utilizing between 60% and 89% of the time series for training, we observe that the MAE decreases and r increases.In this range, a better performance is obtained for the D approach.For a training length inferior to 55%, the CD and HD perform better than D approach.

Iquitos forecasting
The second analyzed city is Iquitos.There we consider 333 weeks of training and 263 of forecasting.We choose this training length (56% of the whole time series) due the fact that those 333 weeks separated two, apparently, symmetric outbreaks.The result is exhibited in Fig. 7(a), where the red and dotted blue lines are related to training and the magenta and green lines to the test region for the CD and D approach, respectively.The black points exhibit the dengue cases (×10).The hyper parameters are the number of trees equal to 200, maximum features equal to 0.6, depth equal to 50 and criteria to choose partitions equal to the Absolute Error.In the CD approach, we consider the following features: Cases, minimum and average temperature, relative and absolute humidity and precipitation.The algorithm weights these features as 0.39, 0.11, 0.10, 0.13, 0.13 and 0.13, respectively, in the learning process.For this, humidity and precipitation are mostly used.In addition, we consider a very long horizon forecast (up to 263 weeks).We obtain one MAE = 4.42 and r = 0.83 in the test range for the CD situation.For D approach, these values change to 4.02 and 0.81, showing that CD performs better than D in this case.The first outbreak in the testing range has a symmetric series of outbreaks.Due to this reason, we choose this training length.Doing this consideration, we observe that RF performs a good forecast, as displayed in the magnification in Fig. 7(b).As in the previous subsection, the algorithm also fails in predicting the highest peaks.The better performance in this situations seems to be in the D approach, however, this is not true.By looking for the MAE along the time series, as exhibited in Fig. 7(c), we verify that the error associated with the CD approach (light magenta line, denoted by e CD ) is less in the outbreaks when compared with the errors associated with the D situation (light green line, denoted by e D ).Also, after the huge peak localized at the 467th week, the D approach exhibits an overestimation of cases (panel (e)) higher than the CD approach (panel (d)).
MAE and r can be improved if we train our algorithm for longer time (Fig. 8).If we take less than 55% of the time series to train the algorithm, MAE increases and r decreases.In addition, for certain values, CD (blue line) performs worse than D (black dotted line) or HD (red dotted line) approaches.For values after week 331 (55% of time series) of training, the results start to improve, in the sense of obtaining a minimum MAE and maximum r, until a threshold value that occurs when we use 90% of the time series for training.At this point, MAE and r are e e respectively equal to 3.57 and 0.88 for CD, 4.06 and 0.77 for D, and 3.53 and 0.85 for HD.We see that CD and HD perform better than D. After 90%, MAE increases and r decreases.Besides, there are enough weeks after this limit (less than 60).In this given range, the algorithm does not perform very well in forecasting.The range captures just the last outbreak of the time series, where the points oscillate very irregularly.

Barranquilla forecasting
The last analyzed city is Barranquilla.For this city our simulations show better results using the HD approach then the D approach.Figure 9(a) displays the real data for dengue cases (×10) by the black points, red and blue lines the training range, and magenta and green the testing curves for CD and HD approaches, respectively.In our simulation, we take 233 weeks for training, corresponding to 76% of the time series length and remaining 74 weeks of forecasting.We use this training length to observe the forecast of the last outbreak.For the hyper parameters the number of trees, maximum features, depth and criteria to choose partitions, we utilize 1500, 0.2, 200 and Poisson.In the CD approach, the ML uses 0.65, 0.10, 0.16 and 0.10 of the features cases, maximum temperature, relative humidity and precipitation, respectively.In the HD features, the algorithm provides 0.82 and 0.17 of importance to cases and relative humidity, respectively.For the first approach, we obtain 7.81 for MAE and 0.92 for r, while for HD, the values are 6.22 and 0.94.The testing range is amplified in Fig. 9(b), where the light magenta and green curves show e CD and e D .It is important to note that the fitting of the peak performs better than in the Natal case.Considering the peak region defined between the 250th until the 262th week, the CD approach returns one MAE equal to 15.90 and HD equal to 13.16.As also observed in Natal, the approach without all considered climate variables performs better in the peak range.From the point where we started the training until 249 the MAE for CD and HD are, respectively, 3.63 and 5.00.After 262 weeks, these values are 7.09 and 4.69.The CD approach performs better before the peak.Another characteristic that emerges after 262 weeks is that the ML overestimates the forecast values, as observed by the results in Fig. 9(c).
A comparison among the approaches CD (blue line), D (black dotted line), and HD (red dotted line) is displayed in Fig. 10, where panel (a) exhibits MAE and panel (b) e e shows r as a function of training length.As observed in the previous results, the algorithm performs better in the testing region, when we consider more than 55% of the time series as training.However, in this case, if we use more than 85% of the time series as training, the error increases and the correlation decreases.It is important to note that if we take the approach D and more than 96% of the time series as training, the correlation becomes negative.The algorithm forecast increases in the cases when the data show a decay.Comparing the methods, D performs better just when we use less than 40% of the time series as training.But in the range where we get a good correlation and small error, i.e., between 55% and 85% the best strategy is consider the HD approach.

Conclusion
In this work, we employ Random Forest (RF) machine learning (ML) technique to forecast dengue infec-r tions based on previous cases and meteorological variables.We test our approach for data from Natal (Brazil), Iquitos (Peru) and Barranquilla (Colombia).In our simulations, we use three approaches: where we consider only dengue cases (D); a combination of climate variables and dengue cases (CD); the combination of humidity and dengue cases (HD).We use the features from the i-week to forecast the cases in (i + 1)-week.We also test different delays in the three cities, however, we do not obtain improvement.For each city, we use a set of climate variables.For Natal, we utilize the average, minimum, maximum temperature, precipitation and relative humidity.We find that the algorithm uses more humidity and the dengue cases in the forecasting process.Nevertheless, for Natal, we get an improvement including the dengue cases in the prediction only.For Iquitos, we consider the relative humidity, minimum and average temperature, specific humidity and precipitation.In this case, the humidity and precipitation are used with more weight by the ML technique.The best forecasting performance for Iquitos occurs, when climate and dengue cases are included.For Barranquilla, we utilize the maximum air temperature, relative humidity and precipitation.Then the most important variable considered by the ML is the relative humidity.For the data from Barranquilla, our results show a better performance, when there is the combination of dengue and humidity data.
A common characteristic that emerges is that the forecast generated by the RF algorithm exhibits a higher error in the peak regions.Moreover, exploring the effects of training length, we find that the algorithm performs better when we use more than 55% of the time series and less than 90%.The optimal region for Natal occurs when we include more than 64% and less than 80% of the time series for training.In this range, the error decreases and the correlation increases.The best performance is generated by the D method.For this approach, we verify correlations varying between 0.917 and 0.949, and MAE between 57.783 and 71.768.The optimal range for Iquitos occurs considering between 79% and 88% of the time series as training length.The method that performs better is CD, having a minimum r equal to 0.850 and maximum equal to 0.887 while MAE oscillates between 2.780 and 4.156.For Barranquilla, the optimal range occurs between 72% and 82% of length training.The method that performs better in this range is HD with minimum r equal to 0.942 and maximum equal to 0.953 while the minimum and maximum MAE are equals to 6.085 and 6.669, respectively.The ML technique reproduces many outbreaks and after that there is no outbreak.Then, the algorithm overestimated some cases, as we observed in Barranquilla.Another situation is that there are only few points in order to do statistical analysis, as in Natal.For Iquitos, the high values of training length select a very noisy region in the time series.The region is characterized by one outbreak and many fluctuations that are not predicted by RF.
Besides the dengue vector needs proper climate conditions.Our work exhibits that there is not a preferable climate variable in the forecasting process by using the RF method.On the other hand, climate variables can increase the error and decrease the correlation in some situations.In addition, also as a novelty, we find that the humidity can be more relevant than other climate variables in the forecasting process.Our results were also tested using the classic algorithms CNN 1D, LSTM, ARIMA and SARIMA regression and also we normalize the humidity data removing the annual component.However, the RF algorithm displays a better result for our goal and the modification in the humidity did not improve our results.
Our work shows that the RF method is useful for the forecasting of new dengue cases based on meteorological characteristics.This methodology can be employed by public health organization to forecast and study control measures, because RF is fast, efficient, robust, and exhibits a high correlation and low mean absolute error in the forecasting.For future works, we plan to explore other characteristics, such as mosquito eggs, and contamination rates to better predict new cases.

Figure 1 :
Figure 1: Structure of the Random Forest algorithm.K random bootstrap samples are extracted from the data and an unpruned decision tree is fitted for each bootstrap sample.At each node, a small subset of the covariance is randomly chosen to optimise the split.The forecasting is obtained by averaging the prediction of all trees.

Figure 4 :
Figure 4: Time series from Barranquilla for (a) dengue cases, (b) precipitation, (c) relative humidity and (d) maximum air temperature.

Figure 5 :
Figure 5: (a) Dengue cases (black points) and ML forecast results for Natal.The red and dotted blue lines show the training range for CD and D approach, respectively.The test range is highlighted by the gray background, where the magenta and green lines display the results by considering the CD and D approaches, respectively.(b) Magnification of testing region.The light magenta and blue curves display the absolute error e among real data and CD and D approach, respectively.(c) Error among real and simulated data. r

Figure 6 :
Figure 6: Comparison between the CD, D and HD approaches as a function of the training length for Natal data.(a) Mean absolute error (MAE) and (b) correlation coefficient (r) as function of the length training.The blue line displays the results for CD approach, while the dotted black and red lines for D and HD, respectively.

Figure 7 :
Figure 7: (a) Iquitos dengue cases (black points) and ML forecast results.The red and dotted blue lines show the training range for the CD and D approaches, respectively.The test range is highlighted by the gray background, where the magenta and green lines display the results by considering the CD and D approaches, respectively.(b) Magnification of testing region.(c) MAE among real and simulated data, where the light magenta line is associated with CD and the light green with D. (d) Absolute error for the CD approach.(e) Absolute error for the D approach.

Figure 8 :
Figure 8: Comparison between the CD, D and HD approaches as a function of the training length for Iquitos data.(a) Mean absolute error (MAE) and (b) correlation coefficient (r) as a function of length training.The blue line displays the results for CD approach, while the dotted black and red lines for D and HD, respectively.

Figure 9 :
Figure 9: (a) Barranquilla dengue cases (black points) and ML forecast results.The red and dotted blue lines show the training range for CD and HD approaches, respectively.The test range is highlighted by the gray background, where the magenta and green lines display the results by considering the CD and HD approaches, respectively.(b) Magnification of testing region.(c) Error among real and simulated data.

Figure 10 :
Figure 10: Comparison between the CD, D and HD approaches as a function of the training length for Barranquilla data.(a) Mean absolute error (MAE) and (b) correlation coefficient (r) as a function of length training.The blue line displays the results for the CD approach, while the dotted black and red lines for D and HD approaches, respectively.

Table 1 :
Test of stationarity using ADF Test for the Natal features.More precisely, the cases has a period equal to 55 ± 3 weeks.The ADF test shows that the times series related to Iquitos are stationary, as observed in Table2.The associated p-values in Table2are setted equal zero because our test report values less than 10 −6 .
o C (black line).The considered range for these data varies from the 28th week of 2001 until the 52th week of 2012, totalling 597 weeks (marked in the top x-axis).By the power spectra density analyze, our results shows a period of 1 year for all the variables.

Table 2 :
Test of stationarity using ADF Test for the Iquitos features.

Table 3 :
Test of stationarity using ADF Test for the Barranquilla features.