[R-meta] Four-Level Meta-Analysis in metafor

[Can Evenkamp]

Hello all,
I am an undergraduate student conducting my first meta-analysis.
My data consists of 187 effect sizes (Pearson’s r) nested in 87 samples nested in 34 studies.
Therefore, I used the following four-level model (LRT indicated better fit than a three-level model):

four.level.model <- rma.mv(yi = fisherz,
                            V = vi,
                            data = data,
                            random = ~ 1 | StudyID/SampleID/EffectID,
                            test = "t",
                            method = "REML")

My question is whether and how I can create a forest plot aggregated at the study level?
I have seen approaches for three-level models (https://doi.org/10.5964/meth.4013), but I am unsure how to apply this logic for four levels. Generally, I find very little information on four-level models (seems to be rather uncommon).
In the end, I would like to have a forest plot showing the aggregated effect sizes (+ confidence intervals) for each study and the pooled estimate. Preferably, these should be back-transformed to Pearson’s r (my model is based on Fisher’s z metric). Is this possible with my model?

I appreciate any advice or suggestions you may have.
Please let me know if you need any additional information or details about my data or model.
Thank you very much for your time and support!

Best regards,
Can