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 <[email protected]>: > 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 > [email protected] > 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) _______________________________________________ R-br mailing list [email protected] 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.
