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

weixin_493202632025-04-13 9:43

示例：

复制代码

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对于结局的影响。