我的模型是时间序列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