I have been using lme to model data with multiple nested random effects and continuous response variables however I also have data with a binary response variable, binomial errors and multiple levels of nesting of random effects (e.g. site/block/quadrat/year), is there a package available which will do this? Jim Lindsey's package "repeated" appears to be only able to cope with 2 levels of nesting. Also I have been looking at the GLMMGibbs package but cannot locate the document "GLMMGibbs: An R Package for Estimating Bayesian Generalised Linear Mixed Models by Gibbs Sampling" which is supposed to be distributed with the package. Does anyone have this or can they tell me where to find it? Thanks, Yvonne Buckley ******************************************************************** Yvonne Buckley Department of Biology Imperial College at Silwood Park Ascot, Berkshire SL5 7PY tel: 0207 5942 364 (office) tel: 01344 891 539 (home) fax: 0207 5942 339 email: y.buckley at ic.ac.uk -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._