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