e "power(-1)" para funções de ligação da distribuição gama. No
os mesmos resultados de ajuste. Mas isso não acontece.
veja o CMR seguinte. Alguém poderia me comentar essa situação?
rm(list=ls(all=TRUE))
a=structure(list(x = c(0, 0, 0, 1, 1, 1, 2, 2, 2, 3, 3, 3, 4, 4,
5, 5, 5), y = c(5, 7, 9, 7, 10, 8, 13, 11, 15, 19, 13, 21, 33,
22, 47, 33, 59)), .Names = c("x", "y"), row.names = c(NA, -17L
), class = "data.frame")
attach(a)
mod1<-glm(y ~ x,family=Gamma(link=power(-1)),data=a)
summary(mod1)
mod2<-glm(tvida~1,family=Gamma(link=inverse),data=a)
summary(mod2)