• 综合主站
  • WinBUGS在统计分析中的应用(第三部分)

[未知用户] 我只看了些此公关于可视化的文章,觉得很不错,不过我发现这人比较固执,有些想法很个人主义。
[未知用户] 哈哈,我同意太云。此公经常抛出一些惊人之语,比如上次R和SAS之战的时候,他直接讲“... it's good to hear that SAS is in trouble. I just hate SAS.”,额滴神呀,气氛本来就颇有些火药味,他这么一讲,还不炸开了锅…… :mrgreen:
[未知用户] Andrew Gelman 是rubin的学生,他们group是统计学的一大势力。。meng xiaoli就是他们的人
[未知用户] 去年我跟我们公司一技术人员讨论 GPU 时,就琢磨:如果能用 GPU 跑 R,那就爽了。当然只是一厢情愿,估计会有很多技术细节:)
7 天 后
我正在用WinBUGS做mcmc。想请问个问题:WinBUGS中反三角函数和圆周率怎么表示啊??
11 天 后
大家好

请教一下winbugs可以用Doodle来建模,这个例子中的b[1:N] ~ car.normal(adj[], weights[], num[], tau)该怎样用Doodle来建模,还有

alpha0 ~ dflat()这个dflat在density中也没有啊,这个该怎么解决.
非常感谢!!
请教一下,如何用winbugs来编一些复杂的分布函数?
1 个月 后
[未知用户] 如果你用word的话,直接ctrl-c + ctrl-v就行了。如果是LaTex的话,最好存储成.ps文件。不过我通常都将结果输出在R,然后利用R强大的图形功能作图。
[未知用户] 你需要下载Blackbox,即BUGS的开发软件,然后自己编写需要的distribution.具体请参见
http://www.winbugs-development.org.uk/wbdev.html
[未知用户] 首先,car.normal是一个用户定义的分布函数,所以你在BUGS自带的函数中无法找到,自然也就无法将它在doodle中画出来。

其次,dflat()在distribution中有,只不过没有在user manual中列出来罢了。但是,强烈建议不要使用这个distribution,因为他的source code并没有对他的边界有一个很好的定义,以至于在一些复杂的model很容易让程序crash. 几个比较好的替代有: dunif(), dnorm(),dbeta().

dunif(a, b)很简单,只要自己定义边界就可以了。
dnorm(mean,precision)其实最实用,只要把precision设得小一些就可以了,比如0.000001,通常我们都会用这个来代替dunif()。
dbeta(1,1)也可以用来代替dunif().
[未知用户] 反三角函数比较难一些,比如arctan(x)用微分可以表示为darctan(x)/dx=1/(1+x^2),你需要在WBDiff中编写一个微分函数来求x,具体WBDiff的编写可以参考
http://www.winbugs-development.org.uk/wbdiff.html

pi的表示很简单,你在你的data中定义一个pi=3.1415926,然后在model中使用pi就可以了。
小弟现正在剑桥mrc-bsu做postdoc,具体的项目就是BUGS的开发以及在生物及医学方面应用。WinBUGS这个软件是我两个老板 David Spiegehalter,Dave Lunn 和其它一些牛人共同开发的。我们现在正在从WinBUGS 转向openBUGS,目的是将它做成open source的软件以应用在更广的领域。我现在正在开发BUGS中的WBDiff部分并将它应用在二型糖尿病的动态系统的数据分析中。有兴趣或有问题的同学可以和我联系:

chen.wei@mrc-bsu.cam.ac.uk

还有我们这里每年会举办3-4次BUGS的培训,2天的课程,在英国或能到英国出差的同学有兴趣的话可以参加,主讲人是David speigehalter 和 Dave Lunn。
楼主
你好:

还有一个问题麻烦你,GeoBUGS可以做简单的空间数据分析,map tool里面的地图只有很少的几个国家,请问一下,如果想做中国的空间数据,应该怎么办?
非常感谢楼主回答。
GeoBUGS 界面只能完成简单及有限的数据分析,要想做自己的数据分析,最简单的方式是以楼主附录中列出的model作为source code模板,根据自己的模型加以改动,然后load进你的data,进行mcmc采样分析。
[未知用户] 可以从ArcInfo中通过工具将txt格式的地图数据导出,然后在GeoBUGS中导入。
[未知用户] 可以,具体可以了解OpenBUGS的开发文档。不过pascal语言过于古老,实在不是很适用于目前很多开放平台。有兴趣的同志可以研究。目前我公司也在着力研究开发有关内容。基本上会吸收OpenBUGS的一部分代码,并在sampling的计算中引入更加高效的计算机制。
[未知用户] 很好啊,希望学成之后可以为国效力开发国人自己的专业统计软件。
[未知用户] 当然了ArcInfo中的中国地图的shp文件可以在网上搜索。如果你要做Spatial Data Analysis,又不想自己开发软件,目前恐怕最好的方式就是WinBUGS+GeoBUGS+ArcInfo来做了。当然了R中的GeoR Package也可以利用,不过不一定能达到你的模型要求,需要自己再开发。当然了其实GeoBUGS中的地图数据相对简单,就是把地理区块分成多边形,每个多边形有边界上的点的经纬度坐标定义。如果你可以得到中国行政区域的边界经纬度(如通过一些GIS软件),直接用记事本编辑一个文本格式的map文件也未尝不是一件一劳永逸的事情。
[未知用户] 之前给出的disease mapping那本书中的Ohio空气污染致使肺病的区域性分布规律的研究,有用GeoBUGS的完整范例,以及不同模型之间的比较性分析,可以参考。如果是自己的研究数据,如针对中国某省份的某种疾病受某种空间因素的影响,那你就需要找到那个区域的各个子区域的边界经纬度值,然后根据GeoBUSG的文档来创建map文件。除非你的模型非常特殊,BUGS应付不了,不然应该能满足通常的需求。
您好:

本人是初学者,对ArcInfo软件不熟悉,可不可以提供ArcInfo中的中国地图txt格式的文件,或者ArcInfo中的中国地图的shp文件的下载地址,我的毕业论文想做一个关于中国各省车险费率的一个比较,想通过GeoBUGS中的作图工具标出不同省市的费率差别,以期达到更直观的目的,恳请赐教。
我的邮箱是hongaf@163.com,如果可能的话可不可以发给我现成的中国地图的txt文件。
非常非常非常感谢。