在R中的n次复制中模拟随机二项式序列

Question

我喜欢使用theta参数从二项分布中绘制nrep次，为每个theta创建一个k长度序列，并在矩阵维nrep x k中构建它们。如何在R？

中创建结果矩阵

下面的代码从不同的theta绘制（nrep * k），即序列不是来自相同theta的k长度。 [我的目的是绘制nrep乘以二项式概率θ长度k。]

### simulate some binary sequence data in matrix ted (1000 x 20)
nrep <- 1000
s <- 7; k <- 20
theta <- rbeta(nrep, shape1=s+1, shape2=k-s+1)

ted <- 0
ted <- matrix(rbinom(k * nrep, 1, theta), ncol = k, nrow = nrep)
hist(ted)

Answer 1

rbinom在参数prob上进行了矢量化，因此您可以使用rep(theta, k)来实现

ted <- matrix(rbinom(k * nrep, 1, rep(theta, k)), ncol = k, nrow = nrep)