metafor 中 rma 对象的样本外预测

Out of sample prediction for rma object in metafor

估计回归模型后,通常会提取预测值。但我不知道如何在 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
      )
    )

其中returns 13个拟合值(数据中用于估计rma对象的行,而不是expand.grid(对象中的25行。

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

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

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