model{
for(t in 1:T){y[t]~dnorm(m[t],tau.e)}
m[1]<-u[1]+mu
for(t in 2:T){m[t]<-rho*y[t-1]+u[t]+thita*u[t-1]}
for(t in 1:T){u[t]~dnorm(0,tau.u)}
mu~dnorm(0,0.001)
rho~dnorm(0,0.001)
thita~dnorm(0,0.001)
tau.e~dgamma(.001,0.001)
tau.u~dgamma(.001,0.001)
}
