Ｒでブートストラップ

単純なものですが、半分お遊びで。

boot.mean <- function(x,times=1000){
	set.seed(101)
	b <- numeric(times)
	for(i in 1:times){
		c <- sample(x,replace=T)
		b[i] <- mean(c)
		}
	cat("回数",times,"\n")
	cat("bootstrap平均",mean(b),"\n")
	cat("bootstrap標準偏差",sd(b),"\n")
	cat("95％信頼区間上限",quantile(b,p=0.975),"\n")
	cat("95％信頼区間下限",quantile(b,p=0.025),"\n")
	}

> boot.mean(KG, 2000)
回数 2000 
bootstrap平均 59.33023 
bootstrap標準偏差 11.14857 
95％信頼区間上限 80.46364 
95％信頼区間下限 36.90909