[R] GLMM in R
nalimilan at club.fr
Sat Jun 1 11:54:48 CEST 2013
Le vendredi 31 mai 2013 à 22:27 -0300, Luis Fernando García Hernández a
> Dear Friends,
> I am new on R so I ask you to excuse me if this question sounds fool. I
> want to see if there is a significativa relationship between the mating
> (response variable) and several explanatory variables such as individual
> number (categorical), leg movemente (continous) and the reuse of
> individuals (categorical). My data looks like this
> 1 means mating happened, 0 means did not occur.
> I am trying to organize the data to apply a GLMM to the data, but have not
> been able to do it.I followed the model proposed by CRawley for binary
> response with pseudorreplication but does not wok. The script goes like this
> model2<-lmer(Mating~Replication+(1 ID),family=binomial,method=PQL)
> Can somebody tell me what part is wrong and how could I fix it? If you know
> a better method, will be really welcomed!
You do not tell us what you mean by "wrong", but at least in the form
above it will not work. Try with (1|ID), and use "REML" and "nAGQ"
arguments instead of the deprecated "method".
> All the best!
> [[alternative HTML version deleted]]
> R-help at r-project.org mailing list
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
More information about the R-help