data(Y, X, G)
# define abiotic part of the model
env.formula = "~ X_1 + X_2"
# Run the model with bottom-up control using stan_glm as fitting method and no penalisation
# (set iter = 1000 to obtain reliable results)
m = trophicSDM(Y, X, G, env.formula, iter = 50,
family = binomial(link = "logit"), penal = NULL,
mode = "prey", method = "stan_glm")
# Plot just the first three species
# \donttest{
plot(m, species = c("Y1","Y2","Y3"))
# }
# If species = NULL (default), all species are plotted.
Run the code above in your browser using DataLab