metafor

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
        ))
)