<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=utf-8"><meta name=Generator content="Microsoft Word 15 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:"Calibri Light";
panose-1:2 15 3 2 2 2 4 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:"Adobe Fan Heiti Std B";
panose-1:2 11 7 0 0 0 0 0 0 0;}
@font-face
{font-family:"\@Adobe Fan Heiti Std B";
panose-1:0 0 0 0 0 0 0 0 0 0;}
@font-face
{font-family:"CMU Sans Serif";
panose-1:2 0 6 3 0 0 0 0 0 0;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman",serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.EstiloDeEmail17
{mso-style-type:personal-reply;
font-family:"Arial",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 3.0cm 70.85pt 3.0cm;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=PT-BR link=blue vlink=purple><div class=WordSection1><p class=MsoNormal><span style='font-family:"Calibri",sans-serif'>Boa noite prezado </span><span style='font-family:"Calibri",sans-serif'>Gilênio.<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Calibri",sans-serif'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-family:"Calibri",sans-serif'>Acrescente as seguintes linhas ao seu programa:<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-family:"Calibri",sans-serif'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-family:"Calibri",sans-serif'>mod3<-glm(y~x,family=Gamma(link=log),data=a)<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Calibri",sans-serif'>summary(mod3)<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Calibri",sans-serif'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-family:"Calibri",sans-serif'>Com isso vc pode verificar que as estimativas do mod1 e do mod3 são idênticas. <o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Calibri",sans-serif'>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.<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Calibri",sans-serif'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-family:"Calibri",sans-serif'>Espero ter ajudado.<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Calibri",sans-serif'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-family:"Calibri",sans-serif'>Att,<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Calibri",sans-serif'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-family:"Calibri",sans-serif'>Davi Butturi-Gomes<o:p></o:p></span></p><p class=MsoNormal><b><span style='font-family:"Calibri Light",sans-serif'><o:p> </o:p></span></b></p><p class=MsoNormal><b><span style='font-family:"Calibri Light",sans-serif'><o:p> </o:p></span></b></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Arial",sans-serif;mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><b><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'>De:</span></b><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'> R-br [mailto:r-br-bounces@listas.c3sl.ufpr.br] <b>Em nome de </b>Gilenio Borges Fernandes<br><b>Enviada em:</b> quinta-feira, 24 de março de 2016 17:04<br><b>Para:</b> r-br@listas.c3sl.ufpr.br<br><b>Assunto:</b> [R-br] Funcoes de ligacao inverse e power(-1) no modelo gama/glm<o:p></o:p></span></p><p class=MsoNormal><o:p> </o:p></p><div><div><p class=MsoNormal>Prezados(as)<o:p></o:p></p></div><div><p class=MsoNormal>Para modelos lineares generalizados, o R oferece opções "inverse" <o:p></o:p></p></div><div><p class=MsoNormal>e "power(-1)" para funções de ligação da distribuição gama. No <o:p></o:p></p></div><div><p class=MsoNormal>meu entendimento, estas duas opções deveriam produzir <o:p></o:p></p></div><div><p class=MsoNormal>os mesmos resultados de ajuste. Mas isso não acontece.<o:p></o:p></p></div><div><p class=MsoNormal>veja o CMR seguinte. Alguém poderia me comentar essa situação?<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><div><p class=MsoNormal>rm(list=ls(all=TRUE))<o:p></o:p></p></div><div><p class=MsoNormal>a=structure(list(x = c(0, 0, 0, 1, 1, 1, 2, 2, 2, 3, 3, 3, 4, 4, <o:p></o:p></p></div><div><p class=MsoNormal>5, 5, 5), y = c(5, 7, 9, 7, 10, 8, 13, 11, 15, 19, 13, 21, 33, <o:p></o:p></p></div><div><p class=MsoNormal>22, 47, 33, 59)), .Names = c("x", "y"), row.names = c(NA, -17L<o:p></o:p></p></div><div><p class=MsoNormal>), class = "data.frame")<o:p></o:p></p></div><div><p class=MsoNormal>attach(a)<o:p></o:p></p></div><div><p class=MsoNormal>mod1<-glm(y ~ x,family=Gamma(link=power(-1)),data=a)<o:p></o:p></p></div><div><p class=MsoNormal>summary(mod1)<o:p></o:p></p></div><div><p class=MsoNormal>mod2<-glm(tvida~1,family=Gamma(link=inverse),data=a)<o:p></o:p></p></div><div><p class=MsoNormal>summary(mod2)<o:p></o:p></p></div></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>Grato<o:p></o:p></p></div><div><div><div><div><div><div><div><div><div><p class=MsoNormal>Gilênio Borges Fernandes<o:p></o:p></p><div><p class=MsoNormal>Professor Associado IV (Aposentado)<o:p></o:p></p></div><div><p class=MsoNormal>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<o:p></o:p></p></div><div><p class=MsoNormal>Skype: gilenio.fernandes<br>Lattes: <a href="http://lattes.cnpq.br/6764860618464860" target="_blank">http://lattes.cnpq.br/6764860618464860</a><o:p></o:p></p></div></div></div></div></div></div></div></div></div></div></div></div></body></html>