第三节　效应指标为比值比的二分类数据贝叶斯Meta分析

［例3.1］以结核病数据为例，以OR为效应量，基于二项-正态分布层次模型分析框架，说明用R2jags包进行随机效应模型贝叶斯Meta分析的具体过程。

［解］第一步，设置工作目录和种子数、加载包，定义模型。

第二步，以data.farme（）函数建立一个名dat的数据集，用于后面分析OR、RR、RD等不同效应量使用。

第三步，以list（）函数建立一个适用R2jags包分析的列表数据集，命名为bcgdat；设置初始值和要监控的参数；以jags（）拟合模型；并以print（）函数显示结果。

计算结果如下：

结果解读：由感兴趣的参数的相应R=1.001可知，马尔可夫链已收敛；研究间异质性标准差后验点估计及95%CI为0.688（0.389，1.156）；OR后验均数点估计及95%CI为0.480（0.293，0.723）。