指定固定效果两级变量rma.mv

Question

我在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))

在函数的帮助部分有一些解释如何处理这个但我认为我缺乏统计培训使我无法使其工作。当你有我的情况时，它似乎也有所不同。

Answer 1

您只需将其添加到通过mods指定的变量。因此，假设该变量被称为MethodOfDetection，那么它将是mods = ~ PeriSam + MethodOfDetection。