为何我这个程序到第三步运行不了

model

{

for ( t in 1:n) { yisigma2[t] < - 1/exp( theta[t] );

y[t] ~ dnorm( 0, yisigma2[t] );

}

mu ~ dnorm(0, 0.1);

phistar ~ dbeta(20, 1.5);

itau2 ~ dgamma(2.5, 0.025);

 

phi < - 2*phistar-1;

tau < - sqrt( 1/itau2 );

theta0 ~ dnorm( mu, itau2 );

thmean[1] < - mu + phi* theta0 ;

theta[1] ~ dnorm( thmean[1], itau2 );

for ( t in 2:n ) { thmean[t] < - mu + phi* theta[t-1] ;

theta[t] ~ dnorm( thmean[t], itau2 );

}

}

list(y=此处有1503个数据，n＝1503）

list(mu=0,phistar=0,itau=2)

对高手来说这是简单问题   解决此问题奖人币100元    电话13667115104  本人急需解决  

解决此问题