glm.negbin = glm.nb(Faramea.occidentalis~Elevation, data=faramea) summary(glm.negbin)