Obrigado Walmes e Benilton foi bem esclarecedor. 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. Redobrados agradecimentos, Alexandre dos Santos Engenheiro Florestal, MSc. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 3829-5122 ________________________________ De: Benilton Carvalho <beniltoncarvalho@gmail.com> Para: r-br@listas.c3sl.ufpr.br; Alexandre Santos <alexandresantosbr@yahoo.com.br> Enviadas: Domingo, 21 de Agosto de 2011 19:19 Assunto: Re: [R-br] Coeficiente de correlação em glm é correto? voce pode, de fato, usar algo como um R^2 marginal: R2M = 1 - soma((y - y.predito)^2)/soma((y-y.barra)^2) preocupa-me, entretanto, sua definicao do modelo. para familia (quasi-)binomial, a variavel resposta e' algo binario ou uma matrix de 2 colunas com o numero de sucessos/fracassos. par/100 definitivamente nao e' binario, a menos que a variavel (nao-descrita) assuma apenas 0 ou 100. b 2011/8/21 Alexandre Santos <alexandresantosbr@yahoo.com.br>:
Boa noite pessoal,
Venho colocar uma questão para os colegas estatísticos, até onde sei os modelos ajustados com GLM não são representados com o R^2. Porém, normalmente em revistas de entomologia (minha área) contam com editores da velha escola estatística que ficam batendo o pé para se colocar o R2 em todos os modelos ajustados, onde coloco a questão esta correto? Bom, tentei realizar o cálculo do R^2 usando aquelas abordagens que inúmeras vezes apareceram nas listas do grupo do R sobre R^2 para modelos não lineares e fiz o seguinte cálculo:
##Ajustando o modelo com distribuição de erros binomial
m.model4<-glm((par/100)~densidade+I(densidade^2)+I(densidade^3), data=dados2,family="quasibinomial")
SQEm1<-summary(m.model4)$deviance^2*summary(m.model4)$df.residual#Soma de quadrados do erro SQEm1 [1] 2387.366 SQTm1<-var(par)*(length(par)-1)#Soma de quadrados total corrigida R2<-1-SQEm1/SQTm1 R2 [1] 0.8679119
Diante do problema apresentado isso esta correto? Obrigado,
Alexandre dos Santos Engenheiro Florestal, MSc. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 3829-5122
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
-- Successful people ask better questions, and as a result, they get better answers. (Tony Robbins)