Мне нравится рисовать nrep
раз из биномиального распределения с тета-параметром, чтобы создать одну последовательность длины k
для каждой тета и построить их в матричном измерении nrep x k
. Как создать матрицу результатов в R?
Приведенный ниже код рисует (nrep * k
) каждый из разных тета, т.е. последовательность не имеет длины 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)