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