metafor 中 rma 对象的样本外预测

Question

估计回归模型后，通常会提取预测值。但我不知道如何在metafor::rma(中做到这一点

library(metafor)

res <- rma(measure="RR", ai=tpos, bi=tneg, ci=cpos, di=cneg,
       mods = ~ ablat + year, 
       data=dat.bcg)

predict(res, 
    newdata = expand.grid(
      year = 1980,
      ablat = 30:55
      )
    )

它返回 13 个拟合值（数据中的行用于估计 rma 对象，而不是 expand.grid( 对象中的 25 行。

如何对新的data.frame 进行样本外预测？

Answer 1

?predict.rma 的帮助文件将参数指定为newmods 而不是newdata，它似乎需要矩阵而不是data.frame。这应该工作

predict(res, 
        newmods = as.matrix(expand.grid(
          ablat = 30:55,
          year = 1980
        ))
)