我想在R中使用metafor进行多级元分析 我有嵌套在出版物(“studyid”)中的样本(“sampleid”)嵌套的效果大小(“esid”)。 我有四个子组(“task.type”)。
每个子组的平均效果大小取决于我是使用task.type作为主持人还是单独为每个子组运行rma.mv命令。
这是我使用task.type作为主持人的代码:
multi.task <- rma.mv(yi=g, V = var.g, data=df, random=list(~ 1 | esid, ~1 | sampleid, ~1 | studyid), mods=~factor(task.type)-1)
这是我在为每个任务计算模型时使用的。单独输入:
task.X <- rma.mv(yi=g, V = var.g, data=df, subset=(task=="X"), random=list(~ 1 | esid, ~ 1 |sampleid, ~ 1 | studyid))
使用第一种方法(task.type作为主持人)与第二种方法(每种task.type的rma.mv)的四种任务类型的含义:
为什么结果如此不同?