R或Python中的层次贝叶斯分析
层次贝叶斯模型通常在市场营销、政治科学和计量经济学中使用。不过,我知道的唯一一个相关的工具包是 bayesm
,它其实是一本书(《贝叶斯统计与市场营销》,作者是Rossi等人)的配套工具。我是不是漏掉了什么?有没有其他的R或Python软件包可以完成这个任务,或者有没有在这些语言中做好的示例?
7 个回答
在 MCMCpack 这个R语言的工具包里,有几种层次模型。根据我的了解,它是许多常见模型类型中速度最快的采样器。(我在里面写了一个[层次项目反应][2]模型。)
[RJAGS][3] 的名字就说明了它的功能。你只需要写一个带有jags风格的.bug模型,提供数据给R,然后从R中调用Jags就可以了。
这里有OpenBUGS和一些R语言的辅助包。你可以去Gelman的网站看看他的书,里面有很多相关的链接:
在Python方面,我只知道PyMC这个工具:
补充:我添加了一个链接,指向Gelman书中适当的附录,里面有使用R和BUGS的例子,可以在线查看。
这里有四本关于层次建模和贝叶斯分析的书,书中都包含了R语言的代码。
《空间数据的层次建模与分析》(统计与应用概率专著)(精装) http://www.amazon.com/gp/product/158488410X
《使用回归和多层/层次模型进行数据分析》(平装) http://www.amazon.com/Analysis-Regression-Multilevel-Hierarchical-Models/dp/052168689X/ref=pd_sim_b_1
《使用R进行贝叶斯计算》(Use R)(平装) http://www.amazon.com/Bayesian-Computation-R-Use/dp/0387922970/ref=pd_bxgy_b_img_c
《环境科学的层次建模:统计方法与应用》(牛津生物学)(平装)(我猜这本书也有R语言的代码,因为两位作者都广泛使用R)
我知道有些关于Python的书会涉及多变量分析(比如《集体智能》),但我还没见过真正深入探讨贝叶斯或层次建模的书。