bayesian - 狄利克雷多项式 WinBUGS 代码

我正在尝试使用 BUGS 编写狄利克雷多项式模型。
基本上我有 18 个区域和每个区域 3 个类别。例如，
区域 1:0.50 属于低，0.30 属于中，0.20 属于高。列表当然以不同的比例继续到区域 18。我得到的唯一代码是这个

`model  {
 for (i in 1:N) {
 x[1:3] ~ dmulti(p[],n[i])
 p[1:3] ~ ddirch(alpha[])
 }
 for (k in 1:3) {
 alpha[k] <- 1
 }
 }
 DATA list(n=c(38483, 2259, 1900),x=c(29256.42719, 1857.431404, 1548.007808, 29256.42719, 1857.431404, 1548.007808, 29256.42719, 1857.431404, 1548.007808), N=3)`

例如，我首先将其缩短为 3 个区域。单击“gen inits”后，它会显示“Dirichlet36”。请帮我编码这个。

最佳答案

这可能会有所帮助( source ):

关于bayesian - 狄利克雷多项式 WinBUGS 代码，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/24349692/