[R-sig-ME] link function with glmer and inverse.gaussian family
Helios de Rosario
helios.derosario at ibv.upv.es
Fri Sep 12 12:16:58 CEST 2014
Hello, I did this:
> mod1 <- glmer(y ~ x + (1|block), data=dat,
family=inverse.gaussian(log))
> getFamily(mod1)
Family: inverse.gaussian
Link function: 1/mu^2
I had expected that the link function would be "log", as I entered it
explicitly. If I fit a glm (without random effects), I do get:
> mod2<- glm(y ~ x, data=dat, family=inverse.gaussian(log))
> getFamily(mod2)
Family: inverse.gaussian
Link function: log
Is there a reason for that different behaviour of the functions?
Thanks
Helios De Rosario
SABEMOS QUÉ HICIMOS EN 2013
365 días de Innovación
ANUARIO DE BIOMECÁNICA
http://www.ibv.org/anuario2013
______________________________
INSTITUTO DE BIOMECÁNICA DE VALENCIA
Universidad Politécnica de Valencia • Edificio 9C
Camino de Vera s/n • 46022 VALENCIA (ESPAÑA)
Tel. +34 96 387 91 60 • Fax +34 96 387 91 69
www.ibv.org
Antes de imprimir este e-mail piense bien si es necesario hacerlo.
En cumplimiento de la Ley Orgánica 15/1999 reguladora de la Protección
de Datos de Carácter Personal, le informamos de que el presente mensaje
contiene información confidencial, siendo para uso exclusivo del
destinatario arriba indicado. En caso de no ser usted el destinatario
del mismo le informamos que su recepción no le autoriza a su divulgación
o reproducción por cualquier medio, debiendo destruirlo de inmediato,
rogándole lo notifique al remitente.
More information about the R-sig-mixed-models
mailing list