为何我这个程序到第三步运行不了
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 本人急需解决
解决此问题