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@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@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