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)