全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 winbugs及其他软件专版
4332 7
2014-01-28
各位同仁,正在用Winbugs做课题。写好的程序如下,但是对数据格式始终不知道该怎么写?矩阵格式和矢量格式到底有什么区别?有没有介绍数据格式的参考文献推荐一下。另外,因变量newsp存在有缺失值,有没有办法处理?

我试着一个个变量加载入模型,虽然可以compile和自动加载初始值,但是迭代开始后离开死机,不知道为什么?还请高人帮忙。

另外,有在北京的同仁吗?如果方便的话,想当面请教。

model
{
for (i in 1:m)
  {
   for (t in 1:T)
    {
#Binomial likelihood for observed counts
newsp[i,t]~dbin(p[i,t],nmzrk[i,t])
logit(p[i,t])<-a[1]+a[2]*gdp[i,t]+a[3]*sickbed[i,t]+a[4]*nmzrk[i,t]+a[5]*ncjmrjcsr[i,t]+a[6]*ncrk[i,t]
nmzrk[i,t]~dnorm(0.0001,0.0001)
ncjmrjcsr[i,t]~dnorm(0.0001,0.0001)
ncrk[i,t]~dnorm(0.0001,0.0001)
gdp[i,t]~dnorm(0.0001,0.0001)
sickbed[i,t]~dnorm(0.00001,0.00001)
}
}
#Other priors
for (i in 1:6){
a[i]~dnorm(0,1.0E-6)
}
}
二维码

扫码加我 拉你入群

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

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

全部回复
2014-1-28 23:00:04
i代表不同的地区,共计2065个地区;T代表时间,共6年。样本量比较大,共有6个变量,12390条记录。数据量比较大,所以数据格式不知道该怎么整理。有没有什么软件或较好的方法整理?
二维码

扫码加我 拉你入群

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

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

2014-1-28 23:23:58
As stated in the Winbugs manual (http://www.mrc-bsu.cam.ac.uk/bugs/winbugs/manual14.pdf) data can be entered in S-plus format or in rectangular format. Rectangular format has the form:

Y[] x[]
5 4.1
3 2.1
4 3.5
END
S-plus format (as an example):

list(Y=c(5,3,4),x=c(4.1,2.1,3.5))
二维码

扫码加我 拉你入群

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

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

2014-1-29 00:24:44
在模型中,i代表地区,t代表时间。如果按照您所说的数据格式,是否需要将相同地区按照时间顺序把变量先排一下序?如果不排序,是否会影响结果的解释?最后我还想得到不同地区的率是否的差异?有电话吗?明天给你打电话。
二维码

扫码加我 拉你入群

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

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

2014-1-29 00:29:39

按照你的建议,我从新修改数据格式后,在Load data后,提示 expceted key work structure.什么原因呢?


list(gdp = c(153.5,
170.2695,
194.2984,
216.4595,
271.1755,
311.9047,
85.4,
99.9269,
121.5332,
130.2346,
131.441,
147.9615,
57,
65.0021,
72.0849,
二维码

扫码加我 拉你入群

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

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

2014-1-29 03:30:43
提示: 作者被禁止或删除 内容自动屏蔽
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群