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]]
You should seek local statistical help. This list is not the right place to seek what appears to be statistical consulting and, indeed, basic statistical instruction. -- Bert On Fri, May 31, 2013 at 6:27 PM, Luis Fernando Garc?a Hern?ndez <luysgarcia at gmail.com> wrote:> 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]] > > > ______________________________________________ > 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. >-- Bert Gunter Genentech Nonclinical Biostatistics Internal Contact Info: Phone: 467-7374 Website: http://pharmadevelopment.roche.com/index/pdb/pdb-functional-groups/pdb-biostatistics/pdb-ncb-home.htm
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.