set.seed(123)
s=.8
n=10000
x <- sort(runif(n))
p <- function(x,s){p = 1/(1+((1/x*(1-x))^(s+1)));return(p)}
dat <- data.frame(pr=x, y=rbinom(n,1,p(x,s)))
cb <- calibration_bands(x=dat$pr, y=dat$y,alpha=0.05, method="round", digits =3)
summary(cb)
print(summary(cb), n=5)
Run the code above in your browser using DataLab