在R中使用for循环的Replicate（）

Question

我想在R中创建大量样本并将其存储在变量中。 我做了一些研究，可能最好的方法是使用copy（）

record <- replicate(5000, sample(c(0,1), size = 60, replace = T,prob=c(0.9,0.1)))

我的问题是我该如何使用for循环？ 我可以使用for循环创建5000个样本，但是如何存储它们呢？

x <- 'something here' #I want to store them in x

for (i in 1:5000)
    {record <- sample(c(0,1), size = 60, replace = T,prob=c(0.9,0.1))
    'x += record'}

编辑：X + =行记录令人困惑。 这是我最好的解释，这是在python中创建的一个列表，该列表中还有5000个其他列表，每个列表包含一个不同的示例

Answer 1

我看不到使用for循环的意义。 也不需要replicate 。

您可以直接使用绘制5000 * 60独立样本

smpl <- sample(c(0, 1), size = 60 * 5000, replace = TRUE, prob = c(0.9, 0.1))

如果要将smpl存储在矩阵中，则可以将向量重铸为matrix ，例如

mat <- matrix(smpl, ncol = 5000)

这将为您提供60 x 5000矩阵，其中每列包含从p(0) = 0.9和p(1) = 0.1的分布中抽取的5000随机样本。

这将比使用for循环或replicate更快。

Answer 2

关于什么

x <- list() # or x <- c() or x <- data.frame()

for (i in 1:5000){
    record <- sample(c(0,1), size = 60, replace = T,prob=c(0.9,0.1))
    x[[i]] <- record # or x <- c(x, record) or x <- rbind(x, record)
}

但是，是的，如上所述，不使用循环可能会更好。