library(mvabund)data(spider)X <- spider$x
abund <- spider$abund
# Example 1: Simple examplemyfamily <-"negative.binomial"# Example 1: Funkier example where Species are assumed to have different distributions# Fit models including all covariates are linear terms, but exclude for bare sandfit0 <- stackedsdm(abund, formula_X =~. -bare.sand, data = X, family = myfamily, ncores=2)fitted(fit0)# Example 2: Funkier example where Species are assumed to have different distributionsabund[,1:3]<-(abund[,1:3]>0)*1# First three columns for presence absencemyfamily <- c(rep(c("binomial"),3), rep(c("negative.binomial"),(ncol(abund)-3)))fit0 <- stackedsdm(abund, formula_X =~ bare.sand, data = X, family = myfamily, ncores=2)fitted(fit0)