[Q & A]Dirichlet-Multinomial WinBUGS?

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

So I'm trying to code a dirichlet-multinomial model using BUGS. Basically I have 18 regions and 3 categories per region. In example, Region 1: 0.50 belongs to Low, 0.30 belongs to Middle, and 0.20 belongs to High. The list goes on to Region 18 of course with varying proportions.The only code I got is as follows. I shortened it to 3 regions first just for example. It states 'Dirichlet36' after clicking 'gen inits'. Please help me to code this.

1. model  {
   
2. for (i in 1:N) {
   
3. x[1:3] ~ dmulti(p[],n[i])
   
4. p[1:3] ~ ddirch(alpha[])
   
5. }
   
6. for (k in 1:3) {
   
7. alpha[k] <- 1
   
8. }
   
9. }
   
10. 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)`

复制代码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：Multinomial Dirichlet winbugs WINBUG nomial example course trying

This isn't clear - if the x are multinomial they should be integers. Also the code you've given doesn't reproduce the error you've stated. It looks like the x and the p need to be matrices, e.g. x[i,1:3], p[i,1:3], and if x is supposed to be a matrix, it needs to be defined in the data as something like structure(.Data=..., .Dim=c()) –