[R] GLMM in R

Bert Gunter gunter.berton at gene.com
Sat Jun 1 07:48:56 CEST 2013


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



More information about the R-help mailing list