我正在使用meta for进行多变量荟萃分析。我想预测新主持人的影响,但对于其中一个主持人,我不想输入值本身,但我希望预测函数对该主持人进行荟萃分析(在下面的示例中,这是grade
)
我在这里使用的是玩具数据集dat.bangertdrowns2004
。
# load metafor package
library(metafor)
dat <- dat.bangertdrowns2004
# fit model
res <- rma.mv(yi, vi, mods = ~ grade + length + pers - 1, random = ~ 1 | ni, data=dat)
res
# predict
gridpred <- expand.grid(length=4:10, pers=1)
pred.res <- predict.rma(object = res , newmods = as.matrix(gridpred))
当我运行这段代码时,我得到一个错误,因为它提到我需要输入grade
的值。
Error in predict.rma(object = res, newmods = as.matrix(gridpred)) :
Argument 'newmods' does not specify values for this variable: grade
有没有一种方法可以使用predict.rma()对length=4:10
和pers=1
的成绩进行荟萃分析?
2条答案
按热度按时间xriantvc1#
如果你确实在寻找这样的“调整后的均值/估计值”,那么应该将
grade
的值保持在其均值。所以语法应该是:如果你想使用
emmeans
,你可以用途:j2qf4p5b2#
谢谢你的建议沃尔夫冈我已经根据您发送的链接(https://www.metafor-project.org/doku.php/tips:computing_adjusted_effects)在“附录:使用emmeans包”。
这是所有年级每个长度值的平均值,这是你的意思吗?