Hi folks, I am looking for the package that will allow me to do a generalized (poisson) linear mixed model with spatial correlation structure. If gls() in nlme does this, I don't understand how to implement different families. If glmmPQL() in MASS does this, I don't understand what correlation models it accepts. It does not appear to accept the same models as gls(), but I may be doing something wrong. Is there a mixed model package or function I am overlooking? Thanks, Hank Martin Henry H. Stevens 338 Pearson Hall Department of Botany Ecology Graduate Program Miami University Oxford, OH 45056 USA tel: 513-529-4206 fax: 513-529-4243 http://www.cas.muohio.edu/botany/bot/henry.html http://www.units.muohio.edu/ecology/ http://www.cas.muohio.edu/botany/