[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 ?

fitstats(res)

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 :)