如何创建一个循环来生成R中的随机样本列表？

Question

我正在尝试创建一个循环，创建一系列包含随机样本的对象，如下所示：

sample <- ceiling(runif(9, min=0, max=20))

（这是圆形制服的一个例子，但它可以用普通，泊松或任何你想要的东西代替）。

因此，我构建了一个循环，用于自动生成各种生成器，目的是将它们包含在数据框中。 然后，我设计的循环是这样的：

N=50
dep=as.vector(N)
count=1
for (i in 1:N){
    dep[count] <- ceiling(runif(9, min=0, max=20))  
    count=count+1
}

但它不起作用！ 对于每个dep [i]我只有一个数字，而不是九个列表。

我应该怎么做？ 如果我想在数据框中包含每个dep [i]？

非常感谢，我希望你能理解我想要的东西。

Answer 1

这是因为你做了dep向量（这些是默认1D），但你要保存的二维物体在里面。

您可以dep断为NULL和rbind （行绑定）将其在loop.Also，注意，而不是使用count在循环你可以使用i ：

dep <- NULL
for (i in 1:N){
    dep <- rbind(dep,  ceiling(runif(9, min=0, max=20)))
}
# if you look at dep now it's a 2D matrix.
# We'll convert to data frame
dep <- as.data.frame(dep)

但是 ，有一种更简单的方法可以做到这一点。 您不必逐行生成dep ，您可以通过制作包含9*N圆形均匀分布数的向量来dep生成它：

dep <- ceiling(runif(9*N,min=0,max=20))

现在， dep是一个长度为9 * N的向量。 让我们把它变成一个Nx9矩阵：

dep <- matrix(dep,nrow=N)

完成！

所以你可以在一行中完成上面的所有代码：

dep <- matrix( ceiling(runif(9*N,min=0,max=20)), nrow=N )

如果你想要，你可以在dep上调用data.frame （在它被放入2D矩阵形式之后）来获取数据帧。

Answer 2

正如@ mathematical.coffee解释的那样。 但是，在你的情况下似乎对于runif ，你可以使用sample代替。 实际上sample.int更可靠。 ......比使用runif快3倍左右）：

N <- 1000000
system.time( dep <- matrix(sample.int(20, 9*N, replace=TRUE), N) )  # 0.16 secs
range(dep) # 1 20

system.time( dep <- matrix(ceiling(runif(9*N, min=0, max=20)), N) ) # 0.45 secs
range(dep) # 1 20