我正在使用 R 中 metafor 包中的 rma.mv 进行荟萃分析。我有效应大小 ( yi
)、方差 ( vi
)、一个调节因子和两个随机效应变量。最近有人指出我应该将变量Method of detection
(两个级别)作为固定效应变量包含在模型中。所以这就是问题所在。我不确定如何将此变量指定为固定。
这是我的代码:
overall <- rma.mv(yi, vi, data = ef.s,
method="REML",
mods = ~ PeriSam,
random = list(~1|species,
~1|study),
R=list(species=vcv1))
在函数的帮助部分有一些解释如何解决这个问题,但我认为我缺乏统计培训让我无法让它发挥作用。此外,当你有主持人时,它似乎也会有所不同,就像我的情况一样。