Joint species distribution modelling with the r-package Hmsc
Permalink
Citation
Tikhonov , G , Opedal , O H , Abrego , N , Lehikoinen , A , de Jonge , M M J , Oksanen , J & Ovaskainen , O 2020 , ' Joint species distribution modelling with the r-package Hmsc ' , Methods in Ecology and Evolution , vol. 11 , no. 3 , pp. 442-447 . https://doi.org/10.1111/2041-210X.13345
| Title: | Joint species distribution modelling with the r-package Hmsc |
| Author: | Tikhonov, Gleb; Opedal, Oystein H.; Abrego, Nerea; Lehikoinen, Aleksi; de Jonge, Melinda M. J.; Oksanen, Jari; Ovaskainen, Otso |
| Contributor organization: | Research Centre for Ecological Change Organismal and Evolutionary Biology Research Programme Department of Agricultural Sciences Plant Production Sciences Spatial Foodweb Ecology Group Helsinki Institute of Sustainability Science (HELSUS) Finnish Museum of Natural History Otso Ovaskainen / Principal Investigator Plant Adaptation and Conservation |
| Date: | 2020-03-01 |
| Language: | eng |
| Number of pages: | 6 |
| Belongs to series: | Methods in Ecology and Evolution |
| ISSN: | 2041-210X |
| DOI: | https://doi.org/10.1111/2041-210X.13345 |
| URI: | http://hdl.handle.net/10138/313907 |
| Abstract: | Joint Species Distribution Modelling (JSDM) is becoming an increasingly popular statistical method for analysing data in community ecology. Hierarchical Modelling of Species Communities (HMSC) is a general and flexible framework for fitting JSDMs. HMSC allows the integration of community ecology data with data on environmental covariates, species traits, phylogenetic relationships and the spatio-temporal context of the study, providing predictive insights into community assembly processes from non-manipulative observational data of species communities. The full range of functionality of HMSC has remained restricted to Matlab users only. To make HMSC accessible to the wider community of ecologists, we introduce Hmsc 3.0, a user-friendly r implementation. We illustrate the use of the package by applying Hmsc 3.0 to a range of case studies on real and simulated data. The real data consist of bird counts in a spatio-temporally structured dataset, environmental covariates, species traits and phylogenetic relationships. Vignettes on simulated data involve single-species models, models of small communities, models of large species communities and models for large spatial data. We demonstrate the estimation of species responses to environmental covariates and how these depend on species traits, as well as the estimation of residual species associations. We demonstrate how to construct and fit models with different types of random effects, how to examine MCMC convergence, how to examine the explanatory and predictive powers of the models, how to assess parameter estimates and how to make predictions. We further demonstrate how Hmsc 3.0 can be applied to normally distributed data, count data and presence-absence data. The package, along with the extended vignettes, makes JSDM fitting and post-processing easily accessible to ecologists familiar with r. |
| Subject: |
community ecology
community modelling community similarity hierarchical modelling of species communities joint species distribution modelling multivariate data species distribution modelling FRAMEWORK 1181 Ecology, evolutionary biology |
| Peer reviewed: | Yes |
| Rights: | cc_by |
| Usage restriction: | openAccess |
| Self-archived version: | publishedVersion |
Files in this item
Total number of downloads: Loading...
| Files | Size | Format | View |
|---|---|---|---|
| 2041_210X.13345.pdf | 1.948Mb |
View/ |
This item appears in the following Collection(s)
-
Articles from TUHAT CRIS [44950]
Related items
Showing items related by title, author, creator and subject.
-
(2019)A large array of species distribution model (SDM) approaches has been developed for explaining and predicting the occurrences of individual species or species assemblages. Given the wealth of existing models, it is unclear which models perform best for interpolation or extrapolation of existing data sets, particularly when one is concerned with species assemblages. We compared the predictive performance of 33 variants of 15 widely applied and recently emerged SDMs in the context of multispecies data, including both joint SDMs that model multiple species together, and stacked SDMs that model each species individually combining the predictions afterward. We offer a comprehensive evaluation of these SDM approaches by examining their performance in predicting withheld empirical validation data of different sizes representing five different taxonomic groups, and for prediction tasks related to both interpolation and extrapolation. We measure predictive performance by 12 measures of accuracy, discrimination power, calibration, and precision of predictions, for the biological levels of species occurrence, species richness, and community composition. Our results show large variation among the models in their predictive performance, especially for communities comprising many species that are rare. The results do not reveal any major trade-offs among measures of model performance; the same models performed generally well in terms of accuracy, discrimination, and calibration, and for the biological levels of individual species, species richness, and community composition. In contrast, the models that gave the most precise predictions were not well calibrated, suggesting that poorly performing models can make overconfident predictions. However, none of the models performed well for all prediction tasks. As a general strategy, we therefore propose that researchers fit a small set of models showing complementary performance, and then apply a cross-validation procedure involving separate data to establish which of these models performs best for the goal of the study.
-
(2017)In this article, we use our results from Flores and Kleban (2015 Commun. Math. Phys. 333 389-434, 2015 Commun. Math. Phys. 333 435-81, 2015 Commun. Math. Phys. 333 597-667, 2015 Commun. Math. Phys. 333 669715) to generalize known formulas for crossing probabilities. Prior crossing results date back to Cardy's prediction of a formula for the probability that a percolation cluster in two dimensions connects the left and right sides of a rectangle at the percolation critical point in the continuum limit (Cardy 1992 J. Phys. A: Math. Gen. 25 L201-6). Here, we predict a new formula for crossing probabilities of a continuum limit loop-gas model on a planar lattice inside a 2N-sided polygon. In this model, boundary loops exit and then re-enter the polygon through its vertices, with exactly one loop passing once through each vertex, and these loops join the vertices pairwise in some specified connectivity through the polygon's exterior. The boundary loops also connect the vertices through the interior, which we regard as a crossing event. For particular values of the loop fugacity, this formula specializes to FK cluster (resp. spin cluster) crossing probabilities of a critical Q-state random cluster (resp. Potts) model on a lattice inside the polygon in the continuum limit. This includes critical percolation as the Q = 1 random cluster model. These latter crossing probabilities are conditioned on a particular side-alternating free/fixed (resp. fluctuating/fixed) boundary condition on the polygon's perimeter, related to how the boundary loops join the polygon's vertices pairwise through the polygon's exterior in the associated loop-gas model. For Q is an element of{2, 3, 4}, we compare our predictions of these random cluster (resp. Potts) model crossing probabilities in a rectangle (N = 2) and in a hexagon (N = 3) with high-precision computer simulation measurements. We find that the measurements agree with our predictions very well for Q is an element of{2, 3} and reasonably well if Q = 4.
-
(Elsevier, 2019)Science of The Total Environment 658 (2019), 1278-1292In Finland, a recent ecological classification of surface waters showed that the rivers and coastal waters need attention to improve their ecological state. We combined eco-hydrological and empirical models to study chlorophyll-a concentration as an indicator of eutrophication in a small agricultural river. We used a modified story-and-simulation method to build three storylines for possible changes in future land use due to climate change and political change. The main objective in the first storyline is to stimulate economic activity but also to promote the sustainable and efficient use of resources. The second storyline is based on the high awareness but poor regulation of environmental protection, and the third is to survive as individual countries instead of being part of a unified Europe. We assumed trade of agricultural products to increase to countries outside Europe. We found that chlorophyll-a concentration in the river depended on total phosphorus concentration. In addition, there was a positive synergistic interaction between total phosphorus and water temperature. In future storylines, chlorophyll-a concentration increased due to land use and climate change. Climate change mainly had an indirect influence via increasing nutrient losses from intensified agriculture. We found that well-designed agri-environmental measures had the potential to decrease nutrient loading from fields, as long as the predicted increase in temperature remained under 2 °C. However, we were not able to achieve the nutrient reduction stated in current water protection targets. In addition, the ecological status of the river deteriorated. The influence of temperature on chlorophyll-a growth indicates that novel measures for shading rivers to decrease water temperature may be needed in the future.
