Modelling zoonotic diseases in humans: comparison of methods for hantavirus in Sweden
 Caroline B Zeimes^{1}Email author,
 Gert E Olsson^{2},
 Clas Ahlm^{3} and
 Sophie O Vanwambeke^{1}
https://doi.org/10.1186/1476072X1139
© Zeimes et al.; licensee BioMed Central Ltd. 2012
Received: 26 June 2012
Accepted: 10 September 2012
Published: 17 September 2012
Abstract
Because their distribution usually depends on the presence of more than one species, modelling zoonotic diseases in humans differs from modelling individual species distribution even though the data are similar in nature. Three approaches can be used to model spatial distributions recorded by points: based on presence/absence, presence/available or presence data. Here, we compared one or two of several existing methods for each of these approaches.
Human cases of hantavirus infection reported by place of infection between 1991 and 1998 in Sweden were used as a case study. Puumala virus (PUUV), the most common hantavirus in Europe, circulates among bank voles (Myodes glareolus). In northern Sweden, it causes nephropathia epidemica (NE) in humans, a mild form of hemorrhagic fever with renal syndrome.
Logistic binomial regression and boosted regression trees were used to model presence and absence data. Presence and available sites (where the disease may occur) were modelled using crossvalidated logistic regression. Finally, the ecological niche model MaxEnt, based on presenceonly data, was used.
In our study, logistic regression had the best predictive power, followed by boosted regression trees, MaxEnt and crossvalidated logistic regression. It is also the most statistically reliable but requires absence data. The crossvalidated method partly avoids the issue of absence data but requires fastidious calculations. MaxEnt accounts for nonlinear responses but the estimators can be complex. The advantages and disadvantages of each method are reviewed.
Keywords
Introduction
Modelling point records of presence of zoonotic disease
Zoonotic diseases are complex to model because pathogen presence in humans results from the interaction between humans, hosts, and the environment. In this way, species distribution modelling may be used, but interpretation of results may differ.
Many ecological and epidemiological spatial records are points. They relate to locationspecific records of discrete units such as organisms or reported disease cases. A number of models allow investigating and predicting presence of organisms and pathogens based on a set of independent variables. These methods address in various ways the issue of confronting (or not) presences with absences. Recording a presence may be interpreted as a probabilistic function that depends on the abundance of the species/disease and on its detectability [1]. Absences, i.e. places where it is undoubted that the organism/pathogen is not present, may be recorded but are often a set of points randomly chosen through the study area. Absences may be interpreted in three ways [2]:
Environmental absences, related to unfavorable environmental and climatic conditions (not in potential or realized distribution),
Contingent absences, located in favorable areas (within the potential but not in the realized distribution) and,
Methodological absences, caused by a bias in the data collection.
If the ability to detect a species is constant across the study area (and differs from zero), then absences are reliable or associated to habitats where prevalence is low [1]. Absences of a zoonotic disease imply the absence of at least one of five elements: animal donor, vector, animal recipient, pathogen and an external environment allowing pathogen circulation [3]. Human records of zoonotic diseases can be assumed to approximate suitably the underlying zoonotic cycle if cases are well reported and human population distribution is relatively continuous. As animal data is challenging to collect, human data may offer a suitable alternative, especially when large areas are studied.
Models based on point data can be classified into three categories with regard to input data : presence and absence, presence and available, or presenceonly. Here, representative methods of each approach were investigated and compared. First, binomial logistic regression, based on presence and absence data, was modeled. It fits a logistic curve between dependent variable and explanatory variables. Second, boosted regression trees (BRT) were tested. It is a decision tree where predictive performance is improved by boosting [4]. Third, crossvalidated logistic regression (CV method), as introduced by Boyce et al. (2002), was computed. Points usually compared to presences may be better considered as undergoing a different intensity of use rather than being strict absences. The CV method considers available points instead of absence points [5]. Finally, an ecological niche model relying on presenceonly data, specifically the MaxEnt model, was used. MaxEnt was chosen because it is frequently used and well documented [6, 7]. Presenceonly approaches use absences implicitly. The probabilities computed by the four models were mapped. The outputs of each model were compared using AUC and the kappa index. The presence and absence approach and presenceonly approach have been compared [1, 7–11] but here, in addition to comparing the predictive power or the goodness of fit, advantages and disadvantages are reviewed. Focus on modelling a zoonotic disease in humans implies the consideration of the preferences of multiple species.
Case study: human hantavirus infections in northern Sweden
Human hantavirus infections were chosen as a zoonotic disease of public health importance in Europe, and a major rodentborne disease [12]. In Sweden, Puumala hantavirus (PUUV), (Bunyaviridae)[13, 14], is the most prevalent hantavirus and the only pathogenic one [15, 16]. Its host is the bank vole (Myodes glareolus) [17]. In humans, PUUV causes nephropathia epidemica (NE) [18], a mild form of hemorrhagic fever with renal syndrome (HFRS) [19]. Transmission to humans may be direct by biting but is mainly indirect by breathing aerosolised urine and feces of infected voles [20]. Human infection often occurs during the cleaning of closed and unaired buildings or while handling firewood [16, 21, 22]. At room temperature (and colder) and away from UV light, the virus remains infective for at least two weeks [20]. The number of recorded cases of HFRS in Europe (and in Sweden) has increased recently, which may be partly related to increasing surveillance and possibly to climatic factors [23–26]. In Scandinavia, the peak of NE occurs from November to December. Cases are however recorded yearround [27]. In Sweden, 90% of all NE cases notified are reported from the four northernmost counties [16].
Previous studies have showed that NE is linked to host abundance [16, 25, 28–30] and human risk activities (forestry, farming, wood cutting, construction work, camping, cleaning and/or redecorating building with rodents’ access …) [16, 31]. Virus prevalence and transmission depend on local environmental, anthropogenic, genetic, behavioral and/or physiological factors [32]. Here we focused on environmental factors related to bank vole habitat, ex vivo virus survival and human presence, that influence the spatial distribution of disease [33].
Materials and methods
Data sources
During this period and in the region, a total of 1,724 cases of NE were notified, and 1,305 persons (76%) responded to the survey. Of these, 862 were confident about the time and location of PUUV exposure but only 217 could provide information detailed enough to link them to such an exact location as an estate. Data are reported by centroid of the land holding where the infection was acquired. Of the 217 cases recorded, some occurred in the same location. Only one record was kept for each location, leaving 212 presence points. 300 isolated dwellings were selected at random from the Lantmäteriet database (Swedish mapping, cadastral and land registration authority). They were used as absence points in the logistic regression and boosted regression trees and as available points in the CV model (Figure 1).
Independent variables and hypothesized relationships with the abundance of bank voles, the ex vivo virus survival and the human presence
Variable in logistic model, BRT and CV model  Variable in MaxEnt model  Abundance of bank voles  Ex vivo virus survival  Human presence  Source 

Area of forests in a 3km radius around the dwelling (m^{2})  Forests  x  SLU Skogskarta  
*Mean volume of spruce per hectare in a 3km radius around the dwelling (m^{3}/ha)  Volume of spruce  x  SLU Skogskarta  
Mean volume of pines per hectare in a 3km radius around the dwelling (m^{3}/ha)  Volume of pine  x  SLU Skogskarta  
*Maximum distance to forests in a 3km radius around the dwelling (m)  x  SLU Skogskarta  
Number of patches of forests 3km radius  x  SLU Skogskarta  
Mean shape index of forests 3km radius  x  SLU Skogskarta  
Mean contiguity index of forests in a 3km radius  x  SLU Skogskarta  
Mean Euclidian nearestneighbor distance between patches of forests in a 3km radius (m)  x  SLU Skogskarta  
*Area of peat bogs in a 3km radius around the dwelling (m^{2})  Peat bogs  x  SVK  
Mean snow depth between 1991 and 1998 (cm)  Snow depth  x  x  SMHI  
Average duration of the snow when it is present for at least 10 days (days)  Snow period  x  x  SMHI  
Majority of grain size of the soil (1 = coarse, 2 = medium, 3 = fine) in a 3km radius  Soil grain size  x  SGU  
*Elevation (m)  Elevation  x  x  x  Aster GDEM 
*Distance to the sea coast (m)  Distance to the sea  x  x  x  SVK 
*Population density (inhabitant/km^{2})  Population density  x  Gridded population of the world  
Total length of public roads in a 3km radius (m)  Roads  x  SVK  
*Distance to holiday homes (m)  Holiday homes  x  Statistiska Centralbyran  
Total length of the water ways in a 3km radius (m)  Water ways  x  Swedish Places 
Ex vivo virus survival depends on humidity and temperature [20]. Soil grain size was used as a proxy for soil humidity: soils with thinner particles will retain more moisture and allow better virus persistence [38, 39]. Data on soil grain size were extracted from the Geological Survey of Sweden (SGU) and classified into coarse, medium and fine particles. A thick snow cover provides high levels of humidity, cold temperature and protection against UV light therefore contributing to better ex vivo virus survival [20]. Snow also affects abundance of bank voles by providing food and shelter preserved against harsh weather and predators [40–43]. Snow depth and average snow duration (when only present for at least 10 days) were computed from interpolated weather station data of the Swedish Meteorological and Hydrological Institute.
Human density presence follows a double gradient: from South to North and from East to West. Population density (Gridded Population of the World from Center for International Earth Science Information Network (CIESIN)) was used to reflect the spatial distribution of humans. Distance to the sea and elevation follow also this gradient. Elevation data were extracted from the Aster GDEM elevation data (Global Digital Elevation Model, Earth Remote Sensing Data Analysis Center (ERSDAC)). These two variables may act as proxies for climate, soil composition and attractiveness of the landscape. Other variables were also chosen to reflect human presence by attractiveness and accessibility of the landscape: distance to the nearest holiday house, extracted from the Central Statistical Bureau data (Satistiska Centralbyran) and, length of water ways and roads in a three kilometer radius (Lantmäteriet).
Independent variables with nonnormal distribution were logtransformed (volume of spruce, distance to forests, area of peat bogs, elevation, distance to sea coast, population density and distance to holiday homes). For logistic regression, BRT and CV model, some variables were expressed as a value in a radius around the infection place, allowing consideration of the landscape encountered around the place of infection. MaxEnt however only allows spatially continuous variables and cannot integrate these variables in a straightforward fashion. Data layers included in MaxEnt were: forests, volume of spruce, volume of pine, peat bogs, soil grain size, snow depth, snow period, elevation, distance to the sea, population density, roads, holiday houses and water ways.
Models
Presence vs. absence: logistic regression
Estimators are calculated by maximum likelihood to maximize the probability of obtaining the observed sample [44, 45]. The intercept determines the position of the logistic curve on the dependent variable [46]. The coefficient of independent variable is the rate of change of the logit function per unit of change of the variable [44]. This estimator shows how fast the curve will increase or decrease.
First, univariate analyses (only one explanatory variable) were carried out with each independent variable. After, variables were selected for the multiple model using a backward stepwise procedure in R (“stats” package). The Akaike Information Criterions (AIC) was used to select the best model and the Variance Inflation Factor (VIF) (“car” package) was checked to avoid multicollinearity issues. Interactions between variables were tested but none was significant at the level of five percent. As quadratic terms decreased the goodnessoffit of the model, they were not included.
Presence vs. absence: boosted regression trees
BRT combines decision trees with boosting to improve the performance (“gbm” package in R) [4]. In a regression tree, a branch leads to several internal nodes or to a terminal node. The path chosen at each internal node depends on the value of the explanatory variables. At a terminal node, a decision is made on presence or absence. Decision trees are built by recursive binary split: initial trees are enlarged by new binary split made on the previous trees. Boosting allows improving the optimization by adding new trees that reduce the most the loss in predictive performance. The procedure is forward and stage wise: after one step, a new tree is fitted on the residual of the previous tree and the new model, with new residuals, contains the previous and the new trees. BRT also include stochasticity defined by the bag fraction, the percentage of data randomly selected at each step. The default bag fraction is 0.5.
Three parameters must be defined: the learning rate, the tree complexity and the number of trees [4]. The learning rate is the contribution of each tree to the model. A low learning rate, which implies a larger number of trees, is advised. The tree complexity represents the number of nodes in a tree. A higher tree complexity implies thus a lower learning rate. Learning rate and tree complexity are chosen based on a visual analysis of graphs. Graphs represent, for a given tree complexity and at different learning rates, the loss of predictive performance (here, predictive deviance) according to the number of trees. A slower learning rate is generally preferable. The optimal number of trees is found when the predictive deviance is lowest.
The final tree is too big to be graphed but the contribution of each variable can be calculated and the effect of the variables (on the probabilities) graphed. Interactions are automatically modeled because the response of one variable depends of the previous responses of the other variables higher in the tree [4]. The relative strengths of interactions are reported and they can be plotted.
Presence vs. availability: crossvalidated logistic regression
Where w(x_{i}) is the midpoint probability of the cluster i and A(x_{i}) is the area of cluster i (here, the number of observations in cluster i).
New predicted presences were calculated by multiplying U(x_{ i }) by the total number of observations for each class. These predicted presences can be compared with observed presences for each class:
 1.
Spearman coefficient (and χ ^{2} test of goodnessoffit) compared predicted and observed values. A high positive correlation is desired.
 2.
A linear regression of predicted cases (x) on observed cases (y) was modeled:
 a.
R^{2} was used to assess the predictive power.
 b.
The intercept was expected to be zero and the estimated regression coefficient was expected to be around one.
 a.
Presence only: MaxEnt model
MaxEnt is a model that maximizes entropy satisfying any constrain on the unknown distribution [49], and that minimizes the relative entropy between two probability densities defined in the covariate space, one estimated from the presence data and one from the landscape [50]. Unlike the two previous models, dependent variables are continuous spatial data layers with identical extent and resolution. Here, the extent was the four northern Swedish counties where NE is recorded and the resolution was one kilometer.
MaxEnt builds an occurrence model starting from a uniform distribution of probability for each cell of the raster [51]. Then, it improves the model iteratively until the gain becomes saturated. The gain is a likelihood statistic which maximizes the probability of presence according to the background data.
All dependent variables are used simultaneously. Collinear variables are usually not considered a problem because if a variable has a significant impact on probabilities, variables correlated to it will have little impact [50].
MaxEnt allows to account for sampling bias by including an additional layer representing the relative survey effort across the landscape. Population density layer was tested as sampling bias. MaxEnt also provides response curves showing the influence of a variable on the probability of presence. Jackknife analyses were used in order to evaluate the contribution of each variable to the model. Fivefold crossvalidation was also used.
Comparison between models
Due to the limited size of the database, all data points were used for model training. No external validation was carried out, and internal indices were used to compare models.
As the logistic and CV models were built with the same variables and dataset, the Akaike’s Information Criteria (AIC) can be used. AIC is a measure of the goodnessoffit of the model.
To evaluate the predictive power of the four models, the area under the curve (AUC), from a receiver operating characteristics (ROC) analysis, was calculated (“PresenceAbsence” package in R) [52, 53]. The rate of true positives is plotted against the rate of false positives at all thresholds of classification into presence and absence. An AUC equals to 0.5 is a random distribution of predictions and an AUC equals to one, a perfect prediction.
Cohen’s kappa statistic, an index of agreement for positive and negative observations, was also calculated [54]. A kappa above 0.75 indicates an excellent agreement; between 0.4 and 0.75, a fair to good agreement and under 0.4, a poor agreement.
Probability maps were created for each model. MaxEnt provides a continuous probability map. For logistic regression, BRT and CV, predicted points probabilities were interpolated by kriging to obtain continuous maps.
False positives and false negatives were also mapped. Even if the CV model usually does not classify probabilities into presences and absences, a map was made for comparison. The probability threshold was chosen at the level where sensitivity (number of true positives divided by the sum of true positives and false negatives) equals specificity (number of true negatives divided by the sum of true negatives and false positives).
As inputs variables vary between methods and as MaxEnt AUC is calculated over the entire study area (presences and background), while the others only considered the set of points, AUCs and kappas were computed on an identical set of points and variables, in order to make an accurate comparison.
Partial analyses with variables related to bank voles, virus and humans
Partial logistic regressions, BRT and MaxEnt models were fitted using variables related to each element. In this manner, the relative importance of bank voles, virus and humans distributions on human infections distribution may be speculated.
Results
Logistic regression
Several variables were significant in the univariate analyses (p < 0.05):

with a positive sign: logarithm of mean volume of spruce, mean volume of pine, logarithm of distance to forests, logarithm of population density and,

with a negative sign: logarithm of area of bogs, snow depth, snow period, logarithm of elevation, logarithm of distance to sea, logarithm of distance to holiday home.
Models obtained by logistic regression and crossvalidated logistic regression method
Estimate of logistic regression  Meanestimate of CV method  

Intercept  −5.371**  −5.447 
Area of forests  8.048*10^{8}***  8.133^{8} 
Log (distance to forests)  1.665***  1.689 
Contiguity  1.198  0.226 
Snow depth  −0.016  −0.016 
Log (distance to sea)  −0.470**  −0.471 
Log (population density)  0.544*  0.109 
AIC  629.74  792.77 
AUC  0.972  0.721 
Boosted regression tree
A learning rate of 0.01 and a tree complexity of five were chosen based on visual analyses of graphs, giving an optimal number of trees of 350 for minimizing deviance.
Variables with the most important contributions were: area of forests (11.45%), distance to holiday homes (11.01%), distance to the sea (9.82%), elevation (8.81%) and mean volume of spruce (8.53%).
Interaction effects were the most important between the sum of roads and area of forests, the snow period and snow thickness, and the elevation and area of forests.
The AUC of 0.92 and kappa of 0.65 indicated a good model and a good agreement. Predicted probabilities generally increased from West to East, with an area of minimum probability in the center (Figure 2). The highest probabilities were found along the sea coast. No false absences were found and only 16 false positives.
Crossvalidated logistic regression
The estimated regression coefficients are found in Table 2. As these are averages, the significance degree was not known but no coefficient was close to zero. The Spearman correlation between observed and predicted values was significant (0.92; pvalue = 0.0013). The linear regression between observed and predicted values had an adjustedR^{2} of 0.84. Predicted values were slightly lower than observed values and not around one as expected (0.303; pvalue < 0.001).
The AUC (0.72) indicated that the predictive power was satisfactory but the kappa (0.32) indicated a poor agreement. Probability of presence of disease decreased from South to North and from East to West. Most false positives were in the East (Figure 2).
MaxEnt model
Results with and without the population density layer as sampling bias were similar. Jackknife analyses showed that elevation brought the highest gain when used in isolation from other variables. Roads decreased the gain most when omitted from the model.
The probability of presence increased with the volume of spruce and decreased with the volume of pines. Approximately after 60 m^{3}/ha of pine, the probability of presence decreased, indicating that habitat is less favorable to bank voles. The two variables are probably complementary. When there are fewer pines, there are more spruces and viceversa.
The AUC of MaxEnt model was good (0.91) but, when calculated on the same points than in logistic and CV model, it decreased to 0.66. Kappa was only calculated on the points and indicated poor agreement (0.19). The highest predicted probabilities were found near the sea coast, roads and water ways (Figure 2). There were many false positives in the East.
Comparison between models
When accounting for the confidence interval, the estimators of logistic regression and CV method were similar, except for contiguity and population density (Table 2). AUC and AIC were best for logistic regression.
The linear pattern which appeared in MaxEnt was related to including spatially detailed data on roads and water ways (Figure 2). As the other models were based on points and then interpolated, such a linear pattern cannot appear, but may appear if probabilities were computed per pixel.
Based on the AUC, the logistic regression produced the best model. If a logistic regression is built with the same variables as MaxEnt and if the AUC of MaxEnt is calculated only on the original data points, both AUC were equal to 0.66 and ROC curve were similar, indicating similar goodnessoffit.
The thresholds identified for classifying probabilities into presences and absences were 0.44 for logistic model, 0.42 for BRT, 0.17 for CV model and 0.49 for MaxEnt model. Except for BRT, all methods overestimated presence. Many false positives were near the sea coast.
Partial analyses with variables related to bank voles, virus and humans
AUC of partial models based on variables related to each element
Rodents  Virus  Humans  

Logistic regression  0.732  0.695  0.684 
Boosted regression trees  0.886  0.8244  0.801 
MaxEnt  0.891  0.893  0.922 
Discussion
Advantages and disadvantages of logistic regression, boosted regression trees, crossvalidated logistic regression and Maxent model
Advantages  Disadvantages  

Logistic regression  Best goodnessoffit and predictive power  Need of real absence points 
Inclusion of variables reflecting the surrounding environment  
BRT  Account for nonlinearity of biological processes  Need of real absence points 
Modelling of interactions  Impossible to see all three at one time  
Inclusion of variables reflecting the surrounding environment  Difficulty to extrapolate  
CV method  Available sites instead of absence sites  Fastidious calculations 
Inclusion of variables reflecting the surrounding environment  Limited value compared to logistic regression  
Maxent model  Ease of use  Complex estimators, difficulty to extrapolate 
Spatially continuous results  Need of spatially continuous data  
Accounts for nonlinearity of biological processes  Limited by the coarsest resolution and the smallest extent of variables 
Input data
Logistic regression and BRT required absence data. Here, absences were identified from accurate data on dwellings but, these absence points may be unidentified cases or just an absence of human hantavirus transmission over the study period. Random locations would have been less appropriate as they would not consider human distribution. working on point data allowed the implementation of variables which reflect the surrounding environment such as the composition and configuration of the landscape. Zoonotic transmission indeed relates to factors extending beyond the place of record.
The CV method considered availability rather than absence, therefore avoiding the issue of unreported cases or absences related to the stochasticity of zoonotic disease transmission to humans. As data were points, independent variables reflecting the surrounding environment could also be included.
In MaxEnt, only presence records were required. Heavy constraints lied on the dependent variables (continuous raster maps of same resolution and geographical extent). Here, several variables at the landscape scale concerned landscape structure. This could not be operationalized as continuous variables in a comprehensive and straightforward fashion. On one hand, measures concerning the landscape surrounding infection sites could no longer be used. Continuous rasters could be constructed to represent the landscape variables, but loss of information is inevitable. On the other hand, the spatial pattern of the input variables, which played an important role for the final model was preserved. In partial analyses, the best MaxEnt fit used variables related to humans. The human population in Northern Sweden is highly structured along the sea coast and inland along roads and rivers. As the final model is calculated on each pixel, the linear pattern was more evident than in the others models.
Ease of use
Logistic regression is widely available in statistical softwares and is easily implemented. The R package “gbm” used for BRT includes a userfriendly tutorial. The CV required fastidious calculations. MaxEnt, a free software with a graphical interface, is very userfriendly.
Goodnessoffit and predictive power
Generally, the logistic regression gave the best results. It had the best AUC and so, the best predictive power, followed by BRT and MaxEnt and, finally by CV model. It should be noted that BRT results were quite heavily influenced by the bag fraction. When the same independent variables were used in logistic and MaxEnt models, AUC and Kappa were comparable. So, the stepwise procedure and the input variables based on the surrounding environment allowed a better fit and predictive power.
False positives for a zoonotic disease can be interpreted as a poor prediction, a nonreported case, or the presence of the pathogen in the wild but its absence in humans. NE is generally under diagnosed, and many PUUV infected humans may go undetected. Indeed, up to seven in eight PUUV infected humans may go unrecognized with subclinical symptoms or symptoms mistaken [55]. In maps of predicted versus observed (Figure 2), CV and MaxEnt had more incorrect predictions, indicating a poorer prediction comparing to logistic regression and BRT. Models only based on presence were most likely to overestimate presence. However, false positives may give indications on the potential distribution, while the others approximate the realized distribution. The use of different sets of explanatory variables may also contribute to this, but tests using identical sets of predictors confirmed the results. All models overestimated presence near the sea coast, but BRT did the least.
Interpretation
Logistic regression, BRT and CV models had higher flexibility for the inclusion of diverse variables. Variables with more straightforward biological interpretations and/or closer proxies could be added in the model. It could be argued that it is an attractive feature for explanatory models. In our case, landscape structure variables (e.g. relating to forest structure and arrangement with respect to human habitat) could be included. MaxEnt found powerful associations with altitude, a variable of little biological significance that proxies several other biologically relevant variables such as temperature, snow cover or population density. Use of MaxEnt may thus be less recommended to build explanatory models.
A major advantage of MaxEnt was the production of a spatially continuous result, allowing finer detail and more visually pleasing output and avoiding the necessity to interpolate results spatially. However, this may come at the price of many false positive pixels. It may still be useful for identifying further study sites. As the interpolated surfaces of the other methods are also uncertain, major risk areas could be first outlined by MaxEnt, then, in these areas, other models could be fitted to gain more detail.
MaxEnt and BRT facilitated the identification and interpretation of nonlinear responses. Contributions and Jackknife analyses showed immediately interesting variables. Response curves showed the variation of probabilities in relation to the dependent variables. These curves brought a lot of information but could be quite complex to understand and a good understanding of the system was required. Moreover, in BRT, variables could be strongly correlated, meaning that curves were not reliable. Also, results may be difficult to generalize: in MaxEnt, nonlinearity generated complex estimators and for the BRT, it was not possible to see all trees at one time.
Previous comparisons between presence/absence and presenceonly methods have highlighted that logistic regression is more appropriate in some cases and MaxEnt are more appropriate in others. If absence data are available, logistic regression is better than MaxEnt to discriminate sites with high disease risk [56]. Penalized logistic regression, which avoids performance problems caused by overfitting, performs similarly to MaxEnt and has been found better than standard logistic regression [11]. Another study shows that MaxEnt is slightly better within the known distribution but logistic regression predicts better outside the data distribution [10].
Similar modeling approaches have been used for other hantaviruses. A study on Juquitiba hantavirus infections in humans in Brazil identified risk areas using MaxEnt [57]. The authors concluded that human data were limited for modeling the virus in host populations. A study in Argentine used reservoir host data and logistic regression to estimate risk areas for humans [58]. Another study on Andes hantavirus in Argentina comparing MaxEnt and logistic regression using rodent data and human infection data found good predictive powers for both methods in predicting rodent distribution, while MaxEnt performed less well on human data [59]. In partial analyses, the importance of bank vole distribution was highlighted in logistic regression and BRT. These models allowed including landscape structure variables that describe the rodent habitat in more detail. In our study, MaxEnt model indicated the importance of human distribution because of its spatial pattern. Modelling the spatial distribution of human hantavirus infections requires thus both environmental conditions and human variables.
Our four models were based on environmental conditions and tried to define the intersection of the spatial distributions of bank voles, humans and virus. Care is however needed when interpreting results, particularly differences between potential and realized distributions [8]. Even if all favorable conditions are present, the disease/species is not necessarily found. False positive results may be the result of nontransmission of the pathogen to humans even if it circulates in wild hosts. Moreover, as bank voles have a wide ecological niche, models are less accurate [60]. Modelling zoonotic diseases in humans is best done using human case data as host data often represents a broader distribution. Zoonoses involve several species as well as humans and their activities.
Further proposals
An option could be to use first BRT or MaxEnt, in order to delineate areas of high probabilities. Variables could then be sliced according to their response curve into several variables or transformed into categorical variable. This way, nonlinear processes could be considered. If absences are available, they can be added in logistic regressions or, if not, in the CV method. Noncontinuous landscape variables can then be added. The final purpose of the model, explicative or predictive, would also direct the choice, as would data availability and specificity of the system at hand.
Human hantavirus infections in Sweden
Disease cases were found at the intersection of the distribution of bank voles, humans and virus. Many factors must be taken into account. Distance to the sea, which was included in the logistic model, and elevation, which brought the highest gain in the MaxEnt model were proxies for different phenomena. These variables reflected a double gradient also represented by different explanatory variables. A milder climate is found near the coast and in the south, there were more spruces than pines, the soil was moister and human density was higher. Even if the correlation was not always strong, all variables were interconnected.
Variables included in the logistic model and BRT represented bank vole habitat and its connectivity, survival of the virus and human distribution. Distance to forests and contiguity were measures for connectivity of forest. The connectivity index must be taken with caution because it is not necessarily functional [61, 62]. The habitat of bank voles and the viruspreserving snow cover were important. Other models of hantavirus infections around the world show the importance of land cover and climate [57–59]. In China, a MaxEnt study based on infected rodents highlighted the importance of land cover and elevation [63]. In the USA, a logistic regression model based on human infections showed the importance of elevation, climate and ecotone [64].
Conclusion
Zoonoses, included the rodentborne hantavirus, can be modelled with diverse methods. The methods presented here differed in what they permit and offer, each of which may be more important depending on the study objectives. Each method has advantages and disadvantages. A solution could be to combine the different methods
Declarations
Aknowledgements
This study was funded by EU grants FP7261504 EDENext and is catalogued by the EDENext Steering Committee as EDENext 054 (http://www.edenext.eu). The contents of this publication are the sole responsibility of the authors and don’t necessarily reflect the views of the European Commission. The initial work on data collection was supported by grants from the Centre for Environmental Research in Umeå, the County of Västerbotten, the Medical Faculty of Umeå University, and the County Councils of Northern Sweden. The authors thank N. Hartemink for reviewing.
Authors’ Affiliations
References
 Brotons L, Thuiller W, Araujo M, Hirzel A: Presenceabsence versus presenceonly modelling methods for predicting bird habitat suitability. Ecography. 2004, 27: 437448. 10.1111/j.09067590.2004.03764.x.View ArticleGoogle Scholar
 Lobo J, JimenezValverde A, Hortal J: The uncertain nature of absences and their importance in species distribution modelling. Ecography. 2010, 33: 103114. 10.1111/j.16000587.2009.06039.x.View ArticleGoogle Scholar
 Lambin E, Tran A, Vanwambeke S, Linard C, Soti V: Pathogenic landscapes: Interactions between land, people, disease vectors, and their animal hosts. Int J Health Geogr. 2010, 9: 5410.1186/1476072X954.PubMed CentralView ArticlePubMedGoogle Scholar
 Elith J, Leathwick J, Hastie T: A working guide to boosted regression trees. J Anim Ecol. 2008, 77: 802813. 10.1111/j.13652656.2008.01390.x.View ArticlePubMedGoogle Scholar
 Boyce M, Vernier P, Nielsen S, Schmiegelow F: Evaluating resource selection functions. Ecol Model. 2002, 157: 281300. 10.1016/S03043800(02)002004.View ArticleGoogle Scholar
 Hernandez P, Graham C, Master L, Albert D: The effect of sample size and species characteristics on performance of different species distribution modeling methods. Ecography. 2006, 29: 773785. 10.1111/j.09067590.2006.04700.x.View ArticleGoogle Scholar
 Elith J, Graham H, Anderson P, Dudík M, Ferrier S, Guisan A, Hijmans J, Huettmann F, Leathwick R, Lehmann A: Novel methods improve prediction of species' distributions from occurrence data. Ecography. 2006, 29: 129151. 10.1111/j.2006.09067590.04596.x.View ArticleGoogle Scholar
 JimenezValverde A, Lobo J, Hortal J: Not as good as they seem: the importance of concepts in species distribution modelling. Divers Distrib. 2008, 14: 885890. 10.1111/j.14724642.2008.00496.x.View ArticleGoogle Scholar
 Fielding A, Bell J: A review of methods for the assessment of prediction errors in conservation presence/absence models. Environ Conserv. 1997, 24: 3849. 10.1017/S0376892997000088.View ArticleGoogle Scholar
 Cleve C, Perrine J, Holzman B, Hines E: Addressing biased occurrence data in predicting potential Sierra Nevada red fox habitat for survey prioritization. Endangered Species Research. 2011, 14: 179191. 10.3354/esr00348.View ArticleGoogle Scholar
 Gastón A, GarcíaViñas J: Modelling species distributions with penalised logistic regressions: A comparison with maximum entropy models. Ecol Model. 2011, 222: 20372041. 10.1016/j.ecolmodel.2011.04.015.View ArticleGoogle Scholar
 Vaheri A, Henttonen H, Voutilainen L, Mustonen J, Sironen T, Vapalahti O: Hantavirus infections in Europe and their impact on public health. Reviews in Medical Virology. 2012, doi: 10.1002/rmv.1722.Google Scholar
 Bishop D, Calisher C, Casals J, Chumakov M, Gaidamovich S, Hannoun C, Lvov D, Marshall I, Okerblom N, Pettersson R: BUNYAVIRIDAE. Intervirology. 1980, 14: 125143. 10.1159/000149174.View ArticlePubMedGoogle Scholar
 Hart C, Bennett M: Hantavirus infections: epidemiology and pathogenesis. Microbes Infect. 1999, 1: 12291237. 10.1016/S12864579(99)002385.View ArticlePubMedGoogle Scholar
 Clement J, Lameire N, Keyaerts E, Maes P, Van Ranst M: Hantavirus infections in Europe. Lancet Infect Dis. 2003, 3: 752753. 10.1016/S14733099(03)008272.View ArticlePubMedGoogle Scholar
 Olsson G, Dalerum F, Hornfeldt B, Elgh F, Palo T, Juto P, Ahlm C: Human hantavirus infections, Sweden. Emerg Infect Dis. 2003, 9: 13951401. 10.3201/eid0911.030275.PubMed CentralView ArticlePubMedGoogle Scholar
 Clement J, Heyman P, McKenna P, Colson P, AvsicZupanc T: The hantaviruses of Europe: From the bedside to the bench. Emerg Infect Dis. 1997, 3: 205211. 10.3201/eid0302.970218.PubMed CentralView ArticlePubMedGoogle Scholar
 Lahdevirta J, Savola J, Brummerkorvenkontio M, Berndt R, Illikainen R, Vaheri A: Clinical and serological diagnosis of nephropathia epidemica, the mild type of hemorrhagicfever with renal syndrome. J Infect. 1984, 9: 230238. 10.1016/S01634453(84)90464X.View ArticlePubMedGoogle Scholar
 Leduc J: Epidemiology of hemorrhagicfever viruses. Reviews of Infectious Diseases. 1989, 11: S730S735. 10.1093/clinids/11.Supplement_4.S730.View ArticlePubMedGoogle Scholar
 Kallio E, Klingstrom J, Gustafsson E, Manni T, Vaheri A, Henttonen H, Vapalahti O, Lundkvist A: Prolonged survival of Puumala hantavirus outside the host: evidence for indirect transmission via the environment. J Gen Virol. 2006, 87: 21272134. 10.1099/vir.0.816430.View ArticlePubMedGoogle Scholar
 Dearing M, Dizney L: Ecology of hantavirus in a changing world. Ann Ny Acad Sci. 2010, 1195: 99112. 10.1111/j.17496632.2010.05452.x.View ArticlePubMedGoogle Scholar
 Schmaljohn C, Hasty S, Dalrymple J, Leduc J, Lee H, Vonbonsdorff C, Brummerkorvenkontio M, Vaheri A, Tsai T, Regnery H: Antigenic and genetic properties of viruses linked to hemorrhagicfever with renal syndrome. Science. 1985, 227: 10411044. 10.1126/science.2858126.View ArticlePubMedGoogle Scholar
 Heyman P, Ceianu C, Christova I, Tordo N, Beersma M, Alves M, Lundkvist A, Hukic M, Papa A, Tenorio A: A fiveyear perspective on the situation of haemorrhagic fever with renal syndrome and status of the hantavirus reservoirs in Europe, 2005–2010. Eurosurveillance. 2011, 16: 1522.Google Scholar
 Pettersson L, Boman J, Juto P, Evander M, Ahlm C: Outbreak of Puumala virus infection, Sweden. Emerg Infect Dis. 2008, 14: 808810. 10.3201/eid1405.071124.PubMed CentralView ArticlePubMedGoogle Scholar
 Olsson G, Hjertqvist M, Lundkvist A, Hornfeldt B: Predicting high risk for human hantavirus infections, Sweden. Emerg Infect Dis. 2009, 15: 104106. 10.3201/eid1501.080502.PubMed CentralView ArticlePubMedGoogle Scholar
 Olsson G, Leirs H, Henttonen H: Hantaviruses and their hosts in Europe: reservoirs here and there, but not everywhere?. VectorBorne Zoonot. 2010, 10: 549561. 10.1089/vbz.2009.0138.View ArticleGoogle Scholar
 BrummerKorvenkontio M, Vapalahti O, Henttonen H, Koskela P, Kuusisto P, Vaheri A: Epidemiological study of nephropathia epidemica in Finland 1989–96. Scand J Infect Dis. 1999, 31: 427435. 10.1080/00365549950163941.View ArticlePubMedGoogle Scholar
 Heyman P, Vervoort T, Escutenaire S, Degrave E, Konings J, Vandenvelde C, Verhagen R: Incidence of hantavirus infections in Belgium. Virus Res. 2001, 77: 7180. 10.1016/S01681702(01)002672.View ArticlePubMedGoogle Scholar
 Mills J, Childs J: Ecologic studies of rodent reservoirs: Their relevance for human health. Emerg Infect Dis. 1998, 4: 529537. 10.3201/eid0404.980403.PubMed CentralView ArticlePubMedGoogle Scholar
 Niklasson B, Hornfeldt B, Lundkvist A, Bjorsten S, Leduc J: Temporal dynamics of puumala virusantibody prevalence in voles and of nephropathiaepidemica incidence in humans. Am J Trop Med Hyg. 1995, 53: 134140.PubMedGoogle Scholar
 Piechotowski I, Brockmann S, Schwarz C, Winter C, Ranft U, Pfaff G: Emergence of hantavirus in South Germany: rodents, climate and human infections. Parasitol Res. 2008, 103: S131S137. 10.1007/s0043600810558.View ArticlePubMedGoogle Scholar
 Mills JN: Regulation of RodentBorne viruses in the natural host: implications for human disease. Infectious Diseases from Nature: Mechanisms of Viral Emergence and Persistence. Edited by: Peters CJ, Calisher CH. 2005, Springer Vienna, Vienna, 4557.View ArticleGoogle Scholar
 Ostfeld R, Glass G, Keesing F: Spatial epidemiology: an emerging (or reemerging) discipline. Trends Ecol Evol. 2005, 20: 328336. 10.1016/j.tree.2005.03.009.View ArticlePubMedGoogle Scholar
 Vapalahti O, Mustonen J, Lundkvist A, Henttonen H, Plyusnin A, Vaheri A: Hantavirus infections in Europe. Lancet Infect Dis. 2003, 3: 753754. 10.1016/S14733099(03)008284.View ArticleGoogle Scholar
 Olsson G, White N, Hjalten J, Ahlm C: Habitat factors associated with bank voles (Clethrionomys glareolus) and concomitant hantavirus in northern Sweden. VectorBorne Zoonot. 2005, 5: 315323. 10.1089/vbz.2005.5.315.View ArticleGoogle Scholar
 Kozakiewicz M, Van Apeldoorn R, Bergers P, Gortat T, Kozakiewicz A: Landscape approach to bank vole ecology. Polish Journal of Ecology. 2000, 48: 149161.Google Scholar
 Vanapeldoorn R, Oostenbrink W, Vanwinden A, Vanderzee F: Effects of habitat fragmentation on the bank vole, clethrionomysglareolus, in an agricultural landscape. Oikos. 1992, 65: 265274. 10.2307/3545018.View ArticleGoogle Scholar
 Linard C, Tersago K, Leirs H, Lambin E: Environmental conditions and Puumala virus transmission in Belgium. Int J Health Geogr. 2007, 6: 5510.1186/1476072X655.PubMed CentralView ArticlePubMedGoogle Scholar
 Sauvage F, Langlais M, Yoccoz N, Pontier D: Modelling hantavirus in fluctuating populations of bank voles: the role of indirect transmission on virus persistence. J Anim Ecol. 2003, 72: 113. 10.1046/j.13652656.2003.00675.x.View ArticleGoogle Scholar
 Hansson L, Henttonen H: Gradients in density variations of small rodents  the importance of latitude and snow cover. Oecologia. 1985, 67: 394402. 10.1007/BF00384946.View ArticleGoogle Scholar
 Hansson L, Henttonen H: Rodent dynamics as community processes. Trends Ecol Evol. 1988, 3: 195200. 10.1016/01695347(88)900067.View ArticlePubMedGoogle Scholar
 Hanski I, Hansson L, Henttonen H: Specialist predators, generalist predators, and the microtine rodent cycle. J Anim Ecol. 1991, 60: 353367. 10.2307/5465.View ArticleGoogle Scholar
 Hanski I, Henttonen H, Korpimaki E, Oksanen L, Turchin P: Smallrodent dynamics and predation. Ecology. 2001, 82: 15051520. 10.1890/00129658(2001)082[1505:SRDAP]2.0.CO;2.View ArticleGoogle Scholar
 Hosmer D, Lemeshow S: Applied logistic regression. 1989, Wiley, New YorkGoogle Scholar
 McCullagh P, Nelder JA: Generalized linear models. 1989, Chapman and Hall, London, New York, 2View ArticleGoogle Scholar
 Rogers D: Models for vectors and vectorborne diseases. Advances in Parasitology, Vol 62: Global Mapping of Infectious Diseases: Methods, Examples and Emerging Applications. Volume 62. Edited by: Hay SI, Graham A, Rogers DJ. 2006, Elsevier Academic Press Inc, San Diego, 135.View ArticleGoogle Scholar
 Johnson C, Nielsen S, Merrill E, McDonald T, Boyce M: Resource selection functions based on useavailability data: Theoretical motivation and evaluation methods. J Wildl Manage. 2006, 70: 347357. 10.2193/0022541X(2006)70[347:RSFBOU]2.0.CO;2.View ArticleGoogle Scholar
 Wiens T, Dale B, Boyce M, Kershaw G: Three way kfold crossvalidation of resource selection functions. Ecol Model. 2008, 212: 244255. 10.1016/j.ecolmodel.2007.10.005.View ArticleGoogle Scholar
 Phillips S, Anderson R, Schapire R: Maximum entropy modeling of species geographic distributions. Ecol Model. 2006, 190: 231259. 10.1016/j.ecolmodel.2005.03.026.View ArticleGoogle Scholar
 Elith J, Phillips S, Hastie T, Dudik M, Chee Y, Yates C: A statistical explanation of MaxEnt for ecologists. Divers Distrib. 2011, 17: 4357. 10.1111/j.14724642.2010.00725.x.View ArticleGoogle Scholar
 Phillips S: A brief tutorial on Maxent. Lessons in Conservation. 2012, 3: 107135.Google Scholar
 Pearce J, Ferrier S: Evaluating the predictive performance of habitat models developed using logistic regression. Ecol Model. 2000, 133: 225245. 10.1016/S03043800(00)003227.View ArticleGoogle Scholar
 Fawcett T: An introduction to ROC analysis. Pattern Recognit Lett. 2006, 27: 861874. 10.1016/j.patrec.2005.10.010.View ArticleGoogle Scholar
 Kirkwood B, Sterne J: Essential Medical Statistics. 2003, Blackwell Scientific, Oxford, 2Google Scholar
 Ahlm C, Linderholm M, Juto P, Stegmayr B, Settergren B: Prevalence of serum IgG antibodies to Puumala virus (haemorrhagic fever with renal syndrome) in Northern Sweden. Epidemiol Infect. 1994, 113: 129136. 10.1017/S0950268800051542.PubMed CentralView ArticlePubMedGoogle Scholar
 La Manna L, Matteucci S, Kitzberger T: Modelling phytophthora disease risk in austrocedrus chilensis forests of patagonia. European Journal of Forest Research. 2012, 131: 323337. 10.1007/s1034201105037.View ArticleGoogle Scholar
 Donalisio M, Peterson A: Environmental factors affecting transmission risk for hantaviruses in forested portions of southern Brazil. Acta Tropica. 2011, 119: 125130. 10.1016/j.actatropica.2011.04.019.View ArticlePubMedGoogle Scholar
 Carbajo A, Pardiñas U: Spatial distribution model of a hantavirus reservoir, the longtailed colilargo (Oligoryzomys longicaudatus), in Argentina. Journal of Mammalogy. 2007, 88: 15551568. 10.1644/06MAMMA183R1.1.View ArticleGoogle Scholar
 Andreo V, Glass G, Shields T, Provensal C, Polop J: Modeling potential distribution of oligoryzomys longicaudatus, the andes virus (Genus: Hantavirus) reservoir, in Argentina. EcoHealth. 2011, 8: 332348. 10.1007/s1039301107195.View ArticlePubMedGoogle Scholar
 Tsoar A, Allouche O, Steinitz O, Rotem D, Kadmon R: A comparative evaluation of presenceonly methods for modelling species distribution. Divers Distrib. 2007, 13: 397405. 10.1111/j.14724642.2007.00346.x.View ArticleGoogle Scholar
 Ewers R, Didham R: Confounding factors in the detection of species responses to habitat fragmentation. Biological Reviews. 2006, 81: 117142.View ArticlePubMedGoogle Scholar
 Chetkiewicz CL, St Clair C, Boyce M: Corridors for conservation: integrating pattern and process. Annual Review of Ecology, Evolution and Systematics. 2006, 37: 317342. 10.1146/annurev.ecolsys.37.091305.110050.View ArticleGoogle Scholar
 Wei L, Qian Q, Wang Z, Glass G, Song S, Zhang W, Li XJ, Yang H, Wang X, Fang L, Cao W: Using geographic information systembased ecologic niche models to forecast the risk of hantavirus infection in Shandong Province, China. Am J Trop Med Hyg. 2011, 84: 497503. 10.4269/ajtmh.2011.100314.PubMed CentralView ArticlePubMedGoogle Scholar
 Eisen R, Glass G, Eisen L, Cheek J, Enscore R, Ettestad P, Gage K: A spatial model of shared risk for plague and hantavirus pulmonary syndrome in the southwestern United States. Am J Trop Med Hyg. 2007, 77: 9991004.PubMedGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.