Двумерная модель и функция "Предсказание.рма ()" (пакет "Метафор")

Question

Предположим, у меня есть следующий результат модели:

> summary(msa_res@objects[[1]])

Model Results:

         estimate      se     zval    pval    ci.lb    ci.ub     
intrcpt    0.9397  0.0667  14.0850  <.0001   0.8090   1.0705  ***
MAT_e     -0.0079  0.0035  -2.2691  0.0233  -0.0147  -0.0011    *
Naddl     -0.0385  0.0133  -2.9005  0.0037  -0.0645  -0.0125   **

Я хочу построить линии регрессии, используя функцию прогнозирования для MAT_e а также Naddl модераторов:

preds_MAT_e <- predict(msa_res@objects[[1]], newmods=c(-5:30))
preds_Naddl <- predict(msa_res@objects[[1]], newmods=c(1:6))

Но я получаю этот тип ошибки:

Error in predict.rma(msa_res@objects[[1]], newmods = c(-5:30)) : 
  Dimensions of 'newmods' do not match dimensions of the model.

Я предполагаю, что это потому, что я не указываю, какой модератор должен рассматриваться функцией предиката (). Обратите внимание, что вышеуказанная функция прекрасно работает для одномерных моделей, например, только с MAT_e.

Answer 1

Если у вас есть два предиктора, вам нужно указать значения для обоих предикторов. Так, например:

predict(msa_res@objects[[1]], newmods=cbind(-5:30, 1))

или же

predict(msa_res@objects[[1]], newmods=cbind(10, 1:6))

Таким образом, вы можете держать одну предикторную константу, изменяя другую.