关于SVCJ模型的winBUGS完整代码

谢谢楼主

谢谢分享~

很感谢楼主的分享。我在应用这个程序时，MODEL检验没有通过，可能是什么问题？有编译成功的吗？

qyqx5202000 发表于 2011-3-16 09:39
我自己加了list和inits，MODEL检验和上传数据都行，就是到了load inits那里总是提示有uninitialized variab ...

可以把您上传数据的语句进行分享吗？另外，您的参数赋初值是不是掉了一个muy?请指教~

我把此程序中的y-y[i-1]改成了r,然后在程序的后面加上了以下语句做测试：
#library(R2WinBUGS)
#data
list(
r= c(-0.008600000, 0.007400000, -0.01240000, -0.02860000, -0.009900000, 0.01080000, 0.005800000, -0.02260000, -0.01140000, -0.01800000, 0.0009000000, -0.007600000
),
N = 12
)
#inits
#N = length(y)
#ksy0=rep(0,N)
#ksv0=rep(0,N)
#J0=rep(0,N)
#V0=rep(0,N)
list(mu = 1, k = 1, kt = 1, tauv =1, rho = 1, lambda = 1, tauy=1,muy=1,
muv=1,rhoj=1,
ksy=c(0,0,0,0,0,0,0,0,0,0,0,0
),
ksv=c(0,0,0,0,0,0,0,0,0,0,0,0
),
J=c(0,0,0,0,0,0,0,0,0,0,0,0
),
v = c(1.405700,1.397100,1.404500,1.392100,1.363500,1.353600,1.364400,1.370200,1.347600,1.336200,1.318200,1.319100
)
)
#parameters< -c("mu","k","kt","tauv","rho","lambda","tauy","muy","muv","rhoj")
#jump.sim< -bugs(datax, inits,parameters,"model.txt",n.chains=1,n.iter=1000)
#save(jump.sim,file="jump.Rdata")
但是程序一直提示“undefined variable”是怎么回事？
多谢各位大虾指教啊~

试了很多遍，能编译，但是V值算出来是错误的，原因在于程序中v 的后验分布并不是常规分布，要用模拟退火才能得到。

谢谢。

谢谢啦

qyqx5202000 发表于 2011-3-16 09:39
我自己加了list和inits，MODEL检验和上传数据都行，就是到了load inits那里总是提示有uninitialized variab ...

我也遇到同样的问题  你解决了么

楼主，我想问下，为什么我做complie是总是显示variable n is not defined.