问高手关于winbugs中不同先验分布比较的问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

小弟求教，不知道坛子里玩winbugs的人多不，最近在写一个需要比较在不同先验分布下后验的结论的模型，本来分别定义了三个先验，然后分别用bugs估计，现在想能否直接在同一个模型中定义一个循环，把三个先验都放在里面，但得出的结论和之前有很大差别，不知是不是有什么地方写错，请高人指点，bugs 的code如下：

model {
   for (n in 1:3) {  #定义有三个先验的循环         
         for (j in 1:K) {
                 for (i in 1:M) {
                          alpha[i,j]~dnorm(mu[n,j], tau[n])
           }
           mu[n,j] <- t[n,j]+a[n,j]   
          t[n,j]~dnorm(m[n],V[n])
          a[n,j]~dnorm(0,sigma[n])
                            }                  
            }
# 先验分布1
tau[1]~dgamma(100,0.0075)   
sigma[1]~dgamma(10,0.00033)
m[1]~dnorm(0,400)
V[1]~dgamma(100,0.0075)

# 先验分布2
tau[2]~dgamma(10,0.03)   
sigma[2]~dgamma(1,0.0013)
m[2]~dnorm(0,100)
V[2]~dgamma(10,0.03)

# 先验分布3
tau[3]~dgamma(0.001,0.001)   
sigma[3]~dgamma(0.001,0.001)
m[3]~dnorm(0.001,0.001)
V[3]~dgamma(0.001,0.001)
}

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：winbugs WINBUG 先验分布 BUGS bug