#计算第四个观测在自助法样本里的概率
store<-rep(NA,10000)
for (i in 1:10000){store[i]<-sum(sample(1:100,rep=TRUE)==4)>0
}
#讨论第四个观测
mean(store)
#计算第n个观测在自助法样本里的概率
boot <- function(n, N){return(1-(1-1/N)^N)
}#计算第n个观测在自助法样本里的概率
boot(1, 100)
#计算第n个观测不在自助法样本里的概率
oob <- function(n, N){return(1-(1-1/N)^N)
}
#计算第n个观测不在自助法样本里的概率
oob(1, 100)
