[R-meta] Metafor modeling question

Yuhang Hu yh342 @end|ng |rom n@u@edu
Sun Mar 26 06:27:40 CEST 2023

Dear All,

I have two questions about the following model:

mm <- rma.mv(y ~ SES, V, random = list(~SES|study, ~SES|outcome), struct =

fit to a dataset structurally similar to:

study outcome SES
1     A       1.5
1     A       1.7
1     B       1.5
1     B       1.7
2     A       2.1
3     A       1.1
3     A       2.3
3     B       1.1
3     B       2.3

Question 1: I know, if I used: "random = list(~1 | study, ~1| outcome)",
then 'study' and 'outcome' would be crossed, random effects. But are
'study' and 'outcome' in model 'mm' still crossed,
random effects? (In other words, can crossed random-effects have varying
coefficients in them)

Question 2: Is it possible (or useful) to plot residuals vs. fitted values
for model 'mm'? If yes, how can we do that in metafor?

Many thanks for your help,

	[[alternative HTML version deleted]]

More information about the R-sig-meta-analysis mailing list