metafor 中的元回归 - 亚组差异
Meta-regression in metafor - subgroup differences
我按照 here 中的示例执行了元回归。
这是我使用的适用于我的数据集的代码：
fit meta-regression model to test for subgroup differences
resMeta <- rma(xi=nphy, ti=ni, mods = ~ pop, data=metaAAS)
'metaAAS' 保存来自 csv 文件的数据，'pop' 是版主。
问题是 'pop' 有 4 个级别，描述为 0、1、2 和 3，主持人测试 (QM) 的输出仅显示一个自由度，而不是 3（级别数 - 1):
主持人测试（系数2）：
QM(df = 1) = 8.7150, p-val = 0.0032
我做错了什么？
谢谢！
使用
resMeta <- rma(measure="IR", xi=nphy, ti=ni, mods = ~ factor(pop), data=metaAAS)
或
resMeta <- rma(measure="IR", xi=nphy, ti=ni, mods = ~ 0 + factor(pop), data=metaAAS)
取决于您希望如何参数化模型。有关差异的讨论，请参阅 https://www.metafor-project.org/doku.php/tips:models_with_or_without_intercept。
我按照 here 中的示例执行了元回归。 这是我使用的适用于我的数据集的代码：
fit meta-regression model to test for subgroup differences
resMeta <- rma(xi=nphy, ti=ni, mods = ~ pop, data=metaAAS)
'metaAAS' 保存来自 csv 文件的数据，'pop' 是版主。 问题是 'pop' 有 4 个级别，描述为 0、1、2 和 3，主持人测试 (QM) 的输出仅显示一个自由度，而不是 3（级别数 - 1):
主持人测试（系数2）： QM(df = 1) = 8.7150, p-val = 0.0032
我做错了什么？
谢谢！
使用
resMeta <- rma(measure="IR", xi=nphy, ti=ni, mods = ~ factor(pop), data=metaAAS)
或
resMeta <- rma(measure="IR", xi=nphy, ti=ni, mods = ~ 0 + factor(pop), data=metaAAS)
取决于您希望如何参数化模型。有关差异的讨论，请参阅 https://www.metafor-project.org/doku.php/tips:models_with_or_without_intercept。