[R-br] (Resolvido) Coeficiente de correlação em glm é correto?

[Alexandre Santos]

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 em gmail.com>
Para: r-br em listas.c3sl.ufpr.br; Alexandre Santos <alexandresantosbr em 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 em 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 em 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)
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110822/e4a0644a/attachment.html>