Site icon R-bloggers

Bootstrap with logistic regression

[This article was first published on r - Brandon Bertelsen, and kindly contributed to R-bloggers]. (You can report issue about the content on this page here)
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.

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")  
return(coef(fit))  
}

boot_fit <- boot(  
   data = your_data, 
   statistic = logit_test, 
   R = 1e5
) 

To leave a comment for the author, please follow the link and comment on their blog: r - Brandon Bertelsen.

R-bloggers.com offers daily e-mail updates about R news and tutorials about learning R and many other topics. Click here if you're looking to post or find an R/data-science job.
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.