求助大神:下面这个程序哪里出错了?
model volatility;
{
#distribution of Ys
###################
for (i in 1:N) {
ysigmadet[i]<-exp(th[i,1]+th[i,2])*(1-rhoep[i]*rhoep[i]);
Yisigma2[i,1,1] <- exp(th[i,2])/ysigmadet[i];
Yisigma2[i,2,2] <- exp(th[i,1])/ysigmadet[i];
Yisigma2[i,1,2]<-rhoep[i]*exp(0.5*th[i,1]+0.5*th[i,2])/ysigmadet[i];
Yisigma2[i,2,1] <- Yisigma2[i,1,2];
Y[i,1:2]~ dmt(muy[],Yisigma2[i,,],v);
c1[i+1]<-Y[i,1:2];
p[i]<-pow((1+exp(delta0+delta1*c1[i+1])),-1);
s[i]~dbern(p[i]);
t[i]<- s[i]+1
}
muy[1]<-0;
muy[2]<-0;
thmean[1,1] <- mu1[1];
thmean[1,2] <- mu2[1];
th[1,1]~dnorm(thmean[1,1],itaua2[1]);
th[1,2]~dnorm(thmean[1,2],itaub2[1]);
sig1[1]<-exp(0.5*th[1,1]);
sig2[1]<-exp(0.5*th[1,2]);
q[1]~dnorm(psi0[1],itau2[1]);
rhoep[1]<-(exp(q[1])-1)/(exp(q[1])+1);
for (i in 2:N) {
thmean[i,1] <- mu1[t[i]] + phi1[t[i]]*(th[i-1,1]-mu1[t[i]]);
thmean[i,2] <- mu2[t[i]] + phi2[t[i]]*(th[i-1,2]-mu2[t[i]]);
th[i,1]~dnorm(thmean[i,1],itaua2[t[i]]);
th[i,2]~dnorm(thmean[i,2],itaub2[t[i]]);
sig1[i]<-exp(0.5*th[i,1]);
sig2[i]<-exp(0.5*th[i,2]);
qmean[i]<-psi0+psi*(q[i-1]-psi0);
q[i]~dnorm(qmean[i],itau2);
rhoep[i]<-(exp(q[i])-1)/(exp(q[i])+1);
}
#distribution of phi, mu, rhoep
###########################
phi1star ~ dbeta(20,1.5);
phi1[1] <- 2*phi1star -1;
phi1[2] <- 2*phi1star -1;
phi2star ~ dbeta(20,1.5);
phi2[1] <- 2*phi2star -1;
phi2[2] <- 2*phi2star -1;
psistar ~ dbeta(20,1.5);
psi <- 2*psistar -1;
vstar~ dchisqr(4);
v<- vstar*2;
itaua2[1] ~ dgamma(2.5,0.025);
itaua2[2] ~ dgamma(2.5,0.025);
taua[1] <- sqrt(1/itaua2[1]);
taua[2] <- sqrt(1/itaua2[2]);
itaub2[1] ~ dgamma(2.5,0.025);
itaub2[2] ~ dgamma(2.5,0.025);
taub[1] <- sqrt(1/itaub2[1]);
taub[2] <- sqrt(1/itaub2[2]);
itau2 ~ dgamma(2.5,0.025);
tau <- sqrt(1/itau2);
mu1[1] ~ dnorm(0,0.04);
mu1[2] ~ dnorm(0,0.04);
mu2[1] ~ dnorm(0,0.04);
mu2[2] ~ dnorm(0,0.04);
psi0~dnorm(0.7,0.1);
delta0~dnorm(0,0.0001);
delta1~dnorm(0,0.0001)
}
list(phi1star=0.99,phi2star=0.99,mu1=0,mu2=0,itaua2=100,itaub2=100,psistar=0.99,psi0=1.9,itau2=100,vstar=4)
第一步模型检查就错了,程序报错:expected left pointing arrow <- or twiddles ,然后在“ qmean[i]<-psi0+psi*(q[i-1]-psi0);”的“*”中出现一个小方框,我看这一步也没什么问题啊,求大神解答。