[疑难杂症]求助ＥＧＡＲＣＨ ＣＯＤＥ！！！

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

請問各位大大誰會EGARCH　Winbugs code！！！
能否借小弟參考一下！！！
感激不盡~

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：疑难杂症 winbugs EGARCH WINBUG GARCH

lz程序编出来没？小弟编了一个，但是貌似没有compile出来。

model{
for(t in 3:T){
         y[t]~dnorm(av[t],p[t])
         LL[t]<--0.5*log(6.28)-0.5*log(h[t])-0.5*pow(y[t]-av[t],2)/h[t]
         }

for(t in 4:T){h[t]<-alpha0+alpha[1]*pow(y[t-1]-av[t-1],2)+beta[1]*h[t-1]
         av[t]<-rho[2]*y[t-2]+rho[3]*y[t-3]
         p[t]<-1/h[t]}

         h3~dgamma(2,1000)
         h[3]<-h3
         p[3]<-1/h3
         av[3]~dnorm(0,0.001)
      
         alpha0~dnorm(0,0.001)I(0,100)
         alpha[1]~dnorm(0,0.001)I(0,1)
         beta[1]~dnorm(0,0.001)I(0,1)
         rho[2]~dnorm(0,0.001)I(0,1)
         rho[3]~dnorm(0,0.001)I(0,1)
}

這個可以
model
{
h0<-exp(w/(1-b))
e0~dnorm(0,1)  
for (t in 1:n)
{
     e[t]~dnorm(0,1)
     }
h[1]<-pow(h0,b)*exp(w+alpha*abs(e0)-alpha*0.79788+theta*e0)
p[1]<-1/h[1]
for (t in 1:n)
    {
     y[t]~dnorm(0,p[t])
                       }
for (t in 2:n)
    {
     h[t]<-pow(h[t-1],b)*exp(w+alpha*abs(e[t-1])-alpha*0.79788+theta*e[t-1])
     p[t]<-1/h[t]
                          }

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

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