[R] GLMM in R

[Milan Bouchet-Valat]

Le vendredi 31 mai 2013 à 22:27 -0300, Luis Fernando García Hernández a
écrit :
> 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
> 
> MatingIDReplicationShaking1M1R1100M1R2140M2R1150M2R2121M3R1140M3R2171M4R1190
> M4R2221M5R1180M5R2161M6R117
> 
> 
> 
> 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".


Regards

> All the best!
> 
> 	[[alternative HTML version deleted]]
> 
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.