ブートパッケージ

Question

私は2つのベクトル間のピアソン相関係数の重要性を計算するRのbootパッケージからboot()機能を使用したい場合は、私はこのようにそれを行う必要があります。

boot(re1, cor, R = 1000) 

場所re1は、これらの2つの観測ベクトルの2列の行列ですか？これらのベクトルのうちcorが0.8なので、私はこの権利を得られないようですが、上記の関数は-0.2をt0として返します。

Answer 1

@caracalはすでに彼のコメントを通じてあなたの質問に答えましたが、Rのブートストラップについての一般的な考え方を強調してください。 bootを使用する場合、行でサンプリングできるデータ構造（通常は行列）が必要です。統計の計算は、通常、このデータ行列を受け取り、再サンプリング後に計算された関心の統計を返す関数で行われます。次に、この機能をRに適用して、構造化された形式で結果を複製して収集することを担当するboot()と呼びます。これらの結果は、boot.ci()を順番に使用して評価することができます。

MASSパッケージでlow birth babyの調査を行った2つの作業例があります。

require(MASS) 
data(birthwt) 
# compute CIs for correlation between mother's weight and birth weight 
cor.boot <- function(data, k) cor(data[k,])[1,2] 
cor.res <- boot(data=with(birthwt, cbind(lwt, bwt)), 
       statistic=cor.boot, R=500) 
cor.res 
boot.ci(cor.res, type="bca") 
# compute CI for a particular regression coefficient, e.g. bwt ~ smoke + ht 
fm <- bwt ~ smoke + ht 
reg.boot <- function(formula, data, k) coef(lm(formula, data[k,])) 
reg.res <- boot(data=birthwt, statistic=reg.boot, 
       R=500, formula=fm) 
boot.ci(reg.res, type="bca", index=2) # smoke