从dirichilet多项式分布示意

Question

我想创建一个1000行和91列的矩阵，填充了由Dirichilet多项式分布生成的整数。

到目前为止，我所取得的成就是可以创建一个矩阵（带91列的1行）

library(extraDistr)
mypar<-c(0.1,0.1,0.2,0.3,0.1,0.2)

mysim <- rdirmnom(91, size=1, mypar)

#changing 1 to item number
mysim[,2] <- ifelse(mysim[,2]==1,2,0)
mysim[,3] <- ifelse(mysim[,3]==1,3,0)
mysim[,4] <- ifelse(mysim[,4]==1,4,0)
mysim[,5] <- ifelse(mysim[,5]==1,5,0)
mysim[,6] <- ifelse(mysim[,6]==1,6,0)

tokeep <- which(mysim>0)
mychoice.sim <- mysim[tokeep]
rep.y <- t(matrix(mychoice.sim))

我将非常感谢有关如何实现预期结果的任何建议。

stef

可能有很多方法可以做到这一点。这是一个：

Answer 1