[R-meta] Questions about interpreting the intercept in a rma.mv model

[Slein, Margaret]

Hi Reza, 
Thank you for the speedy reply, I really appreciate it. I have tried running your code from the second link you sent, but I have no luck in running the emmprep function --> am I missing something?

#install the dev version of metafor
#install.packages("emmeans")
#install.packages("remotes")
#remotes::install_github("wviechtb/metafor")
#remotes::install_github("rvlenth/emmeans", dependencies = TRUE, build_opts = "")

library(metafor)
library(emmeans)

trial <- rma.mv(yi, vi, mods = ~region + impact, method="ML", random = ~1|study_id, data=heatdome_ROM, dfs = "contain")
grd <- metafor::emmprep(trial) #this does not run
grd <- emmeans::emmprep(trial) #this does not run
grd <- emmprep(trial) #this does not run
emmeans(grd, ~region+impact)

Any assistance with this would be greatly appreciated.

Cheers, 
Maggie


On 5/19/23, 6:39 PM, "Reza Norouzian" <rnorouzian using gmail.com <mailto:rnorouzian using gmail.com>> wrote:


[CAUTION: Non-UBC Email]


Dear Maggie,


This question has come up with some frequency on the list. For
example, in the archives
(https://stat.ethz.ch/pipermail/r-sig-meta-analysis/2023-April/subject.html <https://stat.ethz.ch/pipermail/r-sig-meta-analysis/2023-April/subject.html>),
I found the following post:
https://stat.ethz.ch/pipermail/r-sig-meta-analysis/2023-April/004552.html <https://stat.ethz.ch/pipermail/r-sig-meta-analysis/2023-April/004552.html>
which, I believe, exactly answers your questions.


Kind regards,
Reza






On Fri, May 19, 2023 at 6:49 PM Slein, Margaret via
R-sig-meta-analysis <r-sig-meta-analysis using r-project.org <mailto:r-sig-meta-analysis using r-project.org>> wrote:
>
> Hello,
>
> Currently working on a meta-analysis and wondering about how to interpret the intercept in a rma.mv output with lnRR as the effect size and several binary categorical moderators? When I put in the -1 to remove the intercept, it doesn’t provide the full estimates for levels across all categorical moderators? How do I get the model output to provide estimates for all levels of all categorical moderators – is there an additional argument for this?
>
> Any help would be greatly appreciated.
>
> Cheers,
> Maggie
>
> [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-meta-analysis mailing list @ R-sig-meta-analysis using r-project.org <mailto:R-sig-meta-analysis using r-project.org>
> To manage your subscription to this mailing list, go to:
> https://stat.ethz.ch/mailman/listinfo/r-sig-meta-analysis <https://stat.ethz.ch/mailman/listinfo/r-sig-meta-analysis>