Message-ID: <CACgv6yX_5ybyxVU5JyayDegwy=YLt6rL0nz61OnManfGz5w-4A@mail.gmail.com>
Date: 2022-01-13T17:37:30Z
From: Simon Harmel
Subject: [R-meta] Substitute metafor fit for lm fit
In-Reply-To: <CACgv6yVKmvvi-YQ7dNWX3JN1mFdZXmJqm3m_T1UWKJ-DNbTwhw@mail.gmail.com>
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!