# fit a model
mod = lm(charges ~ bmi + sex, data = RESI::insurance)
# run resi with the store.boot = TRUE option
resi_obj = resi(mod, nboot = 100, store.boot = TRUE, alpha = 0.01)
# run summary, specifying a different alpha level if desired
summary(resi_obj, alpha = 0.05)
Run the code above in your browser using DataLab