### Olha aí, um GLM da famíla Poisson
modelo <- glm(dano ~ offset(log(tipo)) + construcao+ mes+operacao , family = poisson(link = "log"), data = dt)
summary(modelo)
####Verificação do ajuste do modelo. Faça isto:
1-pchisq(modelo$deviance,modelo$df.residual) ## Probabilidade p
"Deviance"res <- residuals(modelo,type="deviance")
###Verificação de alguns pressuposto do modelo,. Ah, faltou o envelope, hehe!!!
ordem <- 1:length(res)
plot(modelo$fitted.values, res, xlab='fitted', ylab = 'Desvio Residual')
plot(ordem, res, xlab='ordem', ylab = 'Desvio Residual')
hist(res)
qqnorm(res); qqline(res);
sexta-feira, 9 de janeiro de 2009
Assinar:
Postar comentários (Atom)
Nenhum comentário:
Postar um comentário