[R-br] Funcoes de ligacao inverse e power(-1) no modelo gama/glm

[Gilenio Borges Fernandes]

Prezados(as)
Para modelos lineares generalizados, o R oferece opções "inverse"
e "power(-1)" para funções de ligação da distribuição gama. No
meu entendimento, estas duas opções deveriam produzir
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)


Grato
Gilênio Borges Fernandes
Professor Associado IV (Aposentado)
Professor Adjunto A (Substituto)
Universidade Federal da Bahia
Instituto de Matemática
Departamento de Estatística
Av. Adhemar de Barros, s/n – Ondina.
40.170-110 - Salvador - BA, Brasil
Tel.: (071)3283-6340/6341/6337  Fax:  (071)3283-6336
Skype: gilenio.fernandes
Lattes: http://lattes.cnpq.br/6764860618464860
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160324/b02a5d3f/attachment.html>