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!

All the best!

	[[alternative HTML version deleted]]

