Forecasting Dengue Hotspots Associated With Variation in Meteorological Parameters Using Regression and Time Series Models

For forecasting the spread of dengue, monitoring climate change and its effects specific to the disease is necessary. Dengue is one of the most rapidly spreading vector-borne infectious diseases. This paper proposes a forecasting model for predicting dengue incidences considering climatic variability across nine cities of Maharashtra state of India over 10 years. The work involves the collection of five climatic factors such as mean minimum temperature, mean maximum temperature, relative humidity, rainfall, and mean wind speed for 10 years. Monthly incidences of dengue for the same locations are also collected. Different regression models such as random forest regression, decision trees regression, support vector regress, multiple linear regression, elastic net regression, and polynomial regression are used. Time-series forecasting models such as holt's forecasting, autoregressive, Moving average, ARIMA, SARIMA, and Facebook prophet are implemented and compared to forecast the dengue outbreak accurately. The research shows that humidity and mean maximum temperature are the major climate factors and exhibit strong positive and negative correlation, respectively, with dengue incidences for all locations of Maharashtra state. Mean minimum temperature and rainfall are moderately positively correlated with dengue incidences. Mean wind speed is a less significant factor and is weakly negatively correlated with dengue incidences. Root mean square error (RMSE), mean absolute error (MAE), and R square error (R²) evaluation metrics are used to compare the performance of the prediction model. Random Forest Regression is the best-fit regression model for five out of nine cities, while Support Vector Regression is for two cities. Facebook Prophet Model is the best fit time series forecasting model for six out of nine cities. Based on the prediction, Mumbai, Thane, Nashik, and Pune are the high-risk regions, especially in August, September, and October. The findings exhibit an effective early warning system that would predict the outbreak of other infectious diseases. It will help the relevant authorities to take accurate preventive measures.

Introduction

Climate change is variations in climate variables such as temperature, humidity, precipitation, rainfall, wind speed, etc. Climate Change occurs due to natural activities such as variations in the sun, volcanic explosions, or human activities like the emission of carbon dioxide and other greenhouse gases that cause global warming. Infectious diseases are categorized into foodborne, airborne, waterborne, and vector-borne infectious diseases. Vector-borne infectious diseases are transmitted to humans by a microbe, called vectors, such as mosquitoes, ticks, flies, etc. Dengue is a vector-borne infectious disease carried by mosquito vectors that is most susceptible to meteorological conditions. According to WHO, this pandemic spreads over 128 countries across the globe and increased eight times over the last 20 years affecting 4.2 million people in the year 2020. Figure 1 shows the effects of climate change on disease vectors.

FIGURE 1

Figure 1. Effects of climate change on disease vectors.

For understanding the spread of dengue, studying climate change and its effects specific to the disease is necessary. Temperature, rainfall, humidity, wind speed are the significant meteorological factors for the spread of dengue fever. Identifying the relationship between variation in these climatic factors and dengue incidences helps to predict the disease outbreak more accurately. An association has been found between the climatic parameters and dengue incidences for the selected locations in the proposed system. Machine Learning plays a vital role in developing a predictive model to understand the influx of dengue. Previously different classification and regression techniques were implemented for the prediction of dengue outbreaks for different locations. Considering varied geographic topography with changing climatic conditions and frequent disease outbreaks in the past, there is a need for better and accurate predictive models for early surveillance systems and improved prevention strategies. The following points highlight the paper's significant contributions:

• To collect climate and dengue incidence data for the selected locations for the past 10 years

• To identify the correlation between variations in climatic parameters and dengue incidences

• To implement various predictive models and show a comparative analysis based on different evaluation metrics

• To predict different climatic regions at risk in the future based on its climatic conditions

The following sections of the paper are organized as follows: Section Related Work describes related research work carried out for identifying the relationship between climate factors and vector-borne infectious diseases along with existing predictive modeling techniques available and its limitations. Section Proposed Work narrates the proposed system for dengue forecasting with variation in climate change. Section Methodology exhibits the methodologies used for forecasting dengue disease outbreaks. It includes different subsections such as data collection and integration, data preprocessing, exploratory data analysis, model execution, and evaluation metrics. Section Results and Discussion discusses details of predictive analysis and results. Finally, Section Conclusion, Limitations, and Future Work presents and concludes the author's research work.

Related Work

Significant research has been carried out for understanding the association of meteorological variables with dengue incidences. This section describes the existing work related to the prediction of dengue incidences based on climatic factors using several machine learning techniques with its advantages and limitations.

Salim et al. (1) proposed a study to predict dengue outbreaks based on weekly dengue incidence data for the Selangor state of Malaysia. Several Machine Learning algorithms such as CART, ANN, SVM, and Naive Bayes create a predictive model. It has been found that the support vector machine model (SVM) best predicted dengue outbreaks. This research provides categorical output instead of continuous output. Liu et al. (2) implemented a unique approach for forecasting dengue incidences in Guangzhou, China. They integrated SVM-MLP machine learning approaches (3) with environmental features such as water collection sites, dustbins, etc. It performs better than models based on standard features (Temperature) alone. More standard features in addition to temperature and rainfall could be considered for better training of the ML Model. The SVR-based model Tanawi et al. (4) is proposed to predict dengue incidences in DKI, Jakarta. They concluded that SVR with a linear kernel provides better results than SVR with a radial kernel. Recently Mudele et al. (5) proposed a technique that uses a recurrent neural network (RNN) for forecasting the dengue mosquito vector population. This model is compared with random forest and k nearest neighbor for two Brazilian cities. They proposed that other deep learning models should be considered for the study (6–15). Mohapatra et. al. (16) investigated the effect of climate parameters on malaria outbreak using multilayer Perceptron and J48 classifier using WEKA tool. The results show that J48 is the most suitable model than MLP and has better accuracy and less error (RMSE). Also, temperature and humidity are more significant climate parameters than rainfall, and monsoon and post-monsoon are the peak periods for the outbreak. However, other factors such as demography, immunity within the population, society's socio-economic structure, availability of affordable public health facilities are not considered during the research (17). Cheng et al. (18) proposed distributed lag non-linear model to investigate the association between extreme weather events such as floods, heatwaves, high humidity, and dengue epidemic. The researcher implemented the model on daily dengue incidences and climate factors such as temperature, humidity, and rainfall for different cities of China. The threshold for each climate parameter is calculated, and risk for dengue outbreak is identified for the extreme weather events. The limitation of the research is that other time-variant factors such as changes in mosquito density, population movements and habits, and vector control measures are not considered for the study (19).

Xu et al. (20) analyzed dengue incidences data considering different meteorological factors. They proposed long short-term memory (LSTM) based recurrent neural network predictive model to predict monthly dengue cases using climate data for 20 Chinese cities. LSTM model shows the best performance for forecasting dengue incidences. But it is time-consuming compared to other models such as the backpropagation neural network and gradient boosting machine model. Appice et al. (22) formulated different strategies such as Auto Encoding, Window-based Data Slicing, and Cluster Analysis to discover temporal dynamics in temperature and dengue variables. They proposed a new multi-stage Machine Learning model called AutoTiC-NN (22) to find trend patterns between historical data of temperature and dengue in Mexico. The study proved that the model outperforms both in regression and time series forecasting analysis. Benedum et al. (23) compared machine learning, regression, and time-series models to forecast dengue cases and outbreaks in Peru, Puerto Rico, and Singapore. They concluded that Random Forest regression provided better results than Poisson Regression and ARIMA for short-term predictions while ARIMA was better for long-term forecasts. Nkiruka et al. (24) proposed a malaria incidence classification model (MIC) using climate parameters for six countries of Africa over 28 years. The research used k means clustering for outlier detection and the XGBoost model for classification. The proposed model is compared with other classification models such as ARIMA, SARIMA, SVM and showed the best results compared to other models.

Anno et al. (21) have integrated Spatiotemporal Hotspot analysis, RS Data, and a Machine Learning approach to develop a climate-based forecasting model to deliver early warning messages to the relevant public health authorities in Taiwan. This study uses two climate parameters (Rainfall and Temperature) to predict dengue outbreaks. Stolerman et al. (25) provide a better understanding of the long-term effects of climate conditions on the Aedes Aegypti (dengue causing mosquito) population. They have developed a new data-driven method using SVM algorithms to identify climate signatures that predict Dengue epidemics in Brazil. This research uses the binary threshold to classify epidemics/non-epidemics based on the Brazilian Ministry of Health. Two climate parameters (Frequency of precipitation and average Temperature) are used. Carvajal et al. highlighted the use of time lags of meteorological factors to predict dengue incidences. They concluded that Tree based Machine Learning methods (Random Forest, Gradient Boosting) performed better than conventional statistical techniques (GAM, SAIMAX) to predict a temporal pattern of Dengue incidences in Manila, Philippines. They also suggested that Relative Humidity is one of the most critical climate factors for their RF-LG model. All the variables are trained with keeping lag time in consideration to give an early outbreak prediction. Thus, this model cannot be used to predict an immediate output (17, 19, 26, 27). Despite continuous research, due to the varied topography of India, especially Maharashtra state having different climatic regions, there is a need to develop an accurate and enhanced predictive model for effective forecasting (2, 28–36). It will help the medical researcher and public health department promptly respond to the dengue outbreak and undertake corrective majors.

Proposed Work

Figure 2 shows a schematic overview for dengue forecasting using regression and the time series model. It includes data sources and collection for both Climate parameters and dengue incidences. This is followed by data cleaning and integration in which missing data are imputed using the mean of the month data imputation technique. Exploratory data analysis is performed to find the correlation between climate parameters and dengue incidences. Feature engineering is carried out for feature selection and handling outliers. The impact of climate change includes indirect effects such as rising sea and temperature levels, extreme weather events such as droughts, floods, heatwaves, etc. The direct impact of climate change includes endurance, reproduction, or distribution of disease vectors which may affect human health. The climatic variations help in transmitting disease pathogens that may lead to infectious diseases (5, 16–19, 26, 27, 37–42).

FIGURE 2

Figure 2. Schematic overview of the proposed system.

Furthermore, the data is then split into training and testing data sets where training data is used to train different Machine Learning models–Regression Analysis and Time Series forecasting. These models are evaluated based on three evaluation metrics–Root Mean Square Error, Mean Absolute Error, and R Square Error. The models are compared to determine which models work best for different cities based on their geographic locations. Finally, locations at risk and outbreak period are predicted. Various visualization tools and techniques are used to represent the data and results effectively.

Novelties and Contribution of the Proposed Work

The effect of the variation in climate factors with varied topography on infectious diseases such as dengue is an exciting research area. The proposed work illustrates the detailed analysis of the climate and health data for different locations of Maharashtra state of India. It includes finding a correlation between monthly climate factors such as mean minimum temperature, mean maximum temperature, mean wind speed, relative humidity, etc., with dengue incidences for different locations. These locations have diverse geographic topography and weather conditions. Based on the analysis, forecasting of dengue outbreaks is performed using time series and regression models. The performance of these models is compared using various evaluation metrics and identifies the best suitable models for the study. This research will help design an effective surveillance system that will accurately monitor and control the dengue outbreak in a timely manner.

Methodology

Figure 3 shows the detailed workflow and layered architecture for the construction of the dengue forecasting model. The following sub-section (Data Sources and Collection, Data Cleaning and Integration, Exploratory Data Analysis, Feature Engineering, Model Execution, Model Evaluation) elaborates different data sources and data collection process along with data preprocessing techniques implemented such as data imputation for missing values, climate and health data integration, feature selection, and outlier detection. Exploratory data analysis is performed to identify the correlation between climate parameters and dengue incidences using different visualization techniques such as heat maps, feature plots, etc. It also depicts different time series and regression machine learning models applied along with model evaluation metrics. Finally, dengue outbreak is forecasted for different cities of India for the next 3 years.

FIGURE 3

Figure 3. Summarized flow diagram for the forecasting model.

Data Sources and Collection

Climate Data

Maharashtra has diverse climatic regions like Kokan, Khandesh, Desh, Vidarbha, and Marathwada. Based on the intensity of the disease and varied climatic conditions, nine cities like Mumbai, Thane, Ratnagiri, Pune, Solapur, Satara, Nashik, Nagpur, and Amaravati have been selected for the study. Monthly climate data is collected from Indian Meteorological Department (IMD) for 10 years from 2009 to 2019. The parameters in consideration are attributes such as Monthly Mean Maximum Temperature (MMAX) (°C), Monthly Mean Minimum Temperature (MIN) (°C), Total Monthly Rainfall (TMRF) (mm), Relative Humidity (RH) (%), and Mean Wind Speed (MWS) (km/h). Table 1 shows different region-wise locations of Maharashtra state along with population and climatic conditions.

TABLE 1

Table 1. City wise population and weather conditions.

Health Data

The monthly dengue disease incidence data is collected from the National Vector Borne Disease Control Program (NVBDCP) for targeted cities of Maharashtra state mentioned in the climate data section for 10 years from 2009 to 2019. The data collected is in excel format, having inconsistent and missing values. The climate and disease incidence data are integrated into the CSV file for all the nine targeted cities, and data preprocessing is performed. Figure 4 shows the map of Maharashtra state with region-wise selected cities for the study.

FIGURE 4

Figure 4. Map of Maharashtra state with region wise selected cities.

Data Cleaning and Integration

To create the dataset, climate and dengue incidences data are collected and integrated. The dataset generated had inconsistent values due to the diverse nature of weather and health data. For each targeted city, a few irrelevant attributes are removed from the dataset during integration. The resulting dataset consists of missing values, especially in climate parameters. Data cleaning is performed to identify missing values. The data imputation technique is used to clean the dataset. The missing data were imputed using the mean of the Month Imputation technique. In this method, the missing values are replaced with an average of the previous values of the same month throughout different years. The mean of the month imputation function is given by:

$\begin{array}{ll}\text{Vest}=({\displaystyle \sum }\text{Tj}=1\text{Vij})/\text{T}& (1)\end{array}$

$\begin{array}{ll}\text{Vest}=(\text{Vij}1+\text{Vij}2+\text{Vij}3+\dots ..+\text{VijT})/\text{T}& (2)\end{array}$

The estimated value Vest for the missing attribute is calculated by the averaging sum of values (Vij) of the variable for the ith month of the year j, where T is the number of available data for that year. In the present study, the mean of Maximum temperature “MMAX” for August 2016 was missing in the given dataset from 2009 to 2019. The estimated value is calculated by an average of previous values of the same month throughout different years. This value was treated as a data point in place of the missing value.

Exploratory Data Analysis

Once the dataset is cleaned, exploratory data analysis is performed to analyze attributes and summarize its characteristics using statistical techniques to discover useful patterns and graphical representation. City-wise feature graphs are plotted as shown in Figures 5A–I, and it is determined that each parameter for every city has a lot of variations, and there is no fixed pattern. So the correlation between each climate parameter and dengue incidences is found for all targeted locations to check which parameters are more significant.

FIGURE 5

Figure 5. City Wise Features Plot for climate variables monthly mean minimum temperature, mean maximum temperature, Average rainfall, Relative humidity, Mean wind speed, and monthly dengue incidences for nine selected cities in Maharashtra from 2009 to 2019. (A) Amravati, (B) Mumbai, (C) Nagpur, (D) Nashik, (E) Pune, (F) Ratnagiri, (G) Satara, (H) Solapur, (I) Thane.

Pearson correlation is performed on the dataset to determine the association between climate variables and dengue incidences, and heat maps are generated for each targeted city. Pearson correlation is a parametric test that measures the degree of relationship between two variables. It is the most suitable correlation technique based on the method of covariance and deals with numeric values. The person correlation function is given by Manogaran and Lopez (7):

$\begin{array}{ll}r=\frac{{\displaystyle {\sum }_{i=1}^n}(ai-a)(bi-b)}{\sqrt{\left[{\displaystyle {\sum }_{i=1}^n}{(ai-\overline{a})}^2\right]\left[{\displaystyle {\sum }_{i=1}^n}{(bi-\overline{b})}^2\right]}}& (3)\end{array}$

Here, the Pearson correlation coefficient function is employed to determine the relationship between the climate parameters and the number of dengue cases. Climate variables are monthly mean max temperature (MMAX), mean minimum temperature (MMIN), Rainfall (TMRF), Relative Humidity (RH), Mean Wind Speed (MWS).

Correlation between climate factors and dengue incidences shows that each climate variable affects the dengue incidences differently. The mean maximum temperature (MMAX) is negatively correlated with the incidences of dengue despite the locations. This implies that as the maximum temperature decreases, incidences of dengue have increased. Mean minimum temperature (MMIN) is weakly/moderately positively correlated with dengue incidences except for Nagpur. Relative Humidity (RH) is the primary climate factor and exhibits a strong positive correlation with dengue incidences for all locations of Maharashtra state. Similarly, total monthly rainfall (TMRF) is moderately positively correlated with incidences of dengue. As humidity or rainfall is increased, cases have shown an increase for all selected cities of Maharashtra. Maximum incidences occur between June to September, where the average rainfall is between 150 and 350 mm. Mean Wind speed (MWS) is a less significant climate factor and weakly negatively correlated with dengue incidences.

Figure 6 shows city-wise graphs of the Pearson correlation of each climate parameter with the dengue incidences. These graphs are further used for feature selection based on results generated to identify the significant climates factors.

FIGURE 6

Figure 6. Correlation of monthly mean minimum temperature, mean maximum temperature, total monthly rainfall, relative humidity, mean wind speed climate parameters, and monthly dengue incidences: (A) Amravati, (B) Mumbai, (C) Nagpur, (D) Nashik, (E) Pune, (F) Ratnagiri, (G) Satara, (H) Solapur, (I)Thane.

Feature Engineering

The data quality is of utmost importance for developing a predictive model with better accuracy and faster performance. For this purpose, a few data preprocessing techniques are applied, such as outlier detection and feature selection, to improve the data quality. The meteorological data consist of extreme values for specific periods, such as extreme wind speed, rainfall, and humidity. Outliers in the dataset can reduce predictive modeling performance. So the final dataset was normalized to uniform into the same scale.

Feature or attribute selection is the process of selecting the most relevant attributes in a dataset that helps train the model faster, reduces overfitting, and improves the accuracy of the predictive model. Minimal redundancy maximum relevance feature selection technique is used for attribute selection on the dataset to select attributes with high correlation and low variance. For determining the relevant features, two measures are calculated: redundancy and relevance. The following equation is used to find the mean of logical values of each climate parameter for the selected city in terms of dengue incidences:

$\begin{array}{ll}{\overline{c}}_l=\frac{1}{n_i}{{\displaystyle \sum }}_{k=1}^{n_i}{C}_{ik}& (4)\end{array}$

Where, ${\overline{c}}_i$ : Means of climate parameter i, n_i : Number of climate parameters, C_ik: the kth value of climate parameter i.

The below equation shows a variance of the climate parameters triggered by dengue incidences:

$\begin{array}{ll}{C}_i^2=\frac{1}{n_i}{{\displaystyle \sum }}_{k=1}^{n_i}{(C_{ik}-{\overline{c}}_l)}^2& (5)\end{array}$

The minimal redundancy condition for the redundancy measure can be expressed by:

$\begin{array}{ll}\text{ min }R(D),R=\frac{1}{|d{|}^2}{{\displaystyle \sum }}_{yi;yj\in d}I(yi;yj)& (6)\end{array}$

Where, min R (D) is the minimal redundancy for redundancy measure R, |d| is the number of features in the subset of feature D, and I (yi; yj) is the mutual information between feature i and j.

The maximal relevance condition for the relevance measure can be expressed by:

$\begin{array}{ll}\text{ max }RL(D,a),RL=\frac{1}{|D|}{{\displaystyle \sum }}_{yi\in d}I(yi;a)& (7)\end{array}$

Where max RL (D, a) is the maximal relevance for relevance measure RL and target activity a and I (yi; a) is the mutual information between the feature i and target activity a.

The smaller the value of the redundancy measure, the better the criteria for selection. Similarly, the higher the value of relevance measure, the better the feature selection. After exploratory data analysis is performed on the dataset considering several feature variables such as MMIN, MMAX, TMRF, RH, and MWS, different climate variables with high redundancy and low relevance are dropped for few cities under study, as shown in Table 2.

TABLE 2

Table 2. Dropped variables.

Model Execution

Regression is a supervised learning statistical method used to estimate the relationship between a dependent and one or more independent variables to determine trends in the data. It is used in the prediction of continuous values. The following regression models are implemented to predict dengue incidences across different cities based on climate parameters in the proposed system.

Support Vector Regression

It is the regression technique used to predict continuous ordered values. Some commonly used keywords in SVR are the kernel, hyperplane, boundary line, and support vectors. The primary purpose of SVR is to consider as many data points as possible within the boundary lines, and the hyperplane (best-fit line) must contain as many data points as possible. It is easy to implement and shows high prediction accuracy with excellent generalization capability. It can handle outliers very well.

Multiple Linear Regression

It is an extension of simple linear regression that models a linear relationship between more than one independent variable and a single dependent continuous variable. It is a technique for fitting a regression line through a multidimensional space of data points.

ElasticNet Regression

Elastic net is a type of regularized linear regression that includes two well-known penalties, the L1 and L2 penalty functions. The advantage of the elastic net model is that it permits a balance of both penalties, resulting in a more excellent performance on particular tasks than a model with either one or more penalties.

Polynomial Regression

Polynomial regression is a type of linear regression that estimates the connection as an nth degree polynomial. It is an example of multiple linear regression. Because Polynomial Regression is sensitive to outliers, the existence of one or two of them can have a negative impact on the results.

Decision Tree Regression

It is a regression model that breaks down a dataset into smaller subsets forming a tree with decision nodes and leaf nodes. Decision trees are very easy to visualize and reduce the uncertainty in the prediction. However, overfitting and underfitting are common problems with decision trees. If the hyperparameters are incorrectly set, the decision tree's output can vary dramatically.

Random Forest Regression

Random forest is the most commonly used machine learning technique that gives excellent results in predicting disease incidences based on climate conditions. It comprises many decision trees, each with the same node but different inputs, resulting in various leaves. It combines the results of the average of various decision trees. Overfitting can be avoided in the model by using Random Forest regression to create random subsets of the dataset.

Along with regression, the proposed system also used time series forecasting models to predict dengue incidences. Time series data is a sequence of different data points that measure a specific variable over an ordered period. In this method, time-series data extract meaningful statistics and other data characteristics to generate forecasts of our target variable. Different time series forecasting models are applied as given below:

Holt's Forecasting: It is time series forecasting method that depicts trends and seasonality from historical data. It is simple to implement and evolve with changing business requirements.

Auto-Regressive (AR) and Moving Average (MA): Autoregression forecasting technique predicts future values using previous values in time series. It demonstrates linear relation between future and past values. It is used to forecast recurring patterns in the data. The moving average method uses an average of several past points to predict future points. In this method, short-term fluctuations and the effect of extreme values are reduced.

Auto-Regressive Integrated Moving Average (ARIMA): ARIMA model is a popular time series forecasting model which uses its lags to predict future values. It uses dependent relation between past observation and current observation. It involves subtracting recent observations from previous period observations several times. It is broken down into its subtypes to increase the accuracy of incidence predictions based on climate variability. This model does not support seasonal data.

Seasonal Auto-Regressive Integrated Moving Average (SARIMA): When seasonal components are added to the ARIMA model, then it is called SARIMA. It supports univariate time series data. Additional four seasonal elements in SARIMA are P, D, Q, and m, where P is seasonal autoregressive order, D is seasonal difference order, Q is moving average order, and m is the number of time steps (17).

Facebook Prophet Model

The Facebook prophet is a relatively new time series forecasting model developed in 2017 by the Facebook data science team as open-source software. In this model, irregular observations are permitted in the dataset as it ignores temporal data dependence. It is accurate, fast, and shows excellent performance as compared with other time series forecasting models. The prophet equation is given by:

$\begin{array}{ll}x(t)=c(t)+s(t)+h(t)+u(t)& (8)\end{array}$

Where x(t) is forecast value, c(t) is the trend, i.e., change over a long period, s(t) is the seasonality, h(t) is the effect of the holiday, u(t) is unconditional changes or error. This model works best with time series with substantial seasonal influences and historical data from multiple seasons. It is robust to outliers and handles missing values very well. This model gives the best performance for six out of nine cities to forecast dengue incidences based on climate variations in the proposed system.

Model Evaluation

Once all regression and time series forecasting models are trained, the performance of the models is evaluated using three evaluation metrics: Root Mean Square Error (RMSE), Mean Absolute Error (MAE), and R Square Error (R²). RMSE is the Standard Deviation of predicted errors. Lower RMSE values indicate better models. RMSE is evaluated by the Equation (2):

$\begin{array}{ll}\text{RMSE}=\sqrt{{{\displaystyle \sum }}_{i=1}^{\text{N}}{({\text{x}}_{\text{t}}-{\overline{\text{x}}}_{\text{t}})}^2/\text{N}}& (9)\end{array}$

Here, xt is actual dengue incidences for time t and ${\overline{x}}_t$ is the predicted number of incidences by the model.

Mean Absolute Error (MAE) is the difference between the actual values and the predicted values. Lower MAE values indicate better models. MAE is evaluated by Equation (1):

$\begin{array}{ll}MAE=1/n{{\displaystyle \sum }}_{i=1}^n|{\text{X}}_{\text{i}}-\text{X}|& (10)\end{array}$

R Square Error (R²) is also known as the coefficient of determination. It tells us how well a model fits on a dataset. It indicates how close the regression line is to the actual data. The R² value closest to 1 is considered to be the best value. The equation given below evaluates the value of R² (20):

$\begin{array}{ll}{R}^2=1-\frac{{\text{SS}}_{\text{Regression}}}{{\text{SS}}_{\text{Total}}}& (11)\end{array}$

Results and Discussion

Based on the exploratory data analysis, it was observed that each climate variable affects the dengue incidences differently. The average temp range in Maharashtra state is between 26 and 43°C. As shown in Figures 6A–H, histogram graphs generated after performing Pearson's correlation shows that mean maximum temperature (MMAX) is negatively correlated with the incidences of dengue despite the locations. This implies that as the maximum temperature decreases, incidences of dengue have increased. Mean minimum temperature (MMIN) is weakly/moderately positively correlated with dengue incidences except for Nagpur. Relative Humidity (RH) is the primary climate factor and exhibits a strong positive correlation with dengue incidences for all locations of Maharashtra state. Similarly, total monthly rainfall (TMRF) is moderately positively correlated with incidences of dengue. As humidity or rainfall is increased, cases have shown an increase for all selected cities of Maharashtra. Maximum incidences occur between June and September, where the average rainfall is between 150 and 350 mm. Mean Wind speed (MWS) is a less significant climate factor and weakly negatively correlated with dengue incidences. Table 3 shows rudimentary observations for all five climatic parameters. The performance of all the regression and time series forecasting models for each city is evaluated and compared.

TABLE 3

Table 3. Rudimentary observations.

Tables 4–12 present city-wise performance comparison for all regression and time series forecasting models. The best fit values for each metric are highlighted in bold.

TABLE 4

Table 4. Performance metrics comparison table for Amravati.

Table 5 shows that decision tree regression gives the least values for RMSE, MAE, and R² compared to other regression techniques. At the same time, the Facebook prophet gives the least values for RMSE, MAE, and R² compared to other time series models for Mumbai city of Maharashtra. Similarly, Table 6 shows that the random forest model gives the best values for all performance metrics, whereas the AR model gives the least values for RMSE (5.5) and MAE (4.28) for Nagpur city. Table 7 shows that Random forest demonstrates the best performance for metrics RMSE (24.16), MAE (18.88), and R² (0.21) for Nashik city. Table 8 shows that random forest gives the best performance for metrics RMSE (14.4), MAE (9.44), R² (0.25), and Facebook prophet gives the best performance for metrics RMSE (9.3), MAE (6.7), R² (0.64) for Pune city. From all the performance Tables 4–12 and result analysis, it has been observed that Random Forest Regression is the best-fit regression model working on five out of nine cities, i.e., Nagpur, Nashik, Pune, Ratnagiri, Satara, whereas Support Vector Regression shows the best performance on two cities, Thane and Solapur. Facebook Prophet Model is the best fit time series model that worked on six out of nine cities in time series forecasting. For the rest of the cities, various combinations of ARIMA models worked as the best fit.

TABLE 5

Table 5. Performance metrics comparison table for Mumbai.

TABLE 6

Table 6. Performance metrics comparison table for Nagpur.

TABLE 7

Table 7. Performance metrics comparison table for Nashik.

TABLE 8

Table 8. Performance metrics comparison table for Pune.

TABLE 9

Table 9. Performance metrics comparison table for Ratnagiri.

TABLE 10

Table 10. Performance metrics comparison table for Satara.

TABLE 11

Table 11. Performance metrics comparison table for Solapur.

TABLE 12

Table 12. Performance metrics comparison table for Thane.

Figure 7 shows predictions for nine targeted cities using Random forest regression, and Figure 8 shows the predictions using the Facebook prophet time series model for 36 months from the Year 2021 to 2023. A hot spot map of Maharashtra state is created, as shown in Figure 9, using Tableau to compare the average number of monthly cases across our selected cities to visualize these results. The figures show that Mumbai is the most affected city, with monthly average dengue cases going up to more than 80, while Amravati is the least affected location of Maharashtra. Other cities are ranged between 5 and 35 cases. Also, Thane, Nashik, and Pune are the cities at high risk, especially in August, September, and October.

FIGURE 7

Figure 7. Predictions for random forest regression model (X-axis: Date and Y-axis: Dengue cases).

FIGURE 8

Figure 8. Predictions for Facebook prophet model (X-axis: Date and Y-axis: Dengue cases).

FIGURE 9

Figure 9. Monthly average incidence hotspot map.

Conclusion, Limitations, and Future Work

Conclusion

This research paper proposed a framework that can predict dengue incidences across different cities of Maharashtra based on climate parameters. Different meteorological variables like MMIN, MMAX, RH, TMRF, etc., are given as input, and the number of Dengue incidences is produced as output by the proposed system. Nine cities with varied climatic conditions were selected based on geographic regions. A correlation between meteorological parameters and dengue incidences was found out. The proposed system implemented 12 different regression and time series models for the prediction of dengue outbreaks. The performance of all the models is compared using root mean square error, mean absolute error, and R square error evaluation metrics. The result analysis shows that Random Forest outperforms the other Regression models for five out of nine cities. Facebook Prophet Model is the best fit time series forecasting model for six out of nine cities. The system also predicts the high-risk geographic regions from the year 2021 to 2023. It has been observed that Mumbai, Thane, and Pune are the hot spots in Maharashtra, especially from July to October. The medical researchers, public health departments, and health geography analysts can utilize these research results to take the necessary preventive measures based on these predictions.

Limitations

The study only considers climate factors. Non-climatic factors such as the demography, immunity within the population, society's socio-economic structure, availability of affordable public health facilities, and other environmental modifications initiatives are not considered for the study. Also, there is scope to add additional time-variant factors such as changes in mosquito density, population movements and habits, and vector control measures. The study is limited to a few cities of Maharashtra state of India to analyze monthly climate and dengue incidence data due to the unavailability of weekly or daily reports that could have helped better predictions.

Future Work

The result of the research will be helpful in designing an effective surveillance system that will effectively monitor and control dengue outbreaks. An output platform like a website can be created to assess the latest climate change parameters, disease outbreaks, and future projections. Future work can involve more extreme geographic regions of India along with daily or weekly climate data analysis. Vulnerability groups such as age, gender, health status, occupation of the patients can be considered to enhance the surveillance system for better planning and preparation to avoid a future outbreak.

Data Availability Statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.

Author Contributions

SPat and Span: conceptualization, data collection, interpretation, data curation, methodology, and manuscript writing. All authors contributed to the article and approved the submitted version.

Funding

The research work has been supported by Symbiosis International (Deemed) University.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

1. Salim NAM, Wah YB, Reeves C, Smith M, Yaacob WFW, Mudin RN, et al. Prediction of dengue outbreak in Selangor Malaysia using machine learning techniques. Sci Rep. (2021) 11:939. doi: 10.1038/s41598-020-79193-2

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Semenza JC, Suk JE, Estevez V, Ebi KL, Lindgren E. Mapping climate change vulnerabilities to infectious diseases in Europe. Environ Health Persp. (2012) 120:385–92. doi: 10.1289/ehp.1103805

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Liu K, Yin L, Zhang M, Kang M, Deng AP, Li QL, et al. Facilitating fine-grained intra-urban dengue forecasting by integrating urban environments measured from street-view images. Infect Dis Poverty. (2021) 10:40. doi: 10.1186/s40249-021-00824-5

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Tanawi IN, Vito V, Sarwinda D, Tasman H, Hertono GF. Support vector regression for predicting the number of dengue incidents in DKI Jakarta. Proc Comput Sci. (2021) 179:747–53. doi: 10.1016/j.procs.2021.01.063

CrossRef Full Text | Google Scholar

5. Mudele O, Frery AC. Dengue vector population forecasting using multisource earth observation products and recurrent neural networks. IEEE J Select Topics Appl Earth Observ Remote Sens. (2021) 14:4390–404. doi: 10.1109/JSTARS.2021.3073351

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Teng Y, Bi D, Guo X. Predicting the epidemic potential and global diffusion of mosquito-borne diseases using machine learning. SSRN Electr J. (2018). doi: 10.2139/ssrn.3260785

CrossRef Full Text | Google Scholar

7. Manogaran G, Lopez D. Disease surveillance system for big climate data processing and dengue transmission. Int J Ambient Comput Intell. (2017) 8:88–105. doi: 10.4018/IJACI.2017040106

CrossRef Full Text | Google Scholar

8. Guo P, Liu T, Zhang Q, Wang L, Xiao J, Zhang Q, et al. Developing a dengue forecast model using machine learning: a case study in China. PLoS Negl Trop Dis. (2017) 11:e0005973. doi: 10.1371/journal.pntd.0005973

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Wu X, Lu Y, Zhou S, Chen L, Xu B. Impact of climate change on human infectious diseases: empirical evidence and human adaptation. Environ Int. (2016) 86:14–23. doi: 10.1016/j.envint.2015.09.007

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Amutha D, Juliet M. Impact of Climate Changes on Human Health in India. (2017). doi: 10.2139/ssrn.3071055 Available online at: https://ssrn.com/abstract=3071055 (accessed October 1, 2021).

CrossRef Full Text | Google Scholar

11. Sharma V, Kumar A, Panat L. Malaria outbreak prediction model using machine learning. Int J Adv Res Comp Eng Technol. (2015) 4:1–15. Available online at: https://www.ijert.org/malaria-outbreak-prediction-using-machine-learning

PubMed Abstract | Google Scholar

12. Lopez D, Gunasekaran M, Senthil Murugan B. Spatial big data analytics of influenza epidemic in Vellore, India. In: IEEE International Conference on Big Data. Vellore (2014).

PubMed Abstract | Google Scholar

13. Patz JA, Engelberg D, Last J. The effects of changing weather on public health. Ann Rev Public Health. (2000) 21:271–307. doi: 10.1146/annurev.publhealth.21.1.271

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Du Y-D, Wang X-W, Yang X-F, Ma W-J. Impacts of climate change on human health and adaptation strategies in South China. Adv Clim Change Res. (2013) 4:208–14. doi: 10.3724/SP.J.1248.2013.208

CrossRef Full Text | Google Scholar

15. Portier CJ, Thigpen TK, Carter SR, Dilworth CH. A human health perspective on climate change: a report outlining the research needs on the human health effects of climate change. Environ Health Persp. (2017) 621–625. doi: 10.1289/ehp.1002272

CrossRef Full Text | Google Scholar

16. Mohapatra P, Tripathi NK, Pal I, Shrestha S. Determining suitable machine learning classifier technique for prediction of malaria incidents attributed to climate of Odisha. Int J Environ Health Res. (2021) 30:1–17. doi: 10.1080/09603123.2021.1905782

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Sur A, Sah RP, Pandya S. Milk storage system for remote areas using solar thermal energy and adsorption cooling. Materials Today. (2020) 28:1764–70. doi: 10.1016/j.matpr.2020.05.170

CrossRef Full Text | Google Scholar

18. Cheng J, Bambrick H, Frentiu FD, Devine G, Yakob L, Xu Z, et al. Extreme weather events and dengue outbreaks in Guangzhou, China: a time-series quasi-binomial distributed lag non-linear model. Int J Biometeorol. (2021) 65:1033–42. doi: 10.1007/s00484-021-02085-1

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Ngabo D, Wang D, Iwendi C, Anajemba JH, Ajao LA, Biamba C. Blockchain-based security mechanism for the medical data at fog computing architecture of internet of things. Electronics. (2021) 10:2110. doi: 10.3390/electronics10172110

CrossRef Full Text | Google Scholar

20. Xu J, Xu K, Li Z, Meng F, Tu T, Xu L, et al. Forecast of dengue cases in 20 Chinese cities based on the deep learning method. Int J Environ Res Public Health. (2020) 17:453. doi: 10.3390/ijerph17020453

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Anno S, Hara T, Kai H, Lee MA, Chang Y, Oyoshi K, et al. Spatiotemporal dengue fever hotspots associated with climatic factors in Taiwan, including outbreak predictions based on machine learning. Geospatial Health. (2019) 14:771. doi: 10.4081/gh.2019.771

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Appice A, Gel YR, Iliev I, Lyubchich V, Malerba D. A multi-stage machine learning approach to predict dengue incidence: a case study in Mexico. In: IEEE Access. vol. 8 (2020). p. 52713–25. Available online at: https://aquila.usm.edu/fac_pubs/17904 (accessed October 1, 2021).

Google Scholar

23. Benedum CM, Shea KM, Jenkins HE, Kim LY, Markuzon N. Weekly dengue forecasts in Iquitos, Peru; San Juan, Puerto Rico; and Singapore. PLoS Negl Trop Dis. (2020) 14:e0008710. doi: 10.1371/journal.pntd.0008710

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Nkiruka O, Prasad R, Clement O. Prediction of malaria incidence using climate variability and machine learning. Inf Med. (2020) 22:100508. doi: 10.1016/j.imu.2020.100508

CrossRef Full Text | Google Scholar

25. Stolerman LM, Maia PD, Kutz JN. Forecasting dengue fever in Brazil: an assessment of climate conditions. PLoS ONE. (2019) 14:e0220106. doi: 10.1371/journal.pone.0220106

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Javed AR, Sarwar MU, Khan S, Iwendi C, Mittal M, Kumar N. Analyzing the effectiveness and contribution of each axis of tri-axial accelerometer sensor for accurate activity recognition. Sensors. (2020) 20:2216. doi: 10.3390/s20082216

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Kumar N, Susan S. COVID-19 pandemic prediction using time series forecasting models. In: 2020 11th International Conference on Computing, Communication and Networking Technologies. Kharagpur: IEEE (2020). p. 1–7.

PubMed Abstract | Google Scholar

28. Marinucci GD. Building resilience against climate effects—a novel framework to facilitate climate readiness in public health agencies. IJER Public Health. (2014) 11:6433–58. doi: 10.3390/ijerph110606433

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Baker RE, Mahmud AS, Metcalf CJE. Dynamic response of airborne infections to climate change: predictions for varicella. Climatic Change. (2018) 148:547–60. doi: 10.1007/s10584-018-2204-4

CrossRef Full Text | Google Scholar

30. Hathaway J. Health implications of climate change: a review of the literature about the perception of the public and health professionals. Curr Environ Health Rep. (2018) 5:197–204. doi: 10.1007/s40572-018-0190-3

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Baghdad M. Climate change and simulation of cardiovascular disease mortality: a case study of Mashhad, Iran. Iran J Public Health. (2017) 46:396–407. Available online at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5395536/

PubMed Abstract | Google Scholar

32. Pandya S, Sur A, Kotecha K. Smart epidemic tunnel: IoT-based sensor-fusion assistive technology for COVID-19 disinfection. Int J Pervas Comput Commun. (2020). doi: 10.1108/IJPCC-07-2020-0091

CrossRef Full Text | Google Scholar

33. Garg D, Goel P, Pandya S, Ganatra A, Kotecha K. A deep learning approach for face detection using YOLO. In: 2018 IEEE Punecon. Pune: IEEE (2018). p. 1–4.

PubMed Abstract | Google Scholar

34. Ghayvat H, Pandya S, Shah S, Mukhopadhyay SC, Yap MH, Wandra KH. Advanced AODV approach for efficient detection and mitigation of wormhole attack in MANET. In: 2016 10th International Conference on Sensing Technology (ICST). Nanjing (2016). p. 1–6. doi: 10.1109/ICSensT.2016.7796286

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Bansal S, Chowell G, Simonsen L, Vespignani A, Viboud C. Big data for infectious disease surveillance and modeling. J Infect Dis. (2016) 214:S375–9. doi: 10.1093/infdis/jiw400

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Iwendi C, Khan S, Anajemba JH, Mittal M, Alenezi M, Alazab M. The use of ensemble models for multiple class and binary class classification for improving intrusion detection systems. Sensors. (2020) 20:2559. doi: 10.3390/s20092559

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Munasinghe A, Premaratne HL. Towards an early warning system to combat dengue. Int J Comp Sci Electro Engg. (2013) 1:252–6. Available online at: https://www.researchgate.net/profile/Noel-Fernando/publication/328189956_2nd_International_Conference_on_Computer_Electrical_Electronics_and_Mechanical_Engineering_IMCEEME'13/links/5bbd7e4292851c7fde3760a7/2nd-International-Conference-on-Computer-Electrical-Electronics-and-Mechanical-Engineering-IMCEEME13.pdf

Google Scholar

38. Finkel AM. A “solution-focused” comparative risk assessment of conventional and synthetic biology approaches to control mosquitoes carrying the dengue fever virus. Environ Syst Decis. (2018) 38:177–97. doi: 10.1007/s10669-018-9688-3

CrossRef Full Text | Google Scholar

39. Gomide J. Dengue surveillance based on a computational model of the spatiotemporal locality of Twitter. Proceedings of the 3rd International Web Science Conference. (2011). p. 1–8. doi: 10.1145/2527031.2527049

CrossRef Full Text | Google Scholar

40. Iwendi C, Anajemba JH, Biamba C, Ngabo D. Security of things intrusion detection system for smart healthcare. Electronics. (2021) 10:1375. doi: 10.3390/electronics10121375

CrossRef Full Text | Google Scholar

41. Srivastava A, Jain S, Miranda R, Patil S, Pandya S, Kotecha K. Deep learning based respiratory sound analysis for detection of chronic obstructive pulmonary disease. PeerJ Computer Science. (2021) 7:e369. doi: 10.7717/peerj-cs.369

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Iwendi C, Khan S, Anajemba JH, Bashir AK, Noor F. Realizing an efficient IoMT-assisted patient diet recommendation system through machine learning model. In: IEEE Access. vol. 8 (2020). p. 28462–74. doi: 10.1109/ACCESS.2020.2968537

PubMed Abstract | CrossRef Full Text | Google Scholar