[R-meta] Metafor modeling question

[Yuhang Hu]

Dear All,

I have two questions about the following model:

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

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,
Yuhang

	[[alternative HTML version deleted]]