OPENBUGS的SV-T模型程序运行不起来

0关注
4粉丝

已卖：583份资源

博士生

48%

还不是VIP/贵宾

-

0%

威望: 0 级
论坛币: 919 个
通用积分: 2.3936
学术水平: 25 点
热心指数: 22 点
信用等级: 11 点
经验: 8425 点
帖子: 140
精华: 0
在线时间: 414 小时
注册时间: 2013-10-5
最后登录: 2024-1-5

楼主

sheepyang1992

发表于 2015-5-30 21:25:45 |AI写论文

10论坛币

程序如下按了gen inits出现unable to generate initial values for node

#model specification
model
{
for(i in 1:387)
{y[i]~dt(0,p[i],omega)
p[i]<-exp(-theta[i])
}
for(j in 2:387)
{theta[j]~dnorm(theta2[j],itau2)
theta2[j]<-mu+phi*(theta[j-1]-mu)}
theta[1]~dnorm(mu,itau2)
phi<-2*phi1-1
tau<-sqrt(1/itau2)
mu~dnorm(0,0.01)
phi1~dbeta(20,1.5)
omega~dchisqr(8)
itau2~dgamma(2.5,0.025)
}
list(
mu=0,phi1=0.975,itau2=0.02)
#Data
list(y = c(0.000137324,

复制代码

最佳答案

周某仁查看完整内容

谢谢楼主的分享！

分享0 收藏2 回帖

关键词：openbugs SV-T Open BUGS bug 程序模型 openbugs

相关帖子

沙发

周某仁 发表于 2015-5-30 21:25:46

谢谢楼主的分享！

藤椅

sheepyang1992

发表于 2015-5-30 21:48:24

没有人知道吗

板凳

sheepyang1992

发表于 2015-5-31 15:47:39

谁有WINBUGS或者OPENBUGS的帮我跑一下这个程序看跑得出来吗。。我怀疑我的应用有问题。。

#model specification
model
{
for(i in 1:387)
{y[i]~dt(0,p[i],omega)
p[i]<-exp(-theta[i])
}
theta[1]~dnorm(mu,itau2)
for(j in 2:387)
{theta[j]~dnorm(theta2[j],itau2)
theta2[j]<-mu+phi*(theta[j-1]-mu)}
phi<-2*phi1-1
tau<-sqrt(1/itau2)
mu~dnorm(0,0.01)
itau2~dgamma(2.5,0.025)
phi1~dbeta(20,1.5)
omega~dchisqr(8)
}
#Initial values
list(mu=0,itau2=0.975,phi1=0.02,omega=7)
#Data
list(y = c(0.000137324,
-0.054928576,
-0.03238398,
0.044164173,
-0.016604884,
0.018570836,
0.030838062,
0.018829282,
0.022701908,
-0.057055244,
-0.08430086,
-0.09716902,
0.062557457,
-0.022551668,
-0.028610047,
-0.004989212,
-0.011104211,
-0.081817905,
-0.042817379,
-0.058695963,
-0.038120917,
0.013445692,
-0.121048197,
0.139447046,
0.037340599,
-0.021794857,
0.002967788,
-0.042597164,
-0.011508227,
-0.030663261,
-0.148979376,
-0.01300246,
-0.029861524,
-0.028992566,
0.067605335,
-0.027778184,
0.030738829,
-0.022334048,
-0.072560177,
-0.061372055,
-0.01869144,
-0.00327323,
-0.084801926,
-0.05736115,
-0.002204701,
0.100196586,
-0.136768972,
-0.035575406,
-0.048297705,
-0.062137287,
0.010884618,
0.128037732,
-0.008620242,
-0.051165291,
0.075875366,
-0.032442788,
0.032343708,
-0.086327921,
-0.016725466,
0.045127058,
0.025695192,
0.018362535,
0.09142728,
0.062014123,
-0.025888179,
-0.082282318,
0.051537881,
-0.029693102,
0.069071471,
0.040108221,
0.018915007,
0.010053518,
0.024135346,
-0.022353143,
0.011765892,
0.058050243,
0.007440875,
-0.018181417,
0.013509353,
0.044917163,
-0.003685215,
0.048631176,
0.016430896,
0.053252129,
0.00824215,
0.024053815,
0.055744551,
0.011632254,
-0.045377385,
-0.067791175,
-0.028698285,
-0.034386515,
0.000321549,
0.043818747,
-0.009112262,
-0.042695371,
-0.021149647,
0.046498101,
0.022047816,
0.043139381,
-0.036703171,
0.054621694,
0.007434276,
0.037165607,
-0.066254644,
0.068877887,
-0.021242778,
-0.041957292,
0.008779895,
0.042318447,
-0.025071063,
0.008769322,
-0.030086926,
-0.045546523,
-0.016830422,
0.026431955,
0.01114002,
-0.006865061,
-0.005840065,
0.017872066,
-0.002620985,
0.031602841,
-0.004001078,
-0.004796325,
-0.048018332,
-0.038585964,
-0.065585769,
0.003064064,
-0.042850035,
0.027581821,
-0.039234419,
0.00638234,
-0.022317753,
0.015633831,
-0.068880397,
0.036272318,
-0.019060103,
0.059165022,
0.025136028,
0.007889178,
-0.019635623,
0.013568481,
-0.012019829,
0.016957831,
0.002940626,
-0.02452469,
-0.00275132,
0.024437059,
0.030804749,
0.081454488,
0.001305035,
0.001276479,
0.049339284,
-0.047126108,
-0.03298557,
-0.005857381,
-0.010244869,
-0.000489138,
0.018379601,
-0.02045141,
-0.009597397,
0.010880452,
-0.01685966,
-0.027623006,
0.01370165,
0.016647508,
0.010084884,
0.025305517,
-0.007348151,
0.021904818,
-0.002896476,
-0.00921473,
0.02410435,
-0.003498613,
0.020879702,
0.006746126,
-0.013202525,
-0.033440974,
-0.016491006,
0.00249001,
-0.004387832,
-0.053352834,
0.006645887,
-0.008422404,
-0.023307138,
0.03837525,
0.004776989,
0.013823903,
0.007974496,
-0.017664689,
-0.02524017,
-0.028270606,
-0.012740637,
-0.022939926,
0.03024781,
-0.03264972,
-0.012148903,
-0.006188663,
-0.02001084,
-0.030859768,
0.030123829,
-0.048065116,
0.065207997,
0.02194542,
-0.018849236,
-0.026348908,
-0.015152121,
-0.008251679,
-0.01941493,
-0.039841308,
-0.009057272,
-0.002434041,
-0.016517253,
0.036841985,
0.032669382,
0.004852123,
0.0092176,
0.002208463,
0.034380399,
0.00859541,
-0.008665095,
-0.014334914,
-0.023222225,
-0.037620999,
0.019154336,
0.022552722,
0.02001737,
-0.004388766,
0.022973869,
-0.023533216,
-0.02129419,
-0.004689977,
0.016949666,
-0.039529174,
0.011071752,
-0.02012884,
-0.015803963,
-0.000831646,
-0.017095441,
-0.007922825,
-0.018560592,
0.00189602,
0.01674296,
-0.025175828,
-0.010834455,
-0.021539042,
0.038440478,
-0.001839304,
-0.046826559,
0.028925929,
0.008952363,
0.011041327,
-0.029607528,
0.024307599,
-0.022924379,
-0.02661404,
0.006259128,
-0.02358687,
0.040417089,
0.042181635,
0.00125002,
0.03645173,
0.019396472,
-0.01504014,
0.032489201,
-0.011184114,
0.054243157,
0.005515925,
-0.049831654,
0.0194072,
-0.017486102,
-0.017494426,
0.021656355,
-0.040164,
-0.005078554,
-0.008352805,
0.017010713,
-0.030179451,
0.012588539,
0.018566093,
0.015913089,
0.002476267,
0.005255922,
-0.039765939,
-0.022347452,
-0.042011973,
-0.046342473,
0.014042941,
0.015959061,
-0.023234365,
0.009092107,
0.00919252,
0.011181842,
0.007867655,
-0.005327322,
0.019693406,
0.019635538,
0.043985784,
-0.020041004,
-0.014623826,
0.006754818,
0.024294716,
-0.015545561,
-0.028115397,
0.007752484,
-0.02041103,
0.01400319,
0.027955221,
0.010921845,
0.007452459,
-0.018515455,
-0.052001269,
0.007864275,
-0.008656035,
-0.034101204,
-0.004156044,
0.024359844,
-0.010427082,
0.005601376,
0.034303362,
-0.001021388,
-0.027526984,
0.000782653,
-0.026376075,
0.021363314,
-0.002890451,
0.008350168,
0.034237608,
-0.015510127,
-0.029622873,
-0.005001389,
-0.007539355,
0.007608442,
0.003974327,
0.002277984,
-0.00454639,
0.019880284,
-0.021502321,
0.004843506,
0.011162552,
-0.006044345,
0.005898659,
0.032274738,
0.027223956,
0.004169213,
0.01461181,
0.006303266,
-0.010592269,
0.048089765,
0.002369924,
-0.001072639,
0.007812456,
0.006855462,
0.004503628,
-0.014148659,
-0.01675514,
0.049941981,
-0.000830862,
0.024771585,
0.003210082,
0.075883213,
0.090734007,
0.000176997,
0.056385525,
0.015639781,
0.02411776,
0.015561449,
0.027348216,
-0.007354089,
-0.043102501,
-0.042782288,
0.040746177,
0.013356806,
0.019335045,
-0.021098273,
0.039835327,
0.069957547,
0.020191096,
0.045760284,
0.043150483,
0.060821881,
0.024512257,
0.010856508,
-0.054532872,
0.024140857,
0.07786638,
0.022782553))

复制代码

报纸

sheepyang1992

发表于 2015-5-31 20:08:56

我的问题解决啦用winbugs就跑出来了希望帮到需要的同学！64位电脑下载压缩包网址如下就可以了
http://www.mrc-bsu.cam.ac.uk/software/bugs/the-bugs-project-winbugs/然后按照下载网页上的key（免费）按照指示操作，不然软件不能正常使用

[问答] OPENBUGS的SV-T模型程序运行不起来 [推广有奖]

最佳答案

相关帖子

浏览过的帖子

浏览过的版块

本版微信群