library(mvabund) ## Load a dataset from the mvabund package
data(spider)
y <- spider$abun
spider.fit.nb <- boral(y, family = "negative.binomial", num.lv = 2,
site.eff = TRUE, save.model = FALSE, calc.ics = TRUE)
summary(spider.fit.nb)
Run the code above in your browser using DataLab