Прогнозирование вне образца для объекта 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
        ))
)