[R-meta] Plotting interaction effect from a rma.mv() fit

[Viechtbauer, Wolfgang (SP)]

Dear Simon,

Not aware of any packages that interface with rma.mv, but this is relatively easy to do by hand. Just use predict() with the newmods argument to compute the predicted values as a function of 'time' for each level of 'outcome' and then create a plot with those lines.

Best,
Wolfgang

>-----Original Message-----
>From: R-sig-meta-analysis [mailto:r-sig-meta-analysis-bounces using r-project.org] On
>Behalf Of Simon Harmel
>Sent: Tuesday, 16 March, 2021 5:25
>To: R meta
>Subject: [R-meta] Plotting interaction effect from a rma.mv() fit
>
>Dear All,
>
>Are there plotting packages (e.g., `effects::plot(allEffects())`,
>`sjPlot::plot_model()`, `emmeans::emmip()`  ) to plot the predictions of a
>say rma.mv() model?
>
>For example, is it possible to plot the interaction effect below?
>
>library(metafor)
>dat <- read.csv("https://raw.githubusercontent.com/hkil/m/master/tst.csv")
>fit <- rma.mv(d ~ factor(outcome)*time, V = SE^2, random= ~1|id, data = dat)