Internal function to get prediction intervals from glm model
Usage
calc_pi_nb_analytic(
df,
fit,
quantile_levels = c(0.01, 0.025, seq(0.05, 0.95, by = 0.05), 0.975, 0.99)
)
Arguments
- df
dataframe to get predictions on
- fit
model fit via MASS::glm.nb
- quantile_levels
quantiles that you want prediction intervals at
Value
matrix of prediction intervals for each row of df
Details
This is based on the trending
package, as described here: https://cran.r-project.org/web/packages/trending/vignettes/prediction_intervals.html