 Research article
 Open Access
 Published:
Prediction mapping of human leptospirosis using ANN, GWR, SVM and GLM approaches
BMC Infectious Diseases volume 19, Article number: 971 (2019)
Abstract
Background
Recent reports of the National Ministry of Health and Treatment of Iran (NMHT) show that Gilan has a higher annual incidence rate of leptospirosis than other provinces across the country. Despite several efforts of the government and NMHT to eradicate leptospirosis, it remains a public health problem in this province. Modelling and Prediction of this disease may play an important role in reduction of the prevalence.
Methods
This study aims to model and predict the spatial distribution of leptospirosis utilizing Geographically Weighted Regression (GWR), Generalized Linear Model (GLM), Support Vector Machine (SVM) and Artificial Neural Network (ANN) as capable approaches. Five environmental parameters of precipitation, temperature, humidity, elevation and vegetation are used for modelling and predicting of the disease. Data of 2009 and 2010 are used for training, and 2011 for testing and evaluating the models.
Results
Results indicate that utilized approaches in this study can model and predict leptospirosis with high significance level. To evaluate the efficiency of the approaches, MSE (GWR = 0.050, SVM = 0.137, GLM = 0.118 and ANN = 0.137), MAE (0.012, 0.063, 0.052 and 0.063), MRE (0.011, 0.018, 0.017 and 0.018) and R^{2} (0.85, 0.80, 0.78 and 0.75) are used.
Conclusion
Results indicate the practical usefulness of approaches for spatial modelling and predicting leptospirosis. The efficiency of models is as follow: GWR > SVM > GLM > ANN. In addition, temperature and humidity are investigated as the most influential parameters. Moreover, the suitable habitat of leptospirosis is mostly within the central rural districts of the province.
Background
Since the discovery of leptospira in the body of Japanese mine workers over a hundred years ago, human leptospirosis has been treated as a “neglected tropical disease” worldwide [1]. Reports of World Health Organization show that annual incidence rate of leptospirosis per 100,000 people varies from 0.1 to 1 in temperate regions and 10–100 in humid regions and over 100 in tropical areas. Global report of the disease reveals that over 1 million severe cases take place annually with approximately 60,000 fatalities [2]. As a Zoonotic disease, it occurs in tropical and subtropical areas with high humidity [3]. This disease is caused by leptospira bacteria which live in the urine of mammals such as rodents [4]. Human infection from leptospirosis occurs through direct or indirect contact with infected animals or environment [5]. Several contributing factors are contemplated for the incidence of leptospirosis including geographical location with frequent rainfall and floods, adjacency to mammal reservoirs and human activities [6]. One of the most important reasons of leptospirosis mortality is its resemblance to other diseases such as influenza and dengue fever [7]. Indeed, underestimating its infectiousness and loss of timely diagnosis give rise to fatality [8].
Rafyi and Magami in 1968 confirmed the first report of human leptospirosis in Iran, but no definite report has been made about the current status of human leptospirosis distribution in the country [9]. Human leptospirosis, an endemic disease in Caspian region, is more widespread in Gilan Province because of humid and wet climate [3]. In addition, high population densities of rural districts, farmlands (often paddy fields) and fishing activities help propagate the prevalence of leptospirosis in Gilan. Amongst provinces, the annual incidence rate of leptospirosis in Gilan is always the highest. In this region, most farmers keep domestic animals in their houses and irrigate their farms using river resources, where the population of leptospirosiscontaminated rodents is abundant [9]. Hence, modelling and predicting leptospirosis will help policy makers to better understand the disease, prioritise regions and budget for early prevention or treatment and provide accurate planning. It will help the government policy makers ease the burden of medical and health care expenditure on the province.
Several studies were made on modelling leptospirosis worldwide [10,11,12,13]. Many studies elucidated the effect of drivers such as precipitation [14, 15], temperature [16, 17], humidity [18, 19], elevation [20] and vegetation [10, 21] on the distribution of leptospirosis because its prevalence highly depends on environmental factors. However, most studies focused on clinical aspects of the disease and animal type of leptospirosis. Based on literature review and to the best of our knowledge, papers rarely worked on spatial modelling and predicting human leptospirosis utilising Geographical Information System (GIS) and its approaches [11, 12].
GIS is a powerful tool that its capabilities have been already proven in various fields of studies such as disease [22,23,24] and environment [25,26,27]. In disease problems, GIS can play a major role in showing how the disease propagates and finding the parameters that affect its prevalence [28]. The advantages of GIS have been proven in developed countries, but it is rarely employed for health issues in developing countries such as Iran [29, 30].
Given that the heterogeneity relationship between the disease and effective parameters, some methods should be utilized to consider heterogeneity [31]. Geographically weighted regression (GWR) is a common approach that can solve the heterogeneity by considering variability of coefficients in diverse locations across the study area [32]. An advantage of GWR is considering the location of parameters as input to improve spatial prediction capability and reduce heterogeneity effect. GWR is an efficient approach for modelling in various fields of study [33,34,35], especially disease modelling and predicting. However, GWR is a linear method that cannot consider the nonlinear behaviour of the phenomenon. Owing to high capability in solving nonlinear problems, Artificial Neural Network (ANN), a widely used approach in disease prediction, is selected to predict leptospirosis disease [36,37,38,39]. Another approach used in this study is General Linear Model (GLM), which is a statistical model commonly used in modelling and predicting diseases [40]. It utilizes the polynomial regression to investigate the relationship between dependent and independent variables [41]. Also, SVM, a supervised classifier, is used as a novel machine learning method which can be used for classification and in regression analysis [42]. The SVM classifier takes a set of input dataset and predicts the class of each input data which is used in various medical issues [30, 43] .
This study aims to model and predict human leptospirosis in Gilan Province of Iran, using capabilities of GWR, GLM, SVM and ANN approaches. Background section provides knowledge about leptospirosis and the reasons of its prevalence based on previous studies. Methods section explains how data are prepared and asserts fundamentals about utilized approaches. Results section presents the results of models. Discussion section interprets data ally with analysing the information which can be obtained from the results of the models in detail. The final section describes the conclusions of the study and indicates future work.
Methods
Study area
Gilan, a northern province of Iran, ranks second in rice cultivation. Figure 1 depicts the geographical location of Gilan at 48°53′–50°34′ longitude and 36°34′–38°27′ latitude. It consists of approximately 2.531 million inhabitants, 107 rural districts and 14,042 km^{2} area. It stretches across Alborz Mountains with dense forests in the south (highlands) and Caspian Sea in the north (lowlands). In this study, modelling is performed at the rural district level, and the centroids of the rural districts are considered as the base level for analysis. These centroids are selected as the points to which all parameters are allocated. Notably, the centroids are the geometric centres of polygons of the rural districts. The mean, maximum and minimum area of the rural districts are 129,253,376 m^{2}, 441,566,882 m^{2} and 113,055,500 m^{2}, respectively.
Data acquisition and preparation
The input parameters utilized in this study are disease, climate, topography and vegetation data collected from relevant organisations in Iran (Ministry of Health and Meteorology Agency) from 2009 to 2011. The population data of rural districts used in this study are gathered from the National Centre of Statistics of Iran, and these data are updated every 5 years by this organisation across the country. The latest updated population data at the rural district level of Gilan, which included the population size of different divisions of the country separately, are used in this study. The data in 2009 and 2010 are used for modelling, and the models are assessed by the data in 2011. All data are prepared and integrated using ArcGIS 10.2 and Microsoft Excel 2010 for further analysis. To avoid very large or small weights, the input data are normalised between [0,1] using Eq. (1) [44].
where x_{i} denotes the input parameter; x_{min} and x_{max} are minimum and maximum values of x_{i}, respectively.
Disease data
All villages in Iran are covered by the wellfounded National Health Care Network (NHCN), which is sponsored by National Ministry of Health and Treatment of Iran. The disease data (positive results of ELISA^{Footnote 1} blood test of patients) are gathered from database of NHCN and Health Centres (HC) of Gilan. The spatial distribution of the disease throughout the study area is illustrated in Fig. 2.
Incidence rate measures the frequency of disease occurrence in the population over a specified time. The major advantage of calculating incidence rate is the omission of the effect of population on disease prevalence across the study area. To eliminate the effect of population on results, incidence rate is calculated using Eq. (2):
Climate data
Temperature (degree Celsius), humidity (percentage) and precipitation (millimeter) are gathered from 12 synoptic climate stations of Gilan in Excel format (.xlsx) (Fig. 3).
Given that the climate data are collected from the meteorological stations and the limited number of these stations across the study area, a continuous surface of the climate parameters is produced utilising IDW^{Footnote 2} interpolation method. The obtained maps are demonstrated in Fig. 4.
Topographic and vegetation data
Gilan shows remarkable topographic variations with almost 3700 m altitude difference between the lowest and highest locations and average altitude of 1800 above sea level. Elevation continually decreases from south to north. Owing to the significant variability of elevation, climate and vegetation differ across the study area. The elevation map is obtained from NASA^{Footnote 3}‘s 90 m resolution SRTM^{Footnote 4}data. All parameters such as elevation are assigned to the centroids of rural districts for further analysis. ArcGIS software tool ‘Extract to Points’ is employed, and the elevation data are assigned to the centroids.
Vegetation is another environmental factor which influences leptospirosis vector directly or indirectly [10]. To investigate the effect of vegetation at the rural district level, Normalised Difference Vegetation Index (NDVI) is used in this study. This process is performed using the satellite images of the Gilan and the capabilities of ENVI^{Footnote 5} software, a wellknown software in image processing. Satellite images of MODIS^{Footnote 6} during 2009–2011 are used to extract NDVI via ENVI software. Their period is 16day ally with 250 m spatial resolution. The satellite images are mosaicked, and the NDVI index of study area is subsequently calculated and used as the vegetation parameter in the study. The vegetation values of all rural districts are allocated using the calculated NDVI. Elevation data and variability of NDVI in 2009 and 2011 are presented in Fig. 5.
All parameters and their characteristics are presented in Table 1.
GWR
GWR presented by [45] is the most important regression approach in spatial modelling. The general equation of this approach is expressed as follows (Eq. (3)):
where j = 1,2,…,n shows the number of rural districts, Y_{j} is the incidence rate of leptospirosis in rural district j, (U_{j}, V_{j}) denotes the geographic location of rural district j, B_{i} is the local coefficient of parameter k, X_{jk} is the value of input parameter in rural district j, ε_{j} is the error value, and B_{k} is obtained from minimising Eq. (4):
where W_{jk} is distance decay function for location j. Three distance decay functions are applicable in GWR model, namely, Poisson, Gaussian and Logistic. In this study, Gaussian function (Eq. (5)) is used due to its higher efficiency [46]:
where d_{jk} is the spatial distance between rural district i and k, and b identifies the kernel bandwidth. Three bandwidth selection criteria, including AIC (Akaike Information Criterion), CV (Cross Validation) and BIC (Bayesian Information Criterion), and two kernels (fixed and adaptive) are available in modelling by GWR [45].
GWR model for leptospirosis prediction
To predict leptospirosis, a model is established based on environmental parameters utilising GWR approach. Five parameters, including temperature, precipitation, humidity, elevation and vegetation, in 2009 and 2010 together with disease data are used as inputs of the model. The model is used for predicting of leptospirosis in 2011.
According to the description of methods and input parameters, the GWR model is formulated as Eq. (6):
where Y_{j} denotes the incidence rate of leptospirosis (the independent parameter), B_{temp}, B_{prec}, B_{hum}, B_{elev} and B_{Veg} are the correlation coefficient values of input parameters, X_{1}, X_{2}, X_{3} and X_{4} are the values of dependent parameters in a definite rural district, and (U_{j}, V_{j}) denotes the location of rural district j.
Fixed and adaptive kernel functions are applicable for the GWR model. Fixed kernel considers a constant bandwidth (distance to neighbour in metre) across the study area, which is the main deficiency of this kernel, whereas adaptive kernel applies variable and appropriate bandwidths (number of neighbours) in each rural district according to the number of neighbours [47]. In addition to type of kernel, defining bandwidth selection criteria is necessary in the GWR model. Three bandwidth criteria of AIC, CV and BIC are available. Adaptive kernel and AIC criteria are utilized in this study due to better performance [48]. Notably, all steps are performed using GWR 4.0 software.^{Footnote 7}
ANN
ANN is a nonlinear model that focuses on determination of dependence between input and output parameters by simulating highly connected processing units (neurons) of human nervous system [49]. It consists of three layers including input, hidden and output, and it is composed of weighted connections between the inputs and outputs [50]. A major characteristic of ANN is its capability to learn for solving complex problems [51]. The other advantage of ANN is proper description of nonlinear dependences. However, the black box mechanism is its major shortcoming [52].
A particular form of ANN is Multilayer Perceptron (MLP) which is created by multiple layers of nodes in a directed graph [53]. MLPs are FeedForward Neural Networks (FFNN) that stream information in one direction from the input to the output layer. MLPs are the most popular FNNs due to efficient training processes [54].
In ANN, input data should be normalised before feeding to the model because different data with diverse ranges should be mapped into a similar range. Training data which adjust the weights of neurons and decline the model bias are also important in modelling using ANN [55]. Data training has several algorithms, and Levenberg–Marquardt algorithm is a popular one [56]. After training data, test data should be utilized to evaluate the performance of the network. Figure 6 exhibits the structure of MLP used in this research.
ANN model for leptospirosis prediction
MLP, a class of FFNN is utilized for leptospirosis prediction. MATLAB 2018 is used for MLP implementation. According to the trial and error approach (Additional file 1), one hidden layer is selected to be utilized in this study. The final MLP architecture consists of five nodes in input layer, including temperature, precipitation, humidity, elevation and vegetation, one hidden layer with five nodes and one node in output layer, which presents the incidence rate of leptospirosis. Data of 2009 and 2010 and Leungberg–Marquard algorithm are used for training the model to predict the disease in 2011. Weights are randomly initialised, and the threshold of the training process is considered when the error difference of two consecutive runs of the model is negligible. Notably, after running ANN under such condition (reaching a negligible difference of two consecutive runs), the maximum number of epochs is 36. Total sample points are 969 for 2009 and 2010 in which 290 samples are selected as validation set. The learning rate, which is acquired using trial and error approach, is 0.01.
SVM
SVM, first introduced by Vapnik [57], is a supervised classifier based on the statistical theory. In a linear situation, the basic SVM tries to maximize the distance between closest samples of binary classes by creating optimal hyperplanes [57]. However, most of the problems in real world do not behave in linear manner. In order to deal with nonlinear datasets, SVM utilizes kernel functions to map data into higher dimensional space in which the data is linearly separable [58].
Consider the input data as {x_{i}, x_{2}, …, x_{i}} named vectors and their corresponding labels as y_{i} ∈ {−1, +1}, SVM constructs hyperplanes which separate positive labels from negative ones. Equations (7) and (8) are used to investigate the label of data in nonlinear situation [59]:
Where b is the bias, K(x_{i}x_{j}) is the kernel function and α_{i} denotes the Lagrange’s multiplier which can be calculated by maximizing eq. (9). C is regularization constant which balances the maximization of sample distances and model error [60].
SVM model for leptospirosis prediction
In order to apply SVM model, input data are categorized into 5 classes (very low, low, moderate, high, and very high classes). The Data of 2009 and 2010 is used to train the SVM model and it is utilized to predict leptospirosis in 2011. Because SVM is a binary classifier, it cannot be directly used for a multiclass problem. In order to perform a multiclass classification using a binary classifier, oneagainstall method can be used to divide each multiclass classification into groups of binary classifications [57]. In this study, 5bainry SVMs are constructed (5 is the number of classes) in which, each binary classifier separates one class from the rest of the classes. Another vital step in running an SVM model is the selection of its parameter (C) and the type of kernel function [59, 61]. Leaveoneout crossvalidation method [62] is applied on training dataset to select Parameter C and the value of 2 is obtained as the best value in this study. The most common kernel functions have been used in previous studies are the linear, polynomial, and Radial Basis Function (RBF) [63]. Therefore, in order to determine the best kernel function, these functions are compared in this study and the output result is presented in Table 2. As it is shown in this table, RBF could obtain more accurate result in this study. Java programming language is used to implement SVM in this study.
Sensitivity analysis
ANN and SVM function as a black box, so investigating the relative importance of input parameters is not possible. However, sensitivity analysis can be used to examine the contribution of input parameters in modelling and predicting [64]. To perform sensitivity analysis, one parameter is excluded from the model in each run, and the effect of that parameter on model performance is determined based on the evaluation criteria [65]. A larger decrease indicates greater influence of the parameter.
GLM
The Generalized Linear Model is one the most common statistical approach identified for prediction mapping [66]. GLM assumes a relationship between the dependent variable and different independent variables given by (Eq. (10)):
where E (y) is the value of dependent variable y, X_{j} indicates j^{th} independent variables regarding to p covariates to be estimated and B_{j} is the j^{th} coefficient.
GLM model for leptospirosis prediction
GLM model is established based on the input variables in which the variables do not change locally in spite of GWR model. In this study the following model is used as the GLM model for prediction of leptospirosis:
where Ln (A) is log of disease data, X_{j}, j^{th} independent variables (j = 2, ..., p) and B_{j}, j^{th} coefficients of variables (j = 0, ..., p). The ordinary leastsquares estimates are calculated to obtain Maximumlikelihood estimates for GLM which performs like a multivariate analysis. All implementation of this approach is done using SPSS software version 23.
Spatial autocorrelation
Spatial autocorrelation is useful for analysing and examining randomness of residuals [67]. Moran’s I is commonly used for checking spatial autocorrelation and cluster detection which ranges between − 1 and 1 (Eq. (12)) [67]:
where W_{ij} is the spatial weight between i^{th} and j^{th} provinces; z_{i} and z_{j} are the values of zscore in i^{th} and j^{th} provinces, respectively; Y_{i} is the number of cases for i^{th} province; and S is the sum of all spatial weights. Moran’s I is used to determine the spatial autocorrelation of residuals for investigating the model deficiencies.
Evaluation
To assess the results of approaches, Mean Square Error (MSE), Mean Absolute Error (MAE), Mean Relative Error (MRE) and R^{2} are employed as Eqs. 13–16.
MSE is the most common statistic for regression evaluation and it is defined as follow (Eq. (13)) [68]:
where y_{i} is leptospirosis report and ŷ_{i} is its prediction. For each rural district, it calculates the average square difference between the predictions and actual values. It is useful when we have unexpected values that we should pay attention.
MAE means that all the individual differences are weighted equally in the average. It is calculated using Eq. (14) [68]:
The advantage of this statistic is that it is not sensitive to outliers as MSE. We considered the relative errors in each rural district and calculate the mean of it to obtain MRE value. Equation ((15)) represents this statistic [69]:
Realizing the performance of models is difficult when we use only MSE, MAE and MRE criteria. R^{2} is a metric has the advantage of being scalefree and can solve this issue. Many papers indicate that the range of R^{2} is between 0 and 1. Equation (16) is used to calculate R^{2} [36, 61]:
where y_{i} is the incidence rate in rural district i, \( {\hat{y}}_i \) is the predicted value, and n is the number of rural districts.
Results
According to the database of NHCN, leptospirosis occurs annually in certain months (approximately March to September) and remarkably coincides with the beginning of rice planting and end of harvest season (Fig. 7.a). Reports confirmed that in 2009 (312 cases), 2010 (657 cases) and 2011 (217 cases), 1186 positive cases were reported, and the peak of leptospirosis prevalence occurred in 2010 in Gilan, which is twice as much as last year. Amongst reported cases, 70% of patients were men who are more vulnerable to leptospirosis infection than women (Fig. 7.b).
The calculated correlations of input parameters are shown in Table 3. The table shows the maximum correlation between elevation and temperature (− 0.33), which is consistent with reality (the higher elevation, the lower the temperature); the minimum is between vegetation and temperature (0.11). Variance inflation factor (VIF) is calculated for input parameters, and the results are presented in Table 3. All VIF values are less than 2.71, confirming no severe multicollinearity amongst input parameters.
Minimum, maximum, range and standard deviation obtained from GWR model are presented in Table 4, which shows the variability of each parameter in the spatial modelling of leptospirosis.
Table 5 presents the coefficients of input parameters obtained from GLM model. They clarify the impact of each parameter on modelling leptospirosis distribution.
The output of sensitivity analysis of ANN and SVM are presented in Table 6 and Table 7. Temperature and humidity are utmost effective parameters of leptospirosis prediction because their removal leads to a decrease in the value of four criteria. On the contrary, removing vegetation and precipitation lead to improving the accuracy of prediction, which shows less effect of both parameters in prediction.
Figure 8.a shows the actual number of leptospirosis disease in 2011. Figure 8.b, 8.c, 8.d and 8.e show the results of GWR, ANN, SVM and GLM prediction in 2011, respectively. The disease rarely occurs in the southeast rural districts.
Local variability of GWR model in each rural district is shown in Fig. 9.a. The size of dots in the map illustrates the prediction accuracy of GWR model in different rural districts. Local collinearity of GWR model is examined to evaluate the fitness of model via calculating condition numbers for each rural district (Fig. 9.b).
The coefficients of GWR model are demonstrated in Fig. 10. Similarities are observed between coefficients of temperature and humidity with prediction map of 2011.
Detected clusters at 95% significance level are demonstrated in Fig. 11 for GWR, GLM and ANN models. GWR, GLM and ANN models do not perform well in leptospirosis prediction in several districts.
Discussion
During 2009–2011, reports of leptospirosis in Gilan revealed that it occurs in definite months and disappears for the remainder of the year. This periodic prevalence explains the relationship between leptospirosis cases and paddy season when workers start to work in paddy fields. This phenomenon is due to the fact that in paddy season when workers begin to plant or harvest rice, their contact with contaminated water or soil increases, and the possibility of disease prevalence increases. In Gilan, rice farming and livestock are popular amongst farmers because suitable climate contributes to the fertility of soil which is inevitable for farming, and the existence of many rural regions covered by grasslands and forests facilitates feeding animals. Considering that this job is physically demanding, the ratio of men to women workers is approximately 2 to 1 in 2009–2011, which confirms that men are more vulnerable to this disease and deserve more attention (Fig. 7.b). This fact prompted decision makers to carry out prevention programmes such as boosting the knowledge of workers by explaining the advantages of using gloves during work time or bandaging the wound as soon as it occurs. Knowledge and literacy are at low levels in rural districts, so such programmes led to a great decrease of disease reports (almost 1/3) in 2011 (Fig. 7.a).
Spatial modelling of leptospirosis would better clarify different aspects of this phenomenon. To model the disease, the correlation between input parameters should be investigated using the assumption of independence [70]. Correlation values vary from 0 (no correlation between two parameters) to 1 (maximum correlation between two parameters), and the closer the values are to 0, the more reliable they are as input in the model. Based on statistical studies about the assumption of independence, less than 0.70 correlation is acceptable [71]. Thus, twotailed Pearson correlation as a common approach [72] is used in this study to calculate the correlation amongst all parameters. According to the obtained values, maximum correlation is between elevation and temperature parameters (0.33) with 0.005 significance level, and minimum is between vegetation and temperature (0.11) with 0.1 significance level. The results prove that all values are less than critical threshold (0.70) [71] and can be reliably utilized in spatial modelling of leptospirosis (Table 3).
In addition to assumption of independence, multicollinearity should be considered in spatial modelling [73]. Severe multicollinearity increases the variance estimation of coefficients and decreases the reliability of the model. VIF measures the intensity of multicollinearity amongst independent parameters [74]. Confirmed by statistical studies, VIF values of input parameters that are less than 10 are acceptable for entering the model [75]. Table 3 presents that the maximum calculated VIF values of parameters belong to vegetation parameter (2.71), and the minimum is acquired for precipitation parameter (1.17). All VIF values are less than 10, which proves acceptable multicollinearity amongst input parameters. According to the assumption of independence and VIF values, input parameters can be fed to GWR, GLM, SVM and ANN models for predicting leptospirosis distribution in this study.
The values of coefficients calculated for each parameter using GWR and GLM are presented in Table 4 and Table 5. GWR considers a different model for each rural district, so the coefficients of parameters vary across the study area. Slight changes in the range of elevation (D_{2009} = 0.17, D_{2010} = 0.73 and D_{2011} = 0.13) and vegetation (D_{2009} = 0.09, D_{2010} = 0.14 and D_{2011} = 0.16) reveal almost uniform and constant effect of these parameters. High values of temperature, precipitation and humidity range (1740.69, 321.64 and 812.94, respectively) show inconstant effects on diverse rural districts. Despite GWR and GLM models, ANN and SVM operate as black box. The coefficients of parameters cannot be calculated, but sensitivity analysis can be utilized for this issue. The results of sensitivity analysis are presented in Table 6 and Table 7, which show the effect of parameters on spatial modelling of leptospirosis distribution. According to four evaluation criteria, omission of temperature and humidity parameters decreases the fitness of the models, which confirms their importance in modelling the disease. Temperature and humidity do not directly affect leptospirosis distribution but provide appropriate circumstances for durability of leptospira and indirectly affect the prevalence of leptospirosis. Paddy fields are almost always located in rural districts with higher values of these parameters, and they are more vulnerable to the disease occurrence, as shown in Fig. 10, where coefficients are mapped for better understanding of the effect of parameters on different rural districts. Maps of coefficients of humidity and temperature are closer to prediction maps and reports of leptospirosis data in 2011. This finding proves that these two parameters play more important roles in the modelling and predicting leptospirosis.
Prediction maps of GWR, GLM, SVM and ANN
The models clarify the fact that the disease prevalence occurs more in the central rural districts. The existing remarkable number of paddy fields and livestock activities, which leads people to more contact with the contaminated environment, can be the major reasons of this pattern. Given that leptospirosis is an occupational waterborne disease [76] and no paddy fields are in the southeast area of the province, the probability of the disease prevalence is negligible there. Visual comparison of the prediction maps shows that GWR, SVM and GLM models predict high disease prevalence in the central rural districts while the prediction of ANN model is less consistence with the reported cases of disease across the study area. Although SVM and GLM indicate satisfying results, GWR prediction map in 2011 is more similar to the map of leptospirosis data in 2011. Model predictions are statistically discussed in the “prediction evaluation” section.
A major advantage of GWR model is the presentation of local variability and local collinearity [77] which are not available in modelling with ANN, SVM and GLM. Local variability for each rural district shows the power of the model in different locations across the study area. Figure 9a demonstrates that GWR model performs more accurately on some rural districts with high local R^{2}. The maximum value is 0.96, and the minimum is 0.16, but the overall R^{2} is 0.85 for the entire study area (Fig. 9a). The other issue is local collinearity, which is unavoidable in modelling and it has adverse effects on the estimation of coefficients. According to many studies, local collinearity of more than 30 indicates decreased reliability of results [78]. GWR shows local collinearity by measuring the condition number for each location. Condition numbers over 30 result in serious concern. Condition number measures how much the output value of the model can change for a small variation in the input of the model. Figure 9b indicates that the obtained condition number for each rural district is less than 20, and the local collinearity is negligible for the prediction of leptospirosis.
Prediction evaluation
GWR, GLM, SVM and ANN models are trained by utilising the data of 2009 and 2010 to predict leptospirosis distribution in 2011. The results are compared with observations of leptospirosis (reported cases) in 2011. Four evaluation criteria, including R^{2}, MAE, MSE and MRE, are employed to assess the results (Table 8). The values of R^{2} are 0.85, 0.78, 0.80 and 0.75 for GWR, GLM, SVM and ANN models, respectively. The values of MSE, MAE and MRE are calculated for GWR (0.050, 0.012 and 0.011), GLM (0.118, 0.052 and 0.017), SVM (0.103, 0.037 and 0.015) and ANN (0.137, 0.063 and 0.018). Needless to say, the lower the values of these criteria, the better the efficiency of the model. Hence, the performance of models in prediction of leptospirosis is GWR > SVM > GLM > ANN. This might be attributed to several reasons: The advantage of GWR as a weighted regression in modelling local variability and spatial heterogeneity, the nature of leptospirosis distribution varying across the study area locally, the superiority of SVM, as a supervised learning approach, in dealing with small classified datasets, the structure of GLM considering a polynomial with constant coefficients throughout the region and the shortcoming of ANN in handling small datasets.
Spatial autocorrelation (Moran’s I) of residuals and significance level
Spatial autocorrelation in the residuals of model verifies weakness in some parts of the model [79]. In this study, weak but meaningful spatial autocorrelation is found in residuals. Environmental parameters model and predict the disease carefully, but the power of model is less in some regions. The capability of Moran’s I is verified in the investigation of residuals [80], so it is used in this study.
The results of Moran’s I are presented in Table 9. A greater convergence of Moran’s I to expected index indicates better performance of clustering [81]. In addition, zscore and pvalue are criteria to determine the fitness of models. The lower value of pvalue and the higher value of zscore elucidate that residuals of models are clustered in some rural districts. The values of Moran’s I, zscore and pvalue are (0.2947, 0.3673 and 0.5406), (6.71, 7.63 and 12.01) and (0.0010, 0.0010 and 0.0012) for GWR, GLM and ANN respectively. Moran’s I of GWR is closer to Expected Index (− 0.0093) and the value of zscore is lower than GLM and ANN. It means GWR presents less deficiency in modelling and predicting leptospirosis. The result of spatial autocorrelation on residual for SVM is not presented in this part because SVM works with the label of classes.
Spatial clusters of GWR, GLM and ANN residuals obtained from Moran’s I approach are presented in Fig. 11. It illustrates the performance of models for prediction in various areas. High–High (HH) shows rural districts surrounded by neighbours with high spatial autocorrelation. Low–High (LH) indicates rural districts that have low spatial autocorrelation of residuals, but their neighbours have high values. Low–Low (LL) presents rural districts surrounded by neighbours with low spatial autocorrelation. Given the high spatial autocorrelation in residuals, HH clusters illustrate the rural districts where the models have lower performance in prediction of leptospirosis.
Conclusion
Leptospirosis is predicted in this study utilizing GWR, SVM, GLM and ANN models. Five input parameters, including temperature, precipitation, humidity, elevation and vegetation are used in this study. Model predictions are investigated statistically and visually to understand the efficiency of used approaches. According to the results, the performance of the models is as follow: GWR > SVM > GLM > ANN. Also, spatial autocorrelation of residuals is used to investigate the deficiency of models. The results prove that GWR presents less deficiency in modelling and predicting leptospirosis. Additionally, based on coefficients of GWR and GLM parameters and sensitivity analysis of SVM and ANN, temperature and humidity have greater effects on the leptospirosis distribution. Moreover, analysis of coefficients shows that higher temperature and humidity coincide with higher disease occurrence in central regions. In contrast, the southeast rural districts have the lowest outbreaks due to lack of related occupations conducive to leptospirosis propagation. In a nutshell, utilizing useful approaches for prediction of leptospirosis can provide health managers and governments with sufficient information to set proper measures for controlling the disease prevalence across the study area.
Many researches including our study are limited based on data and model. As an analytical shortcoming of many disease studies, Modified Areal Unit Problem (MAUP) presents that scale of study is crucial in spatial analysis [82]. In this study, the results of leptospirosis prediction are acceptable at the rural district level, but this disease should be examined in other scales for better understanding the fitness of models. Disease data used in this study are based on the address of patients, whereas the exact locations of the disease occurrence are paddy fields. The paddy fields must be considered as the base level for more accurate analysis, but such data are not available in Iran. More social and epidemiologic parameters should be considered for more accurate prediction.
As future work, the model will be developed by considering socioepidemiologic parameters. Time series models such as Autoregressive Integrated Moving Average (ARIMA) and their comparison with geographically temporal weighted regression is also considered as future work.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Notes
EnzymeLinked ImmunoSorbent Assay
Inverse Distance Weighting
National Aeronautics and Space Administration
Shuttle Radar Topography Mission
Environment for Visualizing Images
The Moderate Resolution Imaging Spectroradiometer
Abbreviations
 AIC:

Akaike Information Criterion
 ANN:

Artificial Neural Network
 ARIMA:

Autoregressive Integrated Moving Average
 BIC:

Bayesian Information Criterion
 CV:

Cross Validation
 ENVI:

Environment for Visualizing Images
 FFNN:

FeedForward Neural Networks
 GIS:

Geographical Information System
 GLM:

Generalized Linear Model
 GWR:

Geographically Weighted Regression
 HC:

Health Centre
 IDW:

Inverse Distance Weighting
 MAE:

Mean Absolute Error
 MAUP:

Modified Areal Unit Problem
 MLP:

Multilayer Perceptron
 MODIS:

The Moderate Resolution Imaging Spectroradiometer
 MRE:

Mean Relative Error
 MSE:

Mean Square Error
 NASA:

National Aeronautics and Space Administration
 NDVI:

Normalised Difference Vegetation Index
 NHCN:

National Health Care Network
 NMHT:

National Ministry of Health and Treatment of Iran
 RBF:

Radial Basis Function
 SRTM:

Shuttle Radar Topography Mission
 SVM:

Support Vector Machine
 VIF:

Variance Inflation Factor
References
Ko AI, Goarant C, Picardeau M. Leptospira: the dawn of the molecular genetics era for an emerging zoonotic pathogen. Nat Rev Microbiol. 2009;7(10):736.
de Vries SG, et al. Travelrelated leptospirosis in the Netherlands 2009–2016: an epidemiological report and case series. Travel Med Infect Dis. 2018;24:4450.
Rafiei A, et al. Review of leptospirosis in Iran. J Mazandaran Univ Med Sci. 2012;22(94):102–10.
Saito M, et al. Leptospiraidonii sp. nov., isolated from environmental water. Int J Syst Evol Microbiol. 2013;63(7):2457–62.
Priya SP, et al. Leptospirosis: molecular trial path and immunopathogenesis correlated with dengue, malaria and mimetic hemorrhagic infections. Acta Trop. 2017;176:206–23.
Thayaparan S, et al. Leptospirosis, an emerging zoonotic disease in Malaysia. Malays J Pathol. 2013;35(2):123–32.
Lau CL, et al. Leptospirosis: an important zoonosis acquired through work, play and travel. Aust J Gen Pract. 2018;47(3):105.
Saini KC, et al. Clinical and etiological profile of fever with thrombocytopenia–a tertiary care hospital based study. J Assoc Physicians India. 2018;66:33.
Zakeri S, et al. Molecular epidemiology of leptospirosis in northern Iran by nested polymerase chain reaction/restriction fragment length polymorphism and sequencing methods. Am J Trop Med Hyg. 2010;82(5):899–903.
Rood EJ, et al. Environmental risk of leptospirosis infections in the Netherlands: spatial modelling of environmental risk factors of leptospirosis in the Netherlands. PLoS One. 2017;12(10):e0186987.
Mayfield HJ, et al. Use of geographically weighted logistic regression to quantify spatial variation in the environmental and sociodemographic drivers of leptospirosis in Fiji: a modelling study. The lancet Planetary health. 2018;2(5):e223–32.
Zhao J, et al. Mapping risk of leptospirosis in China using environmental and socioeconomic data. BMC Infect Dis. 2016;16(1):343.
Ledien J, et al. Assessing the performance of remotelysensed flooding indicators and their potential contribution to early warning for leptospirosis in Cambodia. PLoS One. 2017;12(7):e0181044.
Gutiérrez J, MartínezVega R. Spatiotemporal dynamics of human leptospirosis and its relationship with rainfall anomalies in Colombia. Trans R Soc Trop Med Hyg. 2018;112(3):115–23.
Matsushita N, et al. The nonlinear and lagged shortterm relationship between rainfall and leptospirosis and the intermediate role of floods in the Philippines. PLoS Negl Trop Dis. 2018;12(4):e0006331.
Habus J, et al. New trends in human and animal leptospirosis in Croatia, 2009–2014. Acta Trop. 2017;168:1–8.
Sumi A, et al. Effect of temperature, relative humidity and rainfall on dengue fever and leptospirosis infections in Manila, the Philippines. Epidemiol Infect. 2017;145(1):78–86.
Denipitiya, D., et al., Spatial and seasonal analysis of human leptospirosis in the district of Gampaha, Sri Lanka. 2016.
Pawar SD, et al. Seasonality of leptospirosis and its association with rainfall and humidity in Ratnagiri, Maharashtra. Int J Health Allied Sci. 2018;7(1):37.
Ferreira M, Ferreira M, INFLUENCE OF TOPOGRAPHIC AND HYDROGRAPHIC FACTORS ON THE SPATIAL DISTRIBUTION OF LEPTOSPIROSIS DISEASE IN SÃO PAULO COUNTY. Brazil: an approach using GEOSPATIAL TECHNIQUES and GIS analysis. Germany: International Archives of the Photogrammetry, Remote Sensing & Spatial Information Sciences; 2016. p. 41.
Della Rossa P, et al. Environmental factors and public health policy associated with human and rodent infection by leptospirosis: a land coverbased study in Nan province, Thailand. Epidemiol Infect. 2016;144(7):1550–62.
Tewara MA, et al. Smallarea spatial statistical analysis of malaria clusters and hotspots in Cameroon; 2000–2015. BMC Infect Dis. 2018;18(1):636.
Yu H, et al. Scrub typhus in Jiangsu Province, China: epidemiologic features and spatial risk analysis. BMC Infect Dis. 2018;18(1):372.
Mollalo A, et al. A GISbased artificial neural network model for spatial distribution of tuberculosis across the continental United States. Int J Environ Res Public Health. 2019;16(1):157.
Saeidian B, Mesgari MS, Ghodousi M, Optimum allocation of water to the cultivation farms using Genetic Algorithm. International Archives of the Photogrammetry. Germany: Remote Sensing & Spatial Information Sciences; 2015. p. 40.
Ghaemi Z, Farnaghi M. A varied densitybased clustering approach for event detection from heterogeneous twitter data. ISPRS Int J Geo Inf. 2019;8(2):82.
Saeidian B, et al. Optimized locationallocation of earthquake relief centers using PSO and ACO, complemented by GIS, clustering, and TOPSIS. ISPRS Int J Geo Inf. 2018;7(8):292.
Mollalo A, et al. Geographic information systembased analysis of the spatial and spatiotemporal distribution of zoonotic cutaneous leishmaniasis in Golestan Province, northeast of Iran. Zoonoses Public Health. 2015;62(1):18–28.
HanafiBojd A, et al. Spatial analysis and mapping of malaria risk in an endemic area, south of Iran: a GIS based decision making for planning of control. Acta Trop. 2012;122(1):132–7.
Mollalo A, et al. Machine learning approaches in GISbased ecological modeling of the sand fly Phlebotomus papatasi, a vector of zoonotic cutaneous leishmaniasis in Golestan province, Iran. Acta Trop. 2018;188:187–94.
Ihantamalala FA, et al. Spatial and temporal dynamics of malaria in Madagascar. Malar J. 2018;17(1):58.
Du Z, et al. Extending geographically and temporally weighted regression to account for both spatiotemporal heterogeneity and seasonal variations in coastal seas. Eco Inform. 2018;43:185–99.
Liu Y, et al. Spatial distribution of snow depth based on geographically weighted regression kriging in the Bayanbulak Basin of the Tianshan Mountains, China. J Mt Sci. 2018;15(1):33–45.
Chu HJ, Kong SJ, Chang CH. Spatiotemporal water quality mapping from satellite images using geographically and temporally weighted regression. Int J Appl Earth Obs Geoinf. 2018;65:1–11.
Huang Y, et al. A semiparametric geographically weighted (SGWR) approach for modeling spatial distribution of population. Ecol Indic. 2018;85:1022–9.
LaureanoRosario AE, et al. Application of artificial neural networks for dengue fever outbreak predictions in the northwest coast of Yucatan, Mexico and San Juan, Puerto Rico. Trop Med Infect Dis. 2018;3(1):5.
Dande P, Samant P. Acquaintance to artificial neural networks and use of artificial intelligence as a diagnostic tool for tuberculosis: a review. Tuberc. 2018;108:1–9.
Wang J, et al. A remote sensing data based artificial neural network approach for predicting climatesensitive infectious disease outbreaks: a case study of human brucellosis. Remote Sens. 2017;9(10):1018.
Xu W, Wang Q, Chen R. Spatiotemporal prediction of crop disease severity for agricultural emergency management based on recurrent neural networks. GeoInformatica. 2017;22:1–19.
Nelder JA, Wedderburn RW. Generalized linear models. J R Stat Soc Ser A. 1972;135(3):370–84.
Faraway, J.J., Extending the linear model with R: generalized linear, mixed effects and nonparametric regression models. 2016: Chapman and Hall/CRC.
Vapnik VN. An overview of statistical learning theory. IEEE Trans Neural Netw. 1999;10(5):988–99.
Ch S, et al. A support vector machinefirefly algorithm based forecasting model to determine malaria transmission. Neurocomputing. 2014;129:279–88.
Maghsoudi M, et al. Artificial neural network (ANN) method for modeling of sunset yellow dye adsorption using zinc oxide nanorods loaded on activated carbon: kinetic and isotherm study. Spectrochim Acta A Mol Biomol Spectrosc. 2015;134:1–9.
Brunsdon C, Fotheringham S, Charlton M. Geographically weighted regression. J R Stat Soc Ser A. 1998;47(3):431–43.
Bidanset, P.E. and J.R. Lombard, Optimal kernel and bandwidth specifications for geographically weighted regression. Applied Spatial Modelling and Planning, 2017.
Dong G, Nakaya T, Brunsdon C. Geographically weighted regression models for ordinal categorical response variables: an application to georeferenced life satisfaction data. Comput Environ Urban Syst. 2018;70:35–42.
Mohammadinia A, Alimohammadi A, Saeidian B. Efficiency of geographically weighted regression in modeling human leptospirosis based on environmental factors in Gilan Province, Iran. Geosciences. 2017;7(4):136.
Zhang Z. Artificial neural network, in Multivariate Time Series Analysis in Climate and Environmental Research: Springer; 2018. p. 1–35.
Mayfield H, et al. Use of freely available datasets and machine learning methods in predicting deforestation. Environ Model Softw. 2017;87:17–28.
Walczak S. Artificial neural networks, in Encyclopedia of Information Science and Technology, Fourth Edition. Finland: IGI Global; 2018. p. 120–31.
Da Silva IN, et al. Artificial Neural Networks. Switzerland: Springer; 2017.
Moreira MW, et al. In International Conference on Frontier Computing. Singapore: Springer; 2017.
Naresh Babu, K. And D.R. Edla, New algebraic activation function for multilayered feed forward neural networks. IETE J Res, 2017. 63(1): p. 71–79.
Chatterjee S, et al. Cuckoo search coupled artificial neural network in detection of chronic kidney disease. In: Electronics, Materials Engineering and NanoTechnology (IEMENTech), 2017 1st International Conference on. India: IEEE; 2017.
Reddy VR, Reddy VV, Mohan VCJ. Speed control of induction motor drive using artificial neural networksLevenbergMarquardt Backpropogation algorithm. Int J Appl Eng Res. 2018;13(1):80–5.
Vapnik VN. Statistical learning theory, vol. 2. New York: Wiley; 1998.
Yu H, Kim S. SVM Tutorial—Classification, Regression and Ranking, in Handbook of Natural Computing. Germany: Springer; 2012. p. 479–506.
Burges CJ. A tutorial on support vector machines for pattern recognition. Data Min Knowl Disc. 1998;2(2):121–67.
Yeganeh B, et al. Prediction of CO concentrations based on a hybrid partial Least Square and support vector machine model. Atmos Environ. 2012;55:357–65.
Ghaemi Z, Alimohammadi A, Farnaghi M. LaSVMbased big data learning system for dynamic prediction of air pollution in Tehran. Environ Monit Assess. 2018;190(5):300.
Cawley GC, Talbot NL. Fast exact leaveoneout crossvalidation of sparse leastsquares support vector machines. Neural Netw. 2004;17(10):1467–75.
Nieto PG, et al. A SVMbased regression model to study the air quality at local scale in Oviedo urban area (northern Spain): a case study. Appl Math Comput. 2013;219(17):8923–37.
Ruben GB, et al. Application and sensitivity analysis of artificial neural network for prediction of chemical oxygen demand. Water Resour Manag. 2018;32(1):273–83.
Pianosi F, et al. Sensitivity analysis of environmental models: a systematic review with practical workflow. Environ Model Softw. 2016;79:214–32.
Lowe R, et al. Spatiotemporal modelling of climatesensitive disease risk: towards an early warning system for dengue in Brazil. Comput Geosci. 2011;37(3):371–81.
Chen Z, et al. Efficiency of using spatial analysis for Norway spruce progeny tests in Sweden. Ann For Sci. 2018;75(1):2.
Norouzi J, et al. Predicting renal failure progression in chronic kidney disease using integrated intelligent fuzzy expert system. Comput Math Methods Med. 2016;2016:19.
Jain S, et al. Design of microstrip moisture sensor for determination of moisture content in rice with improved mean relative error. Microw Opt Technol Lett. 2019;61(7):1764–8.
Hox JJ, Moerbeek M, van de Schoot R. Multilevel analysis: Techniques and applications. UK: Routledge; 2017.
Petitpierre B, et al. Selecting predictors to maximize the transferability of species distribution models: lessons from crosscontinental plant invasions. Glob Ecol Biogeogr. 2017;26(3):275–87.
Weiss S, et al. Correlation detection strategies in microbial data sets vary widely in sensitivity and precision. ISME J. 2016;10(7):1669.
Fotheringham AS, Oshan TM. Geographically weighted regression and multicollinearity: dispelling the myth. J Geogr Syst. 2016;18(4):303–29.
Leysen M, et al. Illness perceptions explain the variance in functional disability, but not habitual physical activity, in patients with chronic low Back pain: a crosssectional study. Pain Pract. 2018;18(4):523–31.
Gallagher CV, et al. Development and application of a machine learning supported methodology for measurement and verification (M&V) 2.0. Energ Buildings. 2018;167:8–22.
Guernier V, et al. A systematic review of human and animal leptospirosis in the Pacific Islands reveals pathogen and reservoir diversity. PLoS Negl Trop Dis. 2018;12(5):e0006503.
Siyadatpanah A, et al. Spatial distribution of Giardia lamblia infection among general population in Mazandaran Province, north of Iran. J Parasit Dis. 2018;42(2):171–6.
Nguyen QH, Understanding Factors Affecting the Outbreak of Malaria Using LocallyCompensated Ridge Geographically Weighted Regression: Case Study in DakNong, Vietnam. Advances and Applications in Geospatial Technology and Earth Resources: Proceedings of the International Conference on GeoSpatial Technologies and Earth Resources 2017. Vietnam: Springer; 2017.
Liu S, et al. Predicting the outbreak of hand, foot, and mouth disease in Nanjing, China: a timeseries model based on weather variability. Int J Biometeorol. 2017;62:1–10.
Ali M, et al. Identification of burden hotspots and risk factors for cholera in India: an observational study. PLoS One. 2017;12(8):e0183100.
Lee J, Li S. Extending moran's index for measuring spatiotemporal clustering of geographic events. Geogr Anal. 2017;49(1):36–57.
Nouri H, et al. NDVI, scale invariance and the modifiable areal unit problem: an assessment of vegetation in the Adelaide parklands. Sci Total Environ. 2017;584:11–8.
Acknowledgements
The authors are very grateful to the editor and reviewers’ comments and suggestions, which helped us to revise the manuscript.
Funding
This research is funded by the Centre for Advanced Modelling and Geospatial Information Systems (CAMGIS), University of Technology Sydney (UTS) under grant numbers 321740.2232335, 323930, and 321740.2232357.
Author information
Authors and Affiliations
Contributions
AM collected the data and implemented the GWR and GLM approaches, BS implemented ANN models for prediction, ZG performed SVM method, BP edited, revised, improved the manuscript as expert professor in this field and also arranged the funding for the publication fees. Analysis of data were done by all authors and they read the manuscript, revised and approved the final version.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Additional file 1.
The results of trial and error approach for ANN. The results of trial and error approach for finding the optimal numbers of hidden layers and nodes in layers in final MLP architecture were presented.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Mohammadinia, A., Saeidian, B., Pradhan, B. et al. Prediction mapping of human leptospirosis using ANN, GWR, SVM and GLM approaches. BMC Infect Dis 19, 971 (2019). https://doi.org/10.1186/s1287901945804
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s1287901945804
Keywords
 Leptospirosis
 GIS
 ANN
 GWR
 SVM
 GLM
 Machine learning
 Prediction