WinBUGS在统计分析中的应用（第三部分）

[未知用户] 这里也有个Andrew Gelman，是哥伦比亚大学的，不知是不是一个人。此人的博客很是精彩。

[未知用户] 那个是Andrew Thomas，不是一个人，不过Andrew Gelman写的几本书很有名啊。他的博客好像有多个作者参与呢，Aleks Jakulin 在March 28的一篇帖子GPU Supercomputers很有启发性啊。随便写个贴就是一个博士生课题啊。牛！

[未知用户] 我只看了些此公关于可视化的文章，觉得很不错，不过我发现这人比较固执，有些想法很个人主义。

[未知用户] 哈哈，我同意太云。此公经常抛出一些惊人之语，比如上次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，那就爽了。当然只是一厢情愿，估计会有很多技术细节：）

我正在用WinBUGS做mcmc。想请问个问题：WinBUGS中反三角函数和圆周率怎么表示啊？？

大家好

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

alpha0 ~ dflat()这个dflat在density中也没有啊,这个该怎么解决.
非常感谢!!

请教一下，如何用winbugs来编一些复杂的分布函数？

[未知用户] 如果你用word的话，直接ctrl-c + ctrl-v就行了。如果是LaTex的话，最好存储成.ps文件。不过我通常都将结果输出在R，然后利用R强大的图形功能作图。

[未知用户] 首先，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文件也未尝不是一件一劳永逸的事情。