請求各位大神幫我看一下我編寫的Winbugs哪裡出錯了

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

EGARH(1,1)
model
{
for (t in 1:T)
      {
        y[t]~dnorm(0,p[t])
                                     }
for (t in 2:T)
      {
        logsigma2[t]<-omega+beta*logsigma2[t-1]+theta*(y[t-1]/sigma[t-1])+alpha*(sqrt((y[t-1]/sigma[t-1])*(y[t-1]/sigma[t-1]))-sqrt(2/3.14159265359))
        sigma[t]<-exp(logsigma2[t]/2)
        p[t]<-1/sigma[t]
        }

logsigma2[1]<-(omega/(1-beta))
sigma[1]<-exp(logsigma2[1]/2)
p[1]<-1/sigma[1]

# piror
omega~dnorm(0,1)
beta~dbeta(8,1)I(0,1)
theta~dnorm(0,1)
alpha~dnorm(0,1)
}

# initial
list(omega=-0.2, beta=0.98, theta=-0.05, alpha=0.15)

# data
list(T=10, y=c(-0.004935767,
-0.007029977,
-0.004895137,
-0.022116165,
0.001567211,
-0.014067703,
-0.012938876,
0.001584612,
0.010732848,
0.007028073
))               

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：winbugs WINBUG BUGS bug Win

用繁体字，是台湾人吗？给个出错的信息也好帮你找

nuomin 发表于 2014-4-2 20:57
用繁体字，是台湾人吗？给个出错的信息也好帮你找

他出錯的信息是
"sorry something went wrong in procedure Node.Likelihood in module GraphNormal"

garch模型要求方差大于零，考虑这个问题了吗？

nuomin 发表于 2014-4-5 21:08
garch模型要求方差大于零，考虑这个问题了吗？

我修改成exponential的方式，資料數少的時候還能跑，但是資料多的時候就跑不動了
一直找不到問題在哪，資料一多load inits就會當掉了
model
{
for (t in 1:T)
    {
     y[t]~dnorm(0,p[t])
                       }
for (t in 2:T)
    {
     h[t]<-pow(h[t-1],b)*exp(w+alpha*(abs(e[t-1])-0.79788456)+theta*e[t-1])
     p[t]<-1/h[t]
     e[t]<-y[t]/sqrt(h[t])
                          }
h1<-exp(0.2/(1-0.98))
h[1]<-h1
p[1]<-1/h[1]
e1<-y[1]/sqrt(h[1])
e[1]<-e1

# piror
w~dnorm(0,1)
b~dbeta(8,1)I(0,1)
theta~dnorm(0,1)
alpha~dnorm(0,1)
}

# initial
list(w=0.2, b=0.98, theta=0.05, alpha=0.15)

# data
list(T=10, y=c(-0.004935767,
-0.007029977,
-0.004895137,
-0.022116165,
0.001567211,
-0.014067703,
-0.012938876,
0.001584612,
0.010732848,
0.007028073))

用winbugs时要获得随机数，这个随机数有些时候不能满足函数的定义域，把这些处理了，问题也许就解决了。不过OpenBUGS和WinBUGS的可调试性比较差。建议试试MCMC包和Rstan。做得比两个BUGS好一些，也都是开源的

nuomin 发表于 2014-4-6 09:01
用winbugs时要获得随机数，这个随机数有些时候不能满足函数的定义域，把这些处理了，问题也许就解决了。不过 ...

感謝你喔~