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，我不太会设置数据，求大神帮忙，不甚感激！