library(boot) intadd<-read.csv("internet.csv") bf <- function(intadd, index){ dd <- intadd[index,] return(cor(dd$f1, dd$f2)) } set.seed(123) corr_b <- boot(intadd, bf, R=10000) print(corr_b) plot(corr_b) boot.ci(corr_b,conf=0.95,type="bca") summary(corr_b)