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.