Skip to content
Prev 2690 / 5632 Next

[R-meta] rma.mv in metafor - model assumptions

Dear Jessie,

rma.mv() does not assume homogeneity of variance (more commonly referred to as homoscedasticity). One supplies the function with the (approximately) known sampling variances via argument 'V' (or possibly even an entire var-cov matrix) and those sampling variances typically differ across studies.

The types of models that one can fit with rma.mv() do assume normality of the sampling errors and of the random effects in the model. That does imply that the observed residuals should also be approximately normally distributed, but they are not expected to have the same standard errors. One can compute standardized residuals that should behave (approximately) like z-scores.

To some degree, one can check normality assumptions by examining the distributions of the estimates of the random effects (from ranef()) and of the residuals (from rstandard()).

What to do if normality assumptions are violated is difficult to say in general. That's a very broad question.

Best,
Wolfgang