下面这段程序怎么运行不起来，求大神指点迷津



model
{                                                               
for (i in 1:n) {
Ymean[i]<-rho/tau*exp(0.5*theta[i])*(theta[i+1]-mu-phi*(theta[i]-mu));
Yisigma2[i] <- 1/(exp(theta[i])*(1-rho*rho));
y[i]~ dnorm(Ymean[i],Yisigma2[i]);   
    }     
isigma2<-(1-phi*phi)*itau2;                                                        
theta0 ~ dnorm(mu,isigma2);  
thetamean[1] <- mu + phi*(theta0-mu);   
theta[1] ~ dnorm(thetamean[1],itau2);                                 
for (i in 2:(n+1)) {                                                                    
thetamean[i] <- mu + phi*(theta[i-1]-mu);                                   
theta[i]~dnorm(thetamean[i],itau2);}                                 
phi1 ~ dbeta(20,1.5);     
phi <- 2*phi1-1;     
mu ~ dnorm(0,0.04);                                                              
itau2 ~ dgamma(2.5,0.025);                                                      
tau <- sqrt(1/itau2) ;
rho ~ dunif(-1,1)                                                        
}# Data for WinBUGS generated by BAUW,
#a free program by Zhang, Z. and Wang, L. (2006)
list(N=47,T=4,y = structure(Data = c(
-0.0194,
-0.0068,
0.0033,
-0.0012,
-0.0054,
-0.0088,
-0.0055,
-0.0023,
-0.0056,
-0.0002,
-0.0092,
-0.0004,
0.0023,
0.0578,
0.0398,
0.0032,
0.0285,
0.0138,
0.0034,
-0.0002,
0.0257,
0.0711,
0.1241,
0.0913,
0.1021,
0.1873,
0.8003,
-0.0719,
-0.3063,
-0.0374,
0.0727,
0.1644,
0.0869,
0.0241,
-0.0298,
-0.3051,
-0.4306,
-0.2197,
-0.1386,
-0.0791,
-0.1567,
-0.0708,
-0.0455,
-0.0577,
-0.0144,
-0.0005,
-0.0733), .Dim = c(47,4)))

就是那个LOAD inits  怎么亮不了，格式错在哪？

没有初始化参数

list 数据格式错了

给顶起来。。。。。

请问winbugs软件得到参数估计之后，应该如何得到历史波动率呢？也就是条件方差序列？

潜在波动率怎么进行估计？

同问，你做出来了吗