-----Original Message-----
From: R-sig-meta-analysis [mailto:r-sig-meta-analysis-bounces at r-project.org] On
Behalf Of Martineau, Roger
Sent: Thursday, 24 March, 2022 13:51
To: r-sig-meta-analysis at r-project.org
Subject: Re: [R-meta] extract pval from robust function
Hi Wolfgang,
Thanks for your answer, it works.
One more thing, using the dat.konstantopoulos2011 example, I asked for the
fitstats of res. I need to extract the AICc value, i.e. 23.3; how can I do it ?
Also, I assume that those fitstats data are valid using the robust function
results; are they ?
REML
logLik: -7.24
deviance: 14.48
AIC: 22.48
BIC: 30.44
AICc: 23.30
Thanks in advance,
Roger ?
-----Original Message-----
From: Viechtbauer, Wolfgang (SP) <wolfgang.viechtbauer at maastrichtuniversity.nl>
Sent: Thursday, March 24, 2022 8:17 AM
To: Martineau, Roger <roger.martineau at AGR.GC.CA>; r-sig-meta-analysis at r-
project.org
Subject: RE: [R-meta] extract pval from robust function
Hi Roger,
robust(res, cluster=dat$district)$pval[2]
Use str() to find out about the structure of an object, so:
str(robust(res, cluster=dat$district))
Scanning through all that stuff will often lead to the answer.
Best,
Wolfgang
-----Original Message-----
From: R-sig-meta-analysis
[mailto:r-sig-meta-analysis-bounces at r-project.org] On Behalf Of
Martineau, Roger
Sent: Thursday, 24 March, 2022 12:58
To: r-sig-meta-analysis at r-project.org
Subject: [R-meta] extract pval from robust function
Dear list members,
In this example, the p-value of year in the robust function is 0.5885;
what are the codes to extract that p-value ?
dat <- dat.konstantopoulos2011
res <- rma.mv(yi ~ year, vi, random = ~ 1 | district/school, data=dat)
robust(res, cluster=dat$district)
Thanks in advance.
Roger :)