Winbugs 数据格式

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

各位同仁，正在用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)
}
}

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：winbugs WINBUG BUGS 数据格式 bug 参考文献 因变量 北京 程序 模型

i代表不同的地区，共计2065个地区；T代表时间，共6年。样本量比较大，共有6个变量，12390条记录。数据量比较大，所以数据格式不知道该怎么整理。有没有什么软件或较好的方法整理？

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))

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

按照你的建议，我从新修改数据格式后，在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,

谢谢您的帮助。我试着联系TimeT以寻求帮助。如果可以的话，也可以合作课题。