r - 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对象的数据中的行，而不是对象中的 25 行expand.grid(。

我怎样才能对新的样本进行预测data.frame？

Answer 1

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

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