双变量模型和 'predict.rma()' 函数('metafor' 包)

Bivariate model and 'predict.rma()' function ('metafor' package)

假设我有以下模型结果:

> 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_eNaddl 调节器的预测函数绘制回归线:

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.

我想这是因为我没有指定 predict() 函数应该考虑哪个主持人。请注意,上面的函数适用于单变量模型,例如仅 MAT_e。

当您有两个预测变量时,您需要为这两个预测变量指定值。所以,例如:

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

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

因此您可以保持一个预测变量不变,同时改变另一个。