2016-05-08 9 views
2

次のコードはベクトルV1をとり、BV1というブートストラップされたサンプルを1つ作成します。私はそれをi回実行し、すべてのBViベクトルを行列MVに配置したいと思います。可能であれば、forループではなくapplyファミリの機能を使用してください。関数を複数回実行し、結果を含む行列を作成する

V1 <- c(0.18, 0.2, 0.24, 0.35, -0.22, -0.17, 0.28, -0.28, -0.14, 0.03, 0.87, -0.2, 0.06, -0.1, -0.72, 0.18, 0.01, 0.31, -0.36, 0.61, -0.16, -0.07, -0.13, 0.01, -0.09, 0.26, -0.14, 0.08, -0.62, -0.2, 0.3, -0.21, -0.11, 0.05, 0.06, -0.28, -0.27, 0.17, 0.42, -0.05, -0.15, 0.05, -0.07, -0.22, -0.34, 0.16, 0.34, 0.1, -0.12, 0.24, 0.45, 0.37, 0.61, 0.9, -0.25, 0.02) 
BV1 <- sample(V1, length(V1), replace=TRUE) 

この行列を使用して、ブートストラップされた要約統計情報の分布を計算します。ご協力いただきありがとうございます。

答えて

1

replicateを使用してsampleを「n」回繰り返してmatrixとして出力できます。

replicate(4, sample(V1, length(V1), replace=TRUE)) 

我々はreplicate

機能(nは、式expr、= "配列" を簡素化) sapply(整数(n)は、eval.parent(代替(関数(...)を見てみると式expr))、 が簡素化=簡素化)

それはsapplyを使用しています(その機能のapply家族のためにOPのニーズが覆われている)

関連する問題