winbugs计算5种SV模型的程序

求救各位大神！！！
我想对BDI指数建立SV模型。用的是楼主的程序，但是Data和INITIAL VALUES就有点摸不着头脑了。
现在Data这里是：list(N=538,y = c(-0.000680522,0.000476216,-0.002995338,……），y用的是BDI的收益率。
INITIAL VALUES是list(itau2 = 7, mu = 0, phi1= 0.975)。
模型正确，成功导入数据，也compile了，gen inits了。
可是到了inference的sample monitor tool，不能设置node，提示：monitor could not be set。
麻烦各位大神告知，这是怎么回事！

求救各位大神！！！
我想对BDI指数建立SV模型。用的是楼主的程序，但是Data和INITIAL VALUES就有点摸不着头脑了。
现在Data这里是：list(N=538,y = c(-0.000680522,0.000476216,-0.002995338,……），y用的是BDI的收益率。
INITIAL VALUES是list(itau2 = 7, mu = 0, phi1= 0.975)。
模型正确，成功导入数据，也compile了，gen inits了。
可是到了inference的sample monitor tool，不能设置node，提示：monitor could not be set。
麻烦各位大神告知，这是怎么回事！

陈金海 发表于 2013-12-18 20:59
求救各位大神！！！
我想对BDI指数建立SV模型。用的是楼主的程序，但是Data和INITIAL VALUES就有点摸不着头 ...

http://cos.name/2008/12/statistical-analysis-and-winbugs-part-1/
这里有入门，代码有些小错误看楼下评论
我也在做SV模型 希望可以交流下 QQ：2360934263

这段程序怎么运行不了？？求解

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)))

太贵了！！！

稀缺资源，所需的论坛币较多，希望能够有所帮助

我给你钱行不，这个论坛币，他妈的比卖肾还难搞。

creditperson 发表于 2008-10-7 00:23
一分钱一分货，这是参考很多文献编写的程序，并且我对每一句代码都做了仔细检查，我也用这些程序作了很多次 ...

你好我想请问一下，程序里面报告数据的载入代码和初始值赋值代码吗？

请问一下程序里面包不包括数据的载入和初始值赋值呀？

很有用~