5

boot()Rのパッケージの関数を使用して、boot2つのベクトル間のピアソン相関係数の有意性を計算する場合、次のように実行する必要があります。

boot(re1, cor, R = 1000)

re1これらの2つの観測ベクトルの2列行列はどこにありますか？corこれらのベクトルがであるため、これを正しく理解できないようです0.8が、上記の関数はとして返さ-0.2れt0ます。

4

1 に答える 1

7

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

これは、パッケージ内のlow birth baby調査での2つの実用的な例です。MASS

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
于 2011-10-20T15:36:02.307 に答える