utils::data(Pima.te, package = "MASS")
r.logit <- glm(type ~ ., data=Pima.te, family="binomial")
# calculate Brier score with confidence intervals
Cstat(r.logit)
CstatCI(r.logit, R=99) # use higher R in real life!
Run the code above in your browser using DataLab