### PI for quasi-Poisson data
pred_int <- quasi_pois_pi(histdat=ames_HCD,
newoffset=3,
nboot=100,
traceplot = FALSE)
# Return the prediction intervals as a data.frame
as.data.frame(pred_int)
# Please note that nboot was set to 100 in order to decrease computing time
# of the example. For a valid analysis set nboot=10000.
Run the code above in your browser using DataLab