3

私はRを初めて使用し、頻度データの大きなファイルの標準誤差のブートストラップ推定を試みています。ブートストラップは 1 つのデータ ポイントで正常に動作していますが、出力を保存する方法がわかりません。理想的には、標準エラーのみを新しいファイルに書き込みたいと考えています。

これが私がこれまでに試したことです：

x = c(1,2,3,4,5,6,7,8,9,10)
samplemean = function(x, d){return(mean(x[d]))}
b = boot(x, samplemean, R=1000)
b
ORDINARY NONPARAMETRIC BOOTSTRAP

Call:
boot(data = x, statistic = samplemean, R = 1000)

Bootstrap Statistics :
    original  bias    std. error
t1*      5.5 -0.0356   0.9145759
4

2 に答える 2

2

tオブジェクトの(複製) スロットを使用して、標準誤差を計算できbootます。

require(boot)

x <- c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10)

samplemean <- function(x, d) {return(mean(x[d]))}

set.seed(123)
b <- boot(x,samplemean,R=1000)
b
## Bootstrap Statistics :
##     original  bias    std. error
## t1*      5.5 -0.0232     0.90887

## str(b) first...
apply(b$t, 2, sd)
##[1] 0.90887
于 2012-08-11T00:10:59.547 に答える
1

boot:::print.boot 関数は、次の式を使用して、通常のブートストラップについて報告する「標準エラー」を計算します。

sqrt(apply(t, 2L, function(t.st) var(t.st[!is.na(t.st)])))

ファイルに書き込みたい場合は、それに渡しb$tて使用できますcat：

cat( sqrt(apply(b$t, 2L, function(t.st) var(t.st[!is.na(t.st)]))), file="seboot.txt")

(NA 値を取り出すこれらの行に沿って、関数内で以前にいくつかの処理が行われていました):

t <- matrix(boot.out$t[, index], nrow = nrow(boot.out$t))
allNA <- apply(t, 2L, function(t) all(is.na(t)))
t <- matrix(t[, !allNA], nrow = nrow(t))
于 2012-08-11T00:37:58.490 に答える