У меня есть вектор Z
длины K Y с размерами N x S из вектора Z
.
У меня есть другой вектор n.positive
длины S. n.positive[s]
говорит мне, сколько значений в Z
нужно вставить в строку s матрицы Y
(остальные записи равны нулю).Следовательно, sum(n.positive)
равно K.
Я хотел бы сделать это в R.
Например,
y <- matrix(rep(0, N*S), nrow=S)
n.pos <- rbinom(S, N, theta)
Z <- rnorm(sum(n.pos))