如何使用R-JAGS实现具有折叠Gibbs采样的LDA？

Question

我正在尝试使用JAGS实现LDA模型：

model {

    for (i in 1:M) {
      theta[i,1:K] ~ ddirch(alpha)
    }


    for (i in 1:K) {
      phi[i,1:V] ~ ddirch(beta)
    }


    for (i in 1:W) {
      z[i] ~ dcat(theta[D[i],])
      w[i] ~ dcat(phi[z[i],])
    }
   }

该计划已经崩溃，因为对于LDA的Gibbs采样太重了。现在我想用collapsed Gibbs sampling实现它。你知道如何用JAGS以这种方式进行采样吗？提前谢谢！