我正在使用metaforR 进行荟萃分析。我想比较一个因素的 7 个水平（即不同类型的治疗）。

fit <- rma (yi, vi, mods = type_of_treatment - 1, data = dat)
fit

我发现几个网站解释了如何通过使用来比较级别anova()（例如，元回归后的事后分析和http://www.metafor-project.org/doku.php/tips:testing_factors_lincoms?s[]=anova）；但是，因为我有七个级别（即七种不同类型的治疗），我自己做比较是很麻烦的。

我在 multcomp 包中进行了尝试glht()，但以下命令提供了错误。

summary(glht(fit, linfct = mcp(~type_of_treatment = "Tukey"))

如果您能告诉我是否有任何简单的方法可以在级别之间进行成对比较，例如lmglht()或emmeans()lm，我将不胜感激。