search for: hydrovar2

...as a random effect (fAVGNTS). I have also conducted a join count test which suggests that there is some spatial autocorrelation. Consequently I have used the following code: library(MASS) attach(Birds) Birds$x <- Birds$LONGITUDE Birds$y <- Birds$LATITUDE M <- glmmPQL(PRESENCE~ HYDROVAR1 + HYDROVAR2, random= ~ 1|fAVGNTS, correlation = corExp(form = ~ x + y), family = binomial(link = "logit"), data = Birds) The model seems to run fine. However, when I compare the results of this model and the residual spread against the same model but without the correlation function, there is abso...