初学winbugs，我的数据咋么运行不了，求助！

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

model
   {
       for( i in 1 : N ) {
           WSG[i]~dnorm(mu,tau)
           tau~dgamma(r,s)
           mu~dnorm(0.5,1000000)
                  }
        r~dnorm(0.0,1.0E-6)
        s~dnorm(0.0,1.0E-6)
        sigma <- 1 / sqrt(tau)
   
    }

Data
   list(WSG=c(0.976,0.946,0.964,0.886,1.3075), N=5)   
inits
   list(r=0.5,s=100000)

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：winbugs WINBUG BUGS Win bug sigma

你的程序问题很多，以致于我不知如何改，问题如下：
1. tau, mu 在for的循环内，就会使tau, mu重复被定义5遍，导致程序无法运行，所以需要将它们的定义移到循环以外。
2. mu~dnorm(0.5, 1000000)这样的定义Mu，使Mu服从MEAN=0.5，VAR=1E-6的正态分布，跟定义MU为常数0.5没什么区别。当然不影响程序运行。
3. Tau是GAMMA分布，GAMMA分布的参数要求为正数,即r>0, s>0。而你定义r, s 是正态（MEAN=0），可以为负，所以程序运行会出错；
4.s=100000这个初始值与s符合MEAN=0，VAR=1000000（即标准差=1000）的分布不匹配，即符合这样分布的s不可能取到这么大的正数，这个初值设定很奇怪。
5. 就5个数据，你的模型却那么复杂。数据量不太够，估计出来的意义不大。当然，我改了上面的4个毛病后，还是能运行的。

TimeT 发表于 2013-12-11 21:10
你的程序问题很多，以致于我不知如何改，问题如下：
1. tau, mu 在for的循环内，就会使tau, mu重复被定义5 ...

求助TimeT大神！！！
我想对BDI指数建立SV模型。用的是论坛里下载的程序，但是Data和INITIAL VALUES就有点摸不着头脑了。
现在Data这里是：list(N=538,y = c(-0.000680522,0.000476216,-0.002995338,……），y用的是BDI的收益率。
INITIAL VALUES是list(itau2 = 7, mu = 0, phi1= 0.975)。
模型正确，成功导入数据，也compile了，gen inits了。
可是到了inference的sample monitor tool，不能设置node，提示：monitor could not be set。
麻烦大神告知，这是怎么回事！

陈金海 发表于 2013-12-18 21:15
求助TimeT大神！！！
我想对BDI指数建立SV模型。用的是论坛里下载的程序，但是Data和INITIAL VALUES就有 ...

我算不上大神。你的问题不太完整，没想出如何回答。

TimeT 发表于 2013-12-11 21:10
你的程序问题很多，以致于我不知如何改，问题如下：
1. tau, mu 在for的循环内，就会使tau, mu重复被定义5 ...

我按照你的要求，重新修改了下，现在能运行，麻烦你看下还有什么错误没？
model
   {
       for( i in 1 : N ) {
           WSG~dnorm(mu,tau)I(0.15,1.3)
           
                  }
        mu~dnorm(mu0,tau0)I(0.15,1.3)
      
        tau~dgamma(r,s)
        r~dnorm(0.0,1.0E+6)I(0.000001,)
        s~dnorm(0.0,1.0E+6)I(0.000001,)
        sigma <- 1 / sqrt(tau)
   
    }


Data
   list(WSG=c(0.976,0.946,0.964,0.886,1.3075), N=5, mu0=0.5,tau0=1.0E+6)   

inits
   list(tau=1,mu=1,r=0.1,s=0.1)

xuyaozhan 发表于 2013-12-19 12:54
我按照你的要求，重新修改了下，现在能运行，麻烦你看下还有什么错误没？
model
   {

第三行的 WSG~dnorm(... 应该是：WSG后面有 [ i ]

TimeT 发表于 2013-12-18 22:00
我算不上大神。你的问题不太完整，没想出如何回答。

您好，能指点下WinBUGS吗？我做结构方程模型，但是初学

我也有monitor 的问题，请问那个是什么意思啊，为什么用里面的参数就可以set呢？用你的那个类似的wsg就不行

陈金海 发表于 2013-12-18 21:15
求助TimeT大神！！！
我想对BDI指数建立SV模型。用的是论坛里下载的程序，但是Data和INITIAL VALUES就有 ...

您好，请问monitor could not be set问题解决了吗？

Try this


model {
        for (i in 1:N) {
                WSG[i] ~ dnorm (mu,tau[i])
                tau[i] ~ dgamma (r,s)
                }
   tau_mean <- r/s
        tau_var<-r/(s*s)
        mu ~ dnorm(0.5,100)
        r <- abs(r1)
        s <- abs(s1)
        r1 ~ dnorm(1, 100)
   s1 ~ dnorm(1, 100)
}

Data
list(N=5, WSG=c(0.976,0.946,0.964,0.886,1.3075))

Inits
list(mu=0.5, r1=1, s1=1)