[R-meta] Substitute metafor fit for lm fit

Simon Harmel @|m@h@rme| @end|ng |rom gm@||@com
Thu Jan 13 18:37:30 CET 2022


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!



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