2012-02-16 12 views
3

私は新しいRユーザーです。ブートパッケージの使用に問題があります。ブートストラップの信頼区間R

x <- rnorm(100, 1, .5) 

任意のヒント:私がやりたいすべてのような、数のベクトルの平均の周りに信頼区間を生成するためにブートストラップを使用しているのですか?

答えて

9

次のことは十分ではありませんか?

library(boot) 
x <- rnorm(100, 1, .5) 
b <- boot(x, function(u,i) mean(u[i]), R = 999) 
boot.ci(b, type = c("norm", "basic", "perc")) 
+0

大変感謝しています。ちょうど構文を理解できませんでした。 – mike

+0

「u」と「i」はここで何を表していますか? – cryptic0

+2

@ cryptic0: '?boot'で説明されています。 'u'は元のデータ(この例では' x')であり、 'i'はブーストラップサンプルに対応するインデックスのベクトルです(このサンプルは異なるサンプルに対して' R'回と呼ばれます)。データがベクトルの場合、ブートストラップサンプルは 'u [i]'で、data.frameの場合は 'u [i、]'です。 (もちろん、 'u'と' i'は関数の正式な引数なので、この関数の外では値はありません。) –

関連する問題