[R-meta] Pairwise moderator testing in multilevel meta-analysis with CRVE / CIs
Röhl, Sebastian
@eb@@t|@n@roeh| @end|ng |rom un|-tueb|ngen@de
Thu Feb 9 12:29:37 CET 2023
Hi,
I have the following problem:
I am conducting a multilevel meta-analysis using metafor with cluster robust variance estimation and want to test the moderating effect of different kinds of outcomes. Additionally I want to test whether the several outcomes differ significantly from each other.
Here is an example:
out_3 <- rma.mv(zr, V=var, random = ~ 1| Sample_ID / number, mods = ~ -1 + out_intg + out_acad + out_socem,
data = data_int)
out_3_rob <- robust(out_3, Sample_ID, clubSandwich = T)
anova(out_3_rob, X=rbind(c(-1,1,0),c(-1,0,1), c(0,-1,1)))
The robust model result shows C.I. that overlap.
Model Results:
estimate se¹ tval¹ df¹ pval¹ ci.lb¹ ci.ub¹
out_intg 0.2302 0.0231 9.9484 30.84 <.0001 0.1830 0.2773 ***
out_acad 0.1646 0.0220 7.4677 17.36 <.0001 0.1182 0.2111 ***
out_socem 0.2458 0.0278 8.8510 22.27 <.0001 0.1882 0.3034 ***
BUT the anova results show significant differences between 2 outcomes:
Hypotheses:
1: -out_intg + out_acad = 0
2: -out_intg + out_socem = 0
3: -out_acad + out_socem = 0
Results:
estimate se tval df pval
1: -0.0655 0.0322 -2.0349 10.92 0.0669
2: 0.0157 0.0330 0.4742 13.21 0.6431
3: 0.0812 0.0188 4.3264 9.60 0.0016
Do I have a thinking error here about the ANOVA or is this pairwise testing not possible with the CRVE-results?
Perhaps I am also interpreting the C.I.s incorrectly? If I calculate a pairwise comparison with the non-robust model, I also get significant difference although also the non-robust C.I. overlap.
Thank you very much for your help!
Best,
Sebastian
****************************
Dr. Sebastian Röhl
Eberhard Karls Universität Tübingen
Institute for Educational Science
Tübingen School of Education (TüSE)
Wilhelmstraße 31 / Room 302
D-72074 Tübingen
Germany
Phone: +49 7071 29-75527
Fax: +49 7071 29-35309
Email: sebastian.roehl using uni-tuebingen.de<mailto:sebastian.roehl using uni-tuebingen.de>
Twitter: @sebastian_roehl @ResTeacherEdu
[[alternative HTML version deleted]]
More information about the R-sig-meta-analysis
mailing list