data(rotif.env)
# note the 'categories' below are made up, as rotifers are not on yet redlisted
# see Details above for how to get actual values for your species
vulnerability(rotif.env[ , 18], categories = 8)
vulnerability(rotif.env, sp.cols = "Abrigh", categories = 8)
vulnerability(rotif.env, sp.cols = c("Apriod", "Burceo", "Kcochl"), categories = c(8, 16, 2))
# fuzzy vulnerability (from favourability values):
pred <- multGLM(rotif.env, sp.cols = c("Apriod", "Burceo", "Kcochl"), var.cols = 5:17)$predictions
head(pred)
vulnerability(pred, sp.cols = "Apriod_F", categories = 8)
vulnerability(pred, sp.cols = c("Apriod_F", "Burceo_F", "Kcochl_F"), categories = c(8, 16, 2))
Run the code above in your browser using DataLab