[R] glmmPQL and random effects

From: <dbahrd_at_unileon.es>
Date: Sat 11 Feb 2006 - 04:53:01 EST

Hello R users,

I am trying to run a model with a binary response variable (nesting success: 0 failure, 1 success) and 8 fixed terms. Nesting success was examined in 72 cases in 34 territories (TER) during a 6 study years. Territories are nested within 14 patches (PATCH). I want to run a model taking into account these nested factors and repeated observation. To do this, I assume that the best option is to use glmmPQL from MASS package. Am I wrong?

In glmmPQL, I have included the random terms as follow: random=~1|YEAR/PATCH/TER, but I am unclear if this syntax is right for this case (?).

I would greatly appreciate any help!



Hugo Robles
Department of Animal Biology
University of León (Spain)

R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Sat Feb 11 04:57:58 2006

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:42:27 EST