之前看到一个关于winbugs中DIC值的获取问题,我也同样遇到了这样的问题。
DevianceInterface.DICAll [0000058EH]
.dBar REAL 5.241754002842416E-308
.dBarTotal REAL 5.241669125272511E-308
.dHat REAL 1.236473820115811E-303
.dHatTotal REAL 8.649545997602619E-317
.dicTotal REAL 7.98293449523214E+182
.f TextMappers.Formatter Fields
.i INTEGER 2284212
.monitors POINTER [010B5060H]
.num INTEGER 1
.pDTotal REAL 4.244199150914158E-313
.res INTEGER 0
我害怕自己建的模型有错,就用winbugs软件中自带的例子进行验证,验证了好多个,都是出现的这种Trap的文本形式。没有办法获得DIC值。不知道是什么原因?我想软件中自带的算例不该都会有问题的。请高手指点?
2,第二个问题,有的文章中推荐设置链的时候可以多设置几条,当我把链设为两条或者三条时,在对模型初始值设置那一步就通不过,经常出现说gamma分布的形状参数设置过小。有的时候改了的确能通过,但是模拟出的结果,并不太好(比如seeds那个算例,其他算例也有这样的情况),甚至有的模型出现模拟迭代按钮UPdate 按钮拒绝运行,不知道是怎么回事?
本人正在写论文,非常的着急,请高手给予指点,非常的感谢
3我现在想建立一个风险模型,通常我们认为风险服从risk~dlnorm(mu,tau),risk<-p*s,(不知道是否可以放在同一个模型中,虽然在winbugs如何使用的文章中也说过可以,但是我至今没发现有使用过的,这里的risk代表风险,p代表风险概率,s代表事故后果)