[R-sig-ME] Unable to fit glmer with family=binomial(link=identity)

[Doran, Harold]

I may be missing something here, but why are you using an identity link with binomial outcomes? If there is a reason why, you might as well just use lmer since the identity link is for a normal distribution and that is the distributional assumption used for the errors in that function.

 

> -----Original Message-----
> From: r-sig-mixed-models-bounces at r-project.org 
> [mailto:r-sig-mixed-models-bounces at r-project.org] On Behalf 
> Of David Evans
> Sent: Wednesday, August 19, 2009 11:03 AM
> To: r-sig-mixed-models at r-project.org
> Subject: [R-sig-ME] Unable to fit glmer with 
> family=binomial(link=identity)
> 
> Fellow R-users,
> 
> I need to estimate the "risk difference" (in epidemiological
> terminology) for a neighbourhood-level exposure (e.g. 
> presence of parks) with a binary outcome adjusted on 5 or so 
> covariables. Study participants are clusted by neighbourhood 
> .  I was hoping to use a generalized linear mixed model with the code:
> 
> mod  <-  glmer(outcome ~ exposure + (1 | neighbourhood), 
> family=binomial(link=identity), data=rec)
> 
> but the identity link is not available for the binomial 
> family with glmer (or for glmmPQL).  Is there any way to do 
> this with lme4 or is there another package in R which could 
> fix my problem?
> 
> I'd be very grateful for any help.
> 
> David.
> 
> --
> David Evans
> UMR-S 707 Inserm - Université Pierre et Marie Curie - Paris 6 
> Faculté de Médecine Saint-Antoine 27, rue Chaligny
> 75571 Paris cedex 12
> 
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list 
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>