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