利用R做Bootstrap统计分析

介绍一个做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&