Hatena::ブログ(Diary)

極めて個人的なメモ Twitter

Θ・)ノシ Bob#3のメモ帳です。
1972 | 12 |
2003 | 03 | 04 | 05 | 06 | 11 | 12 |
2004 | 01 | 02 | 03 | 04 | 05 | 06 | 09 | 10 | 11 | 12 |
2005 | 01 | 02 | 03 | 04 | 05 | 06 | 07 | 09 | 10 | 11 |
2006 | 03 | 04 | 05 | 06 | 07 | 09 | 10 | 11 | 12 |
2007 | 01 | 02 | 03 | 05 | 06 | 07 | 10 | 11 | 12 |
2008 | 01 | 02 | 04 | 05 | 08 | 09 | 10 | 11 |
2009 | 01 | 02 | 03 | 04 | 05 | 06 | 07 | 08 | 09 | 10 | 11 | 12 |
2010 | 01 | 02 | 03 | 04 | 05 | 06 | 07 | 08 | 09 | 10 | 11 |
2011 | 05 | 06 | 09 | 10 |
2012 | 02 | 03 | 05 | 06 | 08 |
2013 | 01 | 02 | 03 |
2015 | 02 |

2008-10-02

[]Rでブートストラップ Rでブートストラップを含むブックマーク Rでブートストラップのブックマークコメント

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

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 
トラックバック - http://d.hatena.ne.jp/bob3/20081002