model{
for(i in 1:T){
y[i]~dnorm(0,p[i])
p[i]<-exp(-h[i])}
h[1]~dnorm(muh,qh)
muh<-mu1/(1-rho)
qh<-sig*(1-rho*rho)
for(j in 2:T){
h[j]~dnorm(mu2[j],sig)
mu2[j]<-mu1+rho*h[j-1]}
mu1~dnorm(0,0.0001)
sig~dgamma(0.0001,0.0001)
rho~dunif(-0.999,0.999)
}
