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