采样的实现中的R

Question

基本上，我想样品9层独立的实现为均匀分布U（0,1）2578倍，这工作正常，即使用

replicate(2578,{runif(9,0,1)}) 

或

F=c() 
for (i in 1:2578){ 
F[i,]=runif(9,0,1) 
} 

现在我想这个被重复让我们说10次，即创建10个新的2578x9样本。我想要创建一个多维数组，或者更好地将其视觉化为长度为9，高度为2578，宽度为（10,1000,100000，...）的长方体。我怎样才能做到这一点？

Answer 1

我觉得你的模拟数据可以被直接构建到一个数组中受益：这会使他们更容易处理：

dims <- c(2578, 9, 100) 
tmp <- runif(prod(dims)) 
A <- array(tmp, dims)