关于用winbugs做GARCH-t模型的代码求解？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

model
{   alpha0 ~ dnorm( 0.0,0.001)I(0,20)
    alpha1 ~ dnorm( 0.0,0.001)I(0,1)
    beta ~ dnorm( 0.0,0.001)I(0,1)
y0 ~ dnorm(0.0,0.001)
k~dunif(1,100)
for( t in 2 : n ) {
         h[t] <- alpha0 + alpha1 * pow((y[t - 1]-y0),2) + beta * h[t - 1]
         tau[t] <- 1 / h[t]
    y[t] ~ dt( y0,tau[t],k)
      }
   y[1] ~ dt( 0,tau[1],k)
   h[1] <- alpha0
  tau[1] <- 1 / h[1]
}
list(alpha0=0.01,alpha1=0.01,beta=0.01,y0=0.01)
list(n=1479,y=c(0.87212,-1.12427,3.0639……）
后面数据省略，一直提醒this chain contains uninitialized variables，我不太会设置数据，求大神帮忙，不甚感激！

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：winbugs WINBUG GARCH BUGS ARCH winbugs

我看了下，发现你的k没有initialize. 你需要在你的第1个list（我猜那是你的INITIAL VALUE）最后加上k=xxx，例如，那个list变成： list(alpha0=0.01,alpha1=0.01,beta=0.01,y0=0.01, k=50)，你试试吧。

我想问一下，dnorm( 0.0,0.001)I(0,20) 中间的 I 是什么意思？