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

[Davi Butturi-Gomes]

Boa noite prezado Gilênio.

 

Acrescente as seguintes linhas ao seu programa:

 

mod3<-glm(y~x,family=Gamma(link=log),data=a)

summary(mod3)

 

Com isso vc pode verificar que as estimativas do mod1 e do mod3 são idênticas. 

Pelo que entendi, o argumento lambda da função (de ligação) power precisa ser não-negativo, pois o comportamento é o mesmo para todos os valores menores ou iguais a 0, isto é, age com uma função de ligação logarítmica.

 

Espero ter ajudado.

 

Att,

 

Davi Butturi-Gomes

 

 

 

De: R-br [mailto:r-br-bounces em listas.c3sl.ufpr.br] Em nome de Gilenio Borges Fernandes
Enviada em: quinta-feira, 24 de março de 2016 17:04
Para: r-br em listas.c3sl.ufpr.br
Assunto: [R-br] Funcoes de ligacao inverse e power(-1) no modelo gama/glm

 

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/92b31fdc/attachment.html>