[R-sig-ME] [R] [lme4]Coef output with binomial lmer

Doran, Harold HDoran at air.org
Fri Aug 8 20:27:34 CEST 2008


The extractor function for the fixed effects is fixef(), not coef(). Out
of curiosity, why are you using (1|ass%in%pop%in%fam)? This notation is
non-standard and does not define the nesting structure of the data.

I think you want (1|ass/pop/fam)

> -----Original Message-----
> From: r-help-bounces at r-project.org 
> [mailto:r-help-bounces at r-project.org] On Behalf Of Tom Cameron
> Sent: Friday, August 08, 2008 11:47 AM
> To: r-help at r-project.org; r-sig-mixed-models at r-project.org
> Subject: [R] [lme4]Coef output with binomial lmer
> 
> Dear R users  
>  
> I have built the following model
>  
> m1<-lmer(y~harn+foodn+(1|ass%in%pop%in%fam),family = "quasibinomial")
>  
> where y<-cbind(alive,dead)
>  
> where harn and foodn are categorical factors and the random 
> effect is a nested term to represent experimental structure 
> e.g.  Day/Block/Replicate ass= 5 level factor, pop= 2 
> populations per treatment factor in each assay, 7 reps per population
>  
> The model can be family = quasibinomial or binomial
>  
> My complete lack of understanding is in retrieving the 
> coefficients for the fixed effects to back-transform the 
> effects of my factors on proportional survival
>  
> I get the following output:
> > coef(m1)
> $`ass %in% pop %in% fam`
>       (Intercept)      harn1     harn2   foodn2
> FALSE   1.0322375 -0.1939521 0.0310434 0.810084
> TRUE    0.5997679 -0.1939521 0.0310434 0.810084
>  
> Where FALSE and TRUE refer to some attribute of the random effect 
>  
> My hunch is that it refers to the Coefficients with (=TRUE) 
> and without
> (=FALSE) the random effects?
>  
> Any help appreciated
>  
> 
> ..............................................................
> ..........
> ............
> Dr Tom C Cameron
> Genetics, Ecology and Evolution
> IICB, University of Leeds
> Leeds, UK
> Office: +44 (0)113 343 2837
> Lab:    +44 (0)113 343 2854
> Fax:    +44 (0)113 343 2835
> 
> 
> Email: t.c.cameron at leeds.ac.uk
> Webpage: click here
> <http://www.fbs.leeds.ac.uk/staff/profile.php?tag=Cameron_TC> 
> 
>  
> 
> 	[[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.
> 




More information about the R-sig-mixed-models mailing list