A quick example of bootstraping a logistic regression. Nothing special here, example could be extended to any other type of model that has a coef() method. library(boot) logit_test <- function(d,indices) { d <- d[indices,] fit <- glm(your ~ formula, data = d, family = "binomial"