亚组差异的元回归

Question

我按照这里的示例进行了元回归。这是我使用的适合我的数据集的代码：

拟合元回归模型以检验亚组差异

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

我做错了什么？

谢谢！

Answer 1

使用

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以讨论差异。