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

Walmes Zeviani walmeszeviani em gmail.com
Domingo Agosto 21 19:17:05 BRT 2011


Alexandre,

Deviance é uma "soma de quadrados" e não um quadrado médio. Então você não
precisa multiplicar pelo grau de liberdade. Veja

y <- rnorm(20)
m0 <- lm(y~1)
m1 <- glm(y~1, family=gaussian)
deviance(m0)
deviance(m1)
summary(m0)$sigma^2*df.residual(m0)
summary(m1)$deviance

Essas coisas de R² são complicadas. Pessoas da velha escola (ou não) sempre
adorar uma medida em porcentagem. Deve trazer um certo conforto, sei lá.
Penso que você deva argumentar outros aspectos da análise do que tentar
justificar que tá tudo ok porque uma espécia de R² foi alto (ou não). De
qualquer forma, uma media seria a porcentagem de deviance do modelo final
contra o modelo nulo, algo como

x <- 1:30
y <- rpois(x, x/4)
m0 <- glm(y~1, family=poisson)
m1 <- glm(y~x, family=poisson)
1-deviance(m1)/deviance(m0)

À disposição.
Walmes.

==========================================================================
Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
Departamento de Estatística - Universidade Federal do Paraná
fone: (+55) 41 3361 3573
VoIP: (3361 3600) 1053 1173
e-mail: walmes em ufpr.br
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110821/c89726f7/attachment.html>


Mais detalhes sobre a lista de discussão R-br