5

ブートストラップ統計の値（元の値、バイアス、エラー）を別のリストに入れたいのですが、その方法がわかりません。

次に例を示します。

> library(boot)
> set.seed(123)
> mean.fun <- function(data, idx) { mean(data[idx]) }
> data <- boot(data=rnorm(100), statistic=mean.fun, R=999)
> names(data)
 [1] "t0"        "t"         "R"         "data"     
 [5] "seed"      "statistic" "sim"       "call"     
 [9] "stype"     "strata"    "weights"  
> data
ORDINARY NONPARAMETRIC BOOTSTRAP


Call:
boot(data = rnorm(100), statistic = mean.fun, R = 999)


Bootstrap Statistics :
    original   bias    std. error
t1* 0.09040591 0.004751773  0.08823615

ここで、テキストの代わりに実際の値が必要です。どうやらdata$t0「オリジナル」ですが、バイアスとエラーの値を取得する方法がわかりません。

また、関数名を入力するとそのコードが得られるのでboost、Rと入力してソースコードからスニペットをコピーし、ローカルのRインストールで検索しようとしました。しかし、何も見つかりませんでした。なぜ、Rはローカルストレージからそのソースコードを取得するべきではないのですか？

4

1 に答える 1

1

std.errorとbiasは、ブートオブジェクトの一部として保存されません。その場で計算されます（https://stat.ethz.ch/pipermail/r-help/2011-July/284660.htmlを参照）

あなたの場合、試してみてください：

mean(data$t) - data$t0
sd(data$t)
于 2012-12-03T16:31:08.910 に答える