介绍一个做Bootstrap好用的R 函数,inference。
先看看不用inference时的代码,假定取100个bootstrap样本
boot_means = rep(NA, 100)
for(i in 1:100){
boot_sample = sample(gained_clean, n, replace = TRUE)
boot_means[i] = mean(boot_sample)
}
再来看看inference,载入:
source(<a target=_blank href="http://bit.ly/dasi_inference">http://bit.ly/dasi_inference</a>)
国内如果无法下载,点此链接保存到本地再source,http://d396qusza40orc.cloudfront.net/statistics/lab_resources/inference.R
1、检验一个变量
取10,000个 bootstrap samples,方法是百分数法perc,置信区间0.90
inference(nc$gained, type = "ci", method = "simulation&