<html><body><div style="color:#000; background-color:#fff; font-family:arial, helvetica, sans-serif;font-size:12pt"><div><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Arial; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; font-size: 16px;"><span class="Apple-style-span" style="font-family: Monospace; font-size: 13px; line-height: 17px; text-align: -webkit-left; white-space: pre-wrap; "><pre tabindex="0" class="GD40030CLR" style="font-family: Monospace; font-size: 10pt !important; outline-style: none; outline-width: initial; outline-color: initial; border-top-style: none; border-right-style: none; border-bottom-style: none; border-left-style: none; border-width: initial; border-color: initial; white-space: pre-wrap !important; margin-top: 0px;
margin-right: 0px; margin-bottom: 0px; margin-left: 0px; line-height: 1.3; "><span class="GD40030COR ace_keyword" style="white-space: pre; color: blue; "><span style="color: rgb(0, 0, 0);">Boa noite pessoal,</span><br style="color: rgb(0, 0, 0);"><br style="color: rgb(0, 0, 0);"><span style="color: rgb(0, 0, 0);"> Venho colocar uma questão para os colegas estatísticos, até onde sei os modelos ajustados com GLM não </span><br style="color: rgb(0, 0, 0);"><span style="color: rgb(0, 0, 0);">são representados com o R^2. Porém, normalmente em revistas de entomologia (minha área) contam com editores</span><br style="color: rgb(0, 0, 0);"><span style="color: rgb(0, 0, 0);">da velha escola estatística que ficam batendo o pé para se colocar o R2 em todos os modelos ajustados, onde coloco a <br>questão esta correto?</span><span style="color: rgb(0, 0, 0);"> Bom, tentei realizar o cálculo do R^2 usando aquelas abordagens que inúmeras vezes
apareceram <br>nas listas</span><span style="color: rgb(0, 0, 0);"> do grupo do R sobre R^2 para modelos não lineares e fiz o seguinte cálculo:</span><br><br>##Ajustando o modelo com distribuição de erros binomial<br>> </span><span class="GD40030CCR ace_keyword" style="color: blue; ">m.model4<-glm((par/100)~densidade+I(densidade^2)+I(densidade^3), data=dados2,family="quasibinomial")</span><span class="GD40030COR ace_keyword" style="white-space: pre; color: blue; "></span>
<span class="GD40030COR ace_keyword" style="white-space: pre; color: blue; ">> </span><span class="GD40030CCR ace_keyword" style="color: blue; ">SQEm1<-summary(m.model4)$deviance^2*summary(m.model4)$df.residual#Soma de quadrados do erro
</span><span class="GD40030COR ace_keyword" style="white-space: pre; color: blue; ">> </span><span class="GD40030CCR ace_keyword" style="color: blue; ">SQEm1
</span>[1] 2387.366
<span class="GD40030COR ace_keyword" style="white-space: pre; color: blue; ">> </span><span class="GD40030CCR ace_keyword" style="color: blue; ">SQTm1<-var(par)*(length(par)-1)#Soma de quadrados total corrigida
</span><span class="GD40030COR ace_keyword" style="white-space: pre; color: blue; ">> </span><span class="GD40030CCR ace_keyword" style="color: blue; ">R2<-1-SQEm1/SQTm1
</span><span class="GD40030COR ace_keyword" style="white-space: pre; color: blue; ">> </span><span class="GD40030CCR ace_keyword" style="color: blue; ">R2
</span>[1] 0.8679119<br><br>Diante do problema apresentado isso esta correto? Obrigado,<br></pre></span></span></div><div> </div><div align="center"><font size="2"><b>Alexandre dos Santos<br></b>Engenheiro Florestal, MSc.<br>Universidade Federal de Lavras <br>Departamento de Entomologia <br>Laboratório de Entomologia Florestal <br>Caixa Postal 3037 <br>37200-000 - Lavras/MG<br>Fone: +55 (35) 3829-5122</font><br></div></div></body></html>