<html><body><div style="color:#000; background-color:#fff; font-family:arial, helvetica, sans-serif;font-size:10pt"><div><span>Obrigado Walmes e Benilton foi bem esclarecedor.</span></div><div><span><br></span></div><div><span> E Benilton sobre a variável par que utilizei para ajustar o modelo quasebinomial do exemplo eram valores entre 0 e 100 que dividi por 100.</span></div><div><span>Redobrados agradecimentos,<br></span></div><div> </div><div align="center"><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<br></div><div><br></div><div style="font-family: arial, helvetica, sans-serif; font-size: 10pt;"><div style="font-family: times new roman, new york, times, serif; font-size: 12pt;"><font face="Arial" size="2"><hr
size="1"><b><span style="font-weight:bold;">De:</span></b> Benilton Carvalho <beniltoncarvalho@gmail.com><br><b><span style="font-weight: bold;">Para:</span></b> r-br@listas.c3sl.ufpr.br; Alexandre Santos <alexandresantosbr@yahoo.com.br><br><b><span style="font-weight: bold;">Enviadas:</span></b> Domingo, 21 de Agosto de 2011 19:19<br><b><span style="font-weight: bold;">Assunto:</span></b> Re: [R-br] Coeficiente de correlação em glm é correto?<br></font><br>voce pode, de fato, usar algo como um R^2 marginal:<br><br>R2M = 1 - soma((y - y.predito)^2)/soma((y-y.barra)^2)<br><br>preocupa-me, entretanto, sua definicao do modelo.<br><br>para familia (quasi-)binomial, a variavel resposta e' algo binario ou<br>uma matrix de 2 colunas com o numero de sucessos/fracassos.<br><br>par/100 definitivamente nao e' binario, a menos que a variavel<br>(nao-descrita) assuma apenas 0 ou 100.<br><br>b<br><br>2011/8/21 Alexandre Santos <<a
ymailto="mailto:alexandresantosbr@yahoo.com.br" href="mailto:alexandresantosbr@yahoo.com.br">alexandresantosbr@yahoo.com.br</a>>:<br>> Boa noite pessoal,<br>><br>> Venho colocar uma questão para os colegas estatísticos, até onde sei os<br>> modelos ajustados com GLM não<br>> são representados com o R^2. Porém, normalmente em revistas de entomologia<br>> (minha área) contam com editores<br>> da velha escola estatística que ficam batendo o pé para se colocar o R2 em<br>> todos os modelos ajustados, onde coloco a<br>> questão esta correto? Bom, tentei realizar o cálculo do R^2 usando aquelas<br>> abordagens que inúmeras vezes<br>> apareceram<br>> nas listas do grupo do R sobre R^2 para modelos não lineares e fiz o<br>> seguinte cálculo:<br>><br>> ##Ajustando o modelo com distribuição de erros binomial<br>>>
m.model4<-glm((par/100)~densidade+I(densidade^2)+I(densidade^3),<br>>> data=dados2,family="quasibinomial")<br>><br>>> SQEm1<-summary(m.model4)$deviance^2*summary(m.model4)$df.residual#Soma de<br>>> quadrados do erro<br>>> SQEm1<br>> [1] 2387.366<br>>> SQTm1<-var(par)*(length(par)-1)#Soma de quadrados total corrigida<br>>> R2<-1-SQEm1/SQTm1<br>>> R2<br>> [1] 0.8679119<br>><br>> Diante do problema apresentado isso esta correto? Obrigado,<br>><br>><br>> Alexandre dos Santos<br>> 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<br>><br>> _______________________________________________<br>> R-br mailing list<br>> <a ymailto="mailto:R-br@listas.c3sl.ufpr.br"
href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>> <a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>> Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código<br>> mínimo reproduzível.<br>><br><br><br><br>-- <br>Successful people ask better questions, and as a result, they get<br>better answers. (Tony Robbins)<br><br><br></div></div></div></body></html>