### PI for quasi-Poisson data
pred_int <- quasi_pois_pi(histdat=ames_HCD,
newoffset=3,
nboot=100,
traceplot = FALSE)
### Plot the PI
plot(pred_int)
### Since plot.predint is based on ggplot, the grafic can be altered using
# the methodology provided via ggplot2
plot(pred_int)+
theme_classic()
Run the code above in your browser using DataLab