Skip to content
Prev 656 / 5632 Next

[R-meta] Meta-analysis of ANOVA interaction effects

Dear Wolfgang,
 
It was a more difficult issue than I thought!
Thank you for your clarification!
 
Dear Iian,
I am sorry for the confusion!
 
Kyungnam
 
-----Original Message-----
From: "Viechtbauer Wolfgang (SP)"<wolfgang.viechtbauer at maastrichtuniversity.nl>
To: "Iain Hamlin"<iain.hamlin at strath.ac.uk>; "r-sig-meta-analysis at r-project.org"<r-sig-meta-analysis at r-project.org>;
Cc:
Sent: 2018-02-27 (?) 23:00:49
Subject: Re: [R-meta] Meta-analysis of ANOVA interaction effects
 
Just compute the 'interaction effects' and corresponding sampling variances and feed this to your meta-analysis package of choice. Let's say there are 2 studies:

dat <- data.frame(study=c(1,2), ti=c(3.41,2.68), n11i=c(32,12),
                 n12i=c(39,13), n21i=c(35,15), n22i=c(36,14))

ti is the t-statistic for the interaction and the n..i variables indicate the cell sizes of the 2x2 factorial design. Then you can compute the standardized effect with:

dat$yi <- with(dat, ti * sqrt(1/n11i + 1/n12i + 1/n21i + 1/n22i))

The corresponding sampling variances can be obtained with:

dat$vi <- with(dat, 1/n11i + 1/n12i + 1/n21i + 1/n22i + dat$yi^2 / (2*(n11i + n12i + n21i + n22i)))

And finally, using for example the metafor package:

library(metafor)
rma(yi, vi, data=dat)

(leaving aside whether fitting a RE model with k=2 is a wise thing to do).

Best,
Wolfgang
_______________________________________________
R-sig-meta-analysis mailing list
R-sig-meta-analysis at r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-meta-analysis