<div dir="ltr"><div>Prezados(as)</div><div>Para modelos lineares generalizados, o R oferece opções "inverse" </div><div>e "power(-1)" para funções de ligação da distribuição gama. No </div><div>meu entendimento, estas duas opções deveriam produzir </div><div>os mesmos resultados de ajuste. Mas isso não acontece.</div><div>veja o CMR seguinte. Alguém poderia me comentar essa situação?</div><div><br></div><div><br></div><div><div>rm(list=ls(all=TRUE))</div><div>a=structure(list(x = c(0, 0, 0, 1, 1, 1, 2, 2, 2, 3, 3, 3, 4, 4, </div><div>5, 5, 5), y = c(5, 7, 9, 7, 10, 8, 13, 11, 15, 19, 13, 21, 33, </div><div>22, 47, 33, 59)), .Names = c("x", "y"), row.names = c(NA, -17L</div><div>), class = "data.frame")</div><div>attach(a)</div><div>mod1<-glm(y ~ x,family=Gamma(link=power(-1)),data=a)</div><div>summary(mod1)</div><div>mod2<-glm(tvida~1,family=Gamma(link=inverse),data=a)</div><div>summary(mod2)</div></div><div><br></div><div><br></div><div>Grato<br></div><div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr">Gilênio Borges Fernandes<div>Professor Associado IV (Aposentado)</div><div>Professor Adjunto A (Substituto)<br>Universidade Federal da Bahia<br>Instituto de Matemática<br>Departamento de Estatística<br>Av. Adhemar de Barros, s/n – Ondina.<br>40.170-110 - Salvador - BA, Brasil<br>Tel.: (071)3283-6340/6341/6337  Fax:  (071)3283-6336</div><div>Skype: gilenio.fernandes<br>Lattes: <a href="http://lattes.cnpq.br/6764860618464860" target="_blank">http://lattes.cnpq.br/6764860618464860</a><br></div></div></div></div></div></div></div></div></div></div>
</div>