[R] Repeated measure Generalized Linear Mixed Model (glmm) ?

From: Jérôme Lemaître <je_lemaitre_at_hotmail.com>
Date: Tue 04 Oct 2005 - 21:16:33 EST


Dear all,

I want to know the probability of dying in a trap as a function of habitat variables for a metapopulation of voles sampled in 44 stations. My dataset is as follow:

Station		tag number	dead	habitat1	habitat2
1		1		yes	20		26
1		2		no	20		26
2		3		no	15		16

..

As far as I know, I should use a mixed-model as: glmmPQL<-(fixed = dead~habitat1+habitat2, random = ~1|station, family = binomial) in the MASS library.
However, some vole individuals were recaptured (but it is not a study designed for capture-mark-recapture!!!) therefore, I have for example, the 4th vole:

Station		tag number	dead	habitat1	habitat2
3		4		no	10		12
3		4		yes	10		12

Someone suggested me to use glmm with repeated measures with the number of tag as the variable to repeat. Apparently, it is possible to do this in SAS.

QUESTIONS :
1) Is this also possible to do this in R, and if yes, could you please tell me how?
2) Is this a valid way to analyse this data, and if not, could someone please put me in the right way?

Thank you very very much to all in advance

Jérôme Lemaître

Étudiant au doctorat
Université Laval
Québec, QC G1K 7P4



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 Tue Oct 04 21:27:05 2005

This archive was generated by hypermail 2.1.8 : Sun 23 Oct 2005 - 18:17:30 EST