An R model to calculate the environmental optima for a species response curve

This record contains the R code for an R2OpenBugs Bayesian model that fits penalised splines to a species response curve and then estimates the means and 95% credible intervals for the optimum, peak, upper and lower limits, and niche breadth. Two response curves can then be compared and the probability that one curve has an optima greater than the other can then be calculated. Six files are included: Two R2OpenBugs models (one logit transformed to deal with presence-absence data, and one untransformed), the R code for running the models on example data, and 3 files containing the relevant example data (Antarctic mosses). More information can be found in Ashcroft et al. (2016) Ecological Informatics 34: 35–43, which should be cited in any publications using this model.