# NOT RUN {
data(musdata)
coords <- musdata[ ,4:5]
mglm <- glm(musculus ~ pollution + exposure, "poisson", musdata)
aic <- aic.calc(musculus ~ pollution + exposure, "poisson", musdata,
mglm$fitted)
aic$AIC
# }
Run the code above in your browser using DataLab