Hello Wolfgang,
A quick question. Is it possible to substitute the vcov(rma.mv_fit)
into a vcov(lm_fit)?
Thank you,
Simon
# An Example:
library(metafor)
dat2 <- escalc(measure="OR", ai=waward, n1i=wtotal, ci=maward,
n2i=mtotal, data=dat.bornmann2007)
res1 <- rma.mv(yi ~ 0+type, vi, random = ~ 1 | study/obs, data=dat2)
res2 <- lm(yi ~ 0+type, data = dat2)
vcov(res2) <- vcov(res1) ## apparently this won't work!