Skip to content
Back to formatted view

Raw Message

Message-ID: <1c160ae463574e188f39409a9cd4587a@UM-MAIL3214.unimaas.nl>
Date: 2019-09-11T07:17:03Z
From: Wolfgang Viechtbauer
Subject: [R-meta] Residuals for random effect

Hi Lukas,

I assume you are after the BLUPs of the random effects. You can obtain those with the ranef() function. So, just do ranef(res1).

Best,
Wolfgang

-----Original Message-----
From: R-sig-meta-analysis [mailto:r-sig-meta-analysis-bounces at r-project.org] On Behalf Of Lukas Dylewski
Sent: Tuesday, 10 September, 2019 14:30
To: r-sig-meta-analysis at r-project.org
Subject: [R-meta] Residuals for random effect

Dear All,

I would like to calculate residuals for a random effect.  Could I calculate
the residual values after subtracting the predicted influence of the
study?
I would like to get the same residuals value for each study level (in the
model below coded pub).

My model:
res1 <- rma.mv(yi, vi, mods = ~ logmass +
I(logmass^2),random=~1|pub,data=metaan, method="REML")

Regards
Lukasz
-- 
?ukasz Dylewski, M. Sc.
PhD Student

Department of Zoology, Institute of Zoology
Pozna? University of Life Sciences
Wojska Polskiego 71C
60-625 Pozna?, Poland

http://ecology1.wixsite.com/farmlandecology