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)
yp1~dnorm(mu1,tau)
mu1<-alpha+rho*y[T]
yp2~dnorm(mu2,tau)
mu2<-alpha+rho*yp1
}
