r - 在 R 中以矩阵形式生成引导样本的快速函数

Question

我有一个矩阵 A，我想从 A 的每一列中抽取样本并构建新矩阵。例如:

A = matrix(seq(1,9),3,3)

因此，为了获得第一个引导矩阵，我将从 A 的第一列进行替换(3 次)采样，即 1,2,3，从中进行替换(3 次)采样A 的第二列，即 4、5、6，以及来自 A 的第三列的替换样本(3 次)，即 7、8、9。之后，我通过组合三个自举向量重新构造第一个自举矩阵 B1。我将重复此过程 B=199 次，以便引导矩阵 B1,...,B199 可用。

我的问题是，如何让这个程序运行得更快？我应该使用哪个功能？我知道 apply 本质上涉及 for 循环，因此无法保证速度。 do.call 怎么样？谢谢!

Answer 1

您可以使用replicate 和apply:

A <- matrix(seq(1,9),3,3)

B <- 199

replicate(B, apply(A, 2, sample, replace = TRUE))