Skip to contents

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