我的模型是时间序列ARMA(1,1)模型，用openbugs的check model 和load data 都通过了，但是compile的时候显示array index is greater than array upper bound than Y,请教是我输入的数据有问题还是代码有问题。附代码:model
{
for(t in 1:N+1){
    Y[1]~dnorm(mu[T],tau);}
    mu[1]<-a[1]+M
for(t in 2:N){
     mu[t]<-a[t]+theta1*Y[t-1]+phi1*a[t-1];
                   }

M~dnorm(0,1.0E-2)
tau~dgamma(0.1,0.1);

for(t in 1:N+1){
       a[t]~dnorm(0,tau)
                   }
beta1~dbeta(1,1);
theta1<-2*beta1-1
gamma1~dbeta(1,1);
phi1<-2*gamma1-1

}

输入数据的格式就是list(N=2276,
mu=c(0.0·······））


想问一下哪里需要修改