winbugs 计算DIC问题求教

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

最近用winbugs 算SV模型族，在得到DIC值，在计算SV-N时，得到DIC值如下，计算的时候应该主要是看Dbar和Dhat，当然也能直接出来DIC，但是我的值是不是很有问题呀 为什么后面是E-308啊，不知道问题在哪，请大牛指导！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏4 回帖

关键词：winbugs WINBUG 问题求教 BUGS dic 模型

哈哈!刚看到你附件的图档.
你的程序根本没跑
就掉进trap

啊？啊？我把odc文件发个您，您能提供个其它的联系方式么，QQ什么的，您能帮我看看么？ 因为参数当中有个初始值我不知道该如何赋初始值， ，所以我采用自动生成初始值，下面显示能生成初始值，模型能够跑，但是跑得很慢，3-4个小时 40000次迭代吧，最后各个参数的结果都能出来，也都收敛，就是DIC有问题，这是为什么呢？因为现在时间比较急，真的非常需要您指导下，谢谢了哦

减少一些迭代和循环，程序执行就快了，
把你的程序传上来大家帮你看看，或许问题就解决了。
如果涉及到知识产权问题，你可以把你的思想和模型说详细一些，也许能解决！
有时候程序中的一些小bugs，如果没有其他人指导，自己或许一辈子也看不出来。

以上是我算sv-t的模型和初始值取值，数据比较多1437个，我是直接贴在y=c（）里面，初始值的theta也如此，但是不知道初始值theta取值有没有问题，总之模型能够运算，但是就是得不到DIC值，得到的后面20000次各个参数的值如下，请大牛们帮帮忙看看我的程序是否有问题，DIC值要如何得到？非常感谢
node  mean         sd       MC error      2.5%           median         97.5%        start        sample
mu   0.6701     0.1612           0.001671    0.3498           0.6709        0.9862        20001        20000
omega 17.55     5.013            0.1664              9.752              16.95       29.2        20001        20000
phi      0.8962     0.02058   5.073E-4   0.8535      0.897            0.9341        20001        20000
tau     0.5872     0.0574     0.001824   0.4763     0.5863      0.701         20001        20000

以上原始数据和theta的初始值省略了一些

伊甸之城 发表于 2011-12-17 13:54
以上原始数据和theta的初始值省略了一些

程序小修
我先试跑程序burn in=10000,n.iter=20000,(约20分钟)
结果如下:
Inference for Bugs model at "sv_t.bug", fit using WinBUGS,
1 chains, each with 20000 iterations (first 10000 discarded)
n.sims = 10000 iterations saved
             mean     sd     2.5%      25%      50%      75%    97.5%
mu          0.105  0.078   -0.051    0.054    0.106    0.157    0.259
phi          0.947  0.010    0.928    0.941    0.948    0.954    0.966
tau         0.139  0.007    0.127    0.134    0.139    0.144    0.153
beta        1.055  0.041    0.975    1.027    1.054    1.082    1.138
k            2.188  0.116    2.014    2.097    2.175    2.260    2.450
deviance 5961.406 15.290 5932.000 5947.000 5966.000 5972.000 5985.000

DIC info (using the rule, pD = Dbar-Dhat)
pD = 7.4 and DIC = 5968.8
DIC is an estimate of expected predictive error (lower deviance is better).

sv_t.bug(need pw)
   sv_t.rar (6.63 KB)

2011-12-17 16:03:49 上传

您好，运行您的那个sv_t.bug文件后，载入初始值时仍然有this chain contains uninitialized variables错误，我选择自动生成初始值，运行后，速度还是一样慢，我估计可能是我笔记本的原因，我这是09年的上半年的配置吧，CPU T6400的，2.0GHZ,条子是2G的DDR3，是不是这个原因导致慢呢？另外我用的是winbugs1.4.3版本

换了一台本子，这回速度快了很多，四万次迭代大概41分钟，跟您跑的速度差不多了，但是问题还是一样，跑完以后，点击DIC里面的set 然后再点DIC值 看到的还是上面贴出来的那个Trap，这是怎么回事呢？是我的软件版本问题？还是软件的设置问题？还是其它原因呢？谢谢您哦

伊甸之城 发表于 2011-12-18 00:29
换了一台本子，这回速度快了很多，四万次迭代大概41分钟，跟您跑的速度差不多了，但是问题还是一样，跑完以 ...

1.通常不會給足initial values,只是重點給出
  所以this chain contains uninitialized variables
  是個提醒不是錯誤
  按完 load inits
  再按 gen inits 就OK

2.我改个方式讓你執行:
1.在c碟,新建文件夹"Bugs",
    作为working.directory
    放进两文件 sv_t.bug,eden.txt(your data)
2.注意你的winbugs安装路径  若不是"D:/WinBUGS14/"
    请依你所安装修改
3.所有结果都存放在c:\bugs,请注意看log.odc
4.在 R 运行程序(eden_script)
5.结果如下:
  Inference for Bugs model at "sv_t.bug", fit using WinBUGS,
1 chains, each with 40000 iterations (first 30000 discarded)
n.sims = 10000 iterations saved
             mean    sd     2.5%      25%      50%      75%    97.5%
mu         0.222 0.060    0.105    0.182    0.221    0.262    0.340
phi         0.928 0.012    0.904    0.920    0.928    0.936    0.950
tau         0.161 0.009    0.142    0.154    0.161    0.168    0.178
beta       1.118 0.034    1.054    1.095    1.117    1.140    1.185
k            2.528 0.157    2.229    2.417    2.525    2.630    2.847
deviance 5841.791 5.774 5830.000 5838.000 5842.000 5846.000 5852.000

DIC info (using the rule, pD = Dbar-Dhat)
pD = 11.8 and DIC = 5853.6
DIC is an estimate of expected predictive error (lower deviance is better).

eden (need pw)
   eden.rar (166.13 KB)

2011-12-18 10:13:46 上传