# \donttest{
data(lim)
lim$vi <- (1/sqrt(lim$N - 3))^2
lim_MR <- metafor::rma.mv(
yi = yi, V = vi, mods = ~ Phylum - 1,
random = list(~1 | Article, ~1 | Datapoint),
data = lim)
R2 <- r2_ml(lim_MR, data = lim)
# }
Run the code above in your browser using DataLab