[R] glmmPQL and random effects

Spencer Graves spencer.graves at pdf.com
Tue Feb 14 02:45:20 CET 2006


	  With repeated measures, you will want to use the "correlation" 
argument in "glmmPQL".  For help with how to do that, you need to know 
that "glmmPQL" calls "lme" repeatedly.  Therefore, if you have any 
trouble using it, I suggest you consult Pinheiro and Bates (2000) 
Mixed-Effects Models in S and S-PLUS (Springer).  Over the past five 
years, I may have consulted this book more frequently than any other, 
second only perhaps to Venables and Ripley (2002) Modern Applied 
Statistics with S, 4th ed. (Springer).  If you don't have these books, I 
suggest you consider purchasing them -- and / or recommending to your 
university library that they obtain copies.  I see on 
"http://www.unileon.es/servicios/unileoncat/acquire_web_spi.htm" that 
they invite "Sugerencias de Adquisitión".

	  If you can't get "glmmPQL" to work with your favorite model, I 
suggest you try passing the problem directly to "lme", possibly with 
random numbers for a response.  After I was confident I understood the 
"lme" syntax, I'd then return to "glmmPQL".

	  If you still have troubles, please submit another post -- after first 
the posting guide! "www.R-project.org/posting-guide.html".  I believe 
the people who more closely follow the procedure in the Posting Guide 
are on average rewarded by a faster rate of new knowledge acquisition. 
I have only anecdotes to support that conjecture, but anecdotes can be 
valuable in the absence of more systematic data collection and analysis.

	  Espero que estos commentos le ayude.

	  Spencer Graves
	
dbahrd at unileon.es wrote:

> 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!
> 
> Regards,
> 
> Hugo
> 
> Hugo Robles
> Department of Animal Biology
> University of León (Spain)
> 
> ______________________________________________
> R-help at 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




More information about the R-help mailing list