
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@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================