关于winbugs中运行sv模型中遇到问题，请高手指点

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

在winbugs中运行sv-n模型时，
程序如下，在数据导入之后，按compile按钮时，显示array index is grater than array upper boundy for y ，如果数据减少至20程序就可以运行，请高手指点！万分感谢{:soso_e100:}
model{
### likelihood: joint distribution of y
for(i in 1:n)
{p[i]<-1/exp(theta[i])
y[i]~dnorm(0,p[i])
}
### prior distributions*
phi1~dbeta(20,1.5)
phi<-2*phi1-1
mu~dnorm(0,0.01)
itau2~dgamma(2.5,0.025)
tau<-sqrt(1/itau2)
theta0~dnorm(mu,itau2)
thmean[1]<-mu+phi*(theta0-mu)
theta[1]~dnorm(thmean[1],itau2)
for(j in 2:n)
{thmean[j]<-mu+phi*(theta[j-1]-mu)
theta[j]~dnorm(thmean[j],itau2)
}
}

# DATA
list(n=76,y=c(-0.2352,
-3.7769,
-3.5341,
-1.8623,
4.3202,
2.8627,
0.9487,
1.1008,
-1.7228,
0.4690,
-0.9836,
-1.7980,
0.1213,
-0.8333,
0.0567,
-0.3118,
-1.7170,
-1.8562,
-8.6653,
0.1957,
-1.3342,
3.1475,
-1.6588,
-0.5408,
2.7101,
2.3157,
-2.4390,
0.1674,
-4.3934,
-0.5687,
0.9867,
-0.8954,
-1.4493,
3.3308,
-0.8124,
-1.8262,
-0.3651,
1.5883,
-1.3877,
2.5499,
0.2302,
4.4972,
-3.1278,
-2.4211,
-0.6497,
-0.5901,
-1.5454,
0.9628,
-2.5905,
-0.8705,
0.1471,
-1.2321,
0.5959,
-0.0432,
2.2637,
-0.5950,
-2.1195,
-0.5927,
-0.3382,
0.2680,
0.7632,
-1.3312,
-0.1155,
2.1279,
-1.0219,
-0.8355,
-0.9967,
0.3234,
0.1993,
0.2141,
0.0376,
1.4259,
-1.6190,
-0.0172,
1.6469,
3.0316))
Initial value list(phi1=0.975, mu=0, itau2=50)

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：winbugs WINBUG 高手指点 SV模型 BUGS 模型 数据 程序

程序无误.
Inference for Bugs model at "model.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.941 0.204   0.547   0.819   0.940   1.055   1.356
phi1       0.900 0.044   0.810   0.869   0.901   0.932   0.978
tau        0.284 0.043   0.201   0.254   0.288   0.314   0.363
deviance 314.186 3.745 307.600 311.400 314.100 317.100 320.900

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

恩，难道是安装的程序有问题？我再从新装一下，试试，谢谢啦

hizhangzhang 发表于 2012-5-21 08:25
恩，难道是安装的程序有问题？我再从新装一下，试试，谢谢啦

试一下底下程序
若有错误
请把错误信息完整贴上来
modelsv.rar (733 Bytes) 本附件包括：

• modelsv.bug

2012-5-21 18:40:37 上传

epoh 发表于 2012-5-21 18:40
试一下底下程序
若有错误
请把错误信息完整贴上来

操作过程问题.doc (57 KB)

2012-5-21 21:34:40 上传

epoh 发表于 2012-5-21 18:40
试一下底下程序
若有错误
请把错误信息完整贴上来

已经上传了，不知道是不是操作问题了

hizhangzhang 发表于 2012-5-21 21:35

pls Obtaining the key(Keys.ocf) for unrestricted use

http://www.mrc-bsu.cam.ac.uk/bugs/winbugs/contents.shtml

万分感谢，操作了一下，可以了，谢谢啊！如果可以像好好谢谢您！

epoh 发表于 2012-5-21 22:03
pls Obtaining the key(Keys.ocf) for unrestricted use

http://www.mrc-bsu.cam.ac.uk/bugs/winbugs/ ...

万分感谢，已经可以运行了，还有一点像请教一下，您是否知道拟合过程中的残差数据可以在哪里看到？谢谢啦

小女子 弱弱问一句，我已经估计出SV模型的参数，怎么得到残差序列啊