Modelos Lineares Generalizados

### 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);