R语言网状Meta分析---Meta回归（1）（基于gemtc）

示例：

library(gemtc)
help(package="gemtc")
# Fixed effect meta-regression for heart failure prevention
str(hfPrevention)
regressor <- list(coefficient='shared',
                  variable='secondary',
                  control='control')
model <- mtc.model(hfPrevention,
                   type="regression",
                   regressor=regressor,
                   linearModel="fixed")
qw=mtc.run(model)
summary(qw)

 coeficient='shared' 所有其他处理都有一个相对于对照的系数，该系数对于所有处理都是相同的（“shared”），对于所有处理来说是不同的（“unrelated”），或者在处理之间是可交换的（“exchangeable”）。

数据格式如上： treatment包括control和statin两种。secondary是一种协变量，包括0和1。

• 迭代次数：5001 至 25000 次（共 20000 次有效迭代）。
• 链数：4 条链。
• 每条链的样本量：20000。
• 稀疏间隔：1（即每次迭代都保存结果）。
• d.control.statin：
  • 对数优势比估计值为 -0.2584，95% 置信区间为 [-0.3506, -0.16778]，显著为负。
  • 表明 d.control.statin 的增加会降低因变量的对数优势。
• B：
  • 对数优势比估计值为 -0.2052，95% 置信区间为 [-0.4203, 0.01118]，包含 0。
  • 表明 B 的影响可能不显著。

说明，secondary不影响treatment对于结局的影响。