急急急！ 用R2winbugs来做SV模型初始值设定问题

用winbugs做SV模型后，DIC无法计算，于是决定采用 R中R2winbugs程序包来做SV模型，但是初始值设定遇到了问题.
比如，变量THETA的初始值，我想要随机的生成，但是他有很多有N个，而且是和thetao是循环随机生成的如：
hetamean[i] <- mu + phi*(theta[i-1]-mu);
theta[i] ~ dnorm(thetamean[i],itau2); }
这样子的，按照参考文献用winbugs作的话，只要点gen inits 就可以了 ，现在我换成R语言应该怎么写呢？难道是要我先编程生成随机数，然后再list theta=c()出来吗？如果是这样的话，程序应该怎么写呢？求大侠指导，急急急！万分感谢！

自顶，千万别沉了，求解答

自顶

问题已经解决了，谢谢各位的关注

梦苏＾.＾ 发表于 2013-5-16 09:55
问题已经解决了，谢谢各位的关注

新的问题，解决有赏：
      用winbugs估计出来的波动率在哪里查看呢？是THETA的平均值吗？还是model 选项下面SAVE STATE 选项下面的state of sampler 中的theta值呢？求大虾解答

已解决

你好！我有问题请教你，可以帮我一下吗？我qq是862415257，有点着急，万分感谢！

你好，请问怎么计算波动率来进行预测

你好。我在用R2winbugs做SV模型的时候也遇到了这个问题。用winbugs作的话，只要点gen inits 就可以了 ，现在我换成R语言应该怎么写？大侠求助。跪谢！！

梦苏＾.＾ 发表于 2013-5-16 09:55
问题已经解决了，谢谢各位的关注

用winbugs做SV模型后，DIC无法计算，于是决定采用 R中R2winbugs程序包来做SV模型，但是初始值设定遇到了问题.
比如，变量THETA的初始值，我想要随机的生成，但是他有很多有N个，而且是和thetao是循环随机生成的如：
hetamean <- mu + phi*(theta[i-1]-mu);
theta ~ dnorm(thetamean,itau2); }
这样子的，按照参考文献用winbugs作的话，只要点gen inits 就可以了 ，现在我换成R语言应该怎么写呢？
这个问题怎么解决 就是初始值没法初始化