search for: glmmmpql

Hi, I am trying to use glmmPQL package for Generalized linear mixed models. This package works by repeated calls to lme. lme uses by default REML method for estimation. Then, does glmmmPQL use REML too? In contrast, how can I change it? I have tried it, writing : method="REML", but the program says: invalid method REML. If somebody can answer me....thanks, Sonja

Dear all, I am experiencing problems with glmmmPQL. I am trying to analyze binomial data with some spatial autocorrelation. Here is my code and some of the outputs > colnames(d.glmm) [1] "BV" "Longitude" "Latitude" "nb_pc_02" "nb_expr_02" [6] "pc_02" "nb_pc_...