[R-meta] Reasonable rma.mv() syntax for a longitudinal meta-analysis
Simon Harmel
@|m@h@rme| @end|ng |rom gm@||@com
Wed May 5 02:17:21 CEST 2021
Hello List Members,
I'm doing a meta-analysis where the effect sizes (yi) are computed as the
difference in pre-post changes of a treatment group and that in a control
group. That is:
yi = (standardized mean change of a treatment) - (standardized mean change
of a control)
As a result, I have a categorical variable (post_id) that indicates the
pre-post interval associated with each "yi" (`1` = pre to post 1; `2` =
pre to post 2, . . .).
Almost all studies have more than 2 treatment groups, a couple have more
than one control group, and about 15 studies have used more than 1 outcome
but these related outcomes within studies don't represent the same thing
across these 15 studies. In each study, a "yi" for each unique situation
(e.g., having 2 control groups, 2 or more outcomes etc.) has been computed.
I was wondering given my data structure, is my syntax reasonable?
Many thanks, Simon
#==============================
library(metafor)
library(clubSandwich)
dat <- read.csv("https://raw.githubusercontent.com/ilzl/i/master/pr.csv")
dat <- transform(dat, post_id = factor(post_id))
V <- impute_covariance_matrix(vi = dat$vi, cluster = dat$studyID, r = .7)
mv1 <- rma.mv(yi ~ 0 + post_id + x1 + x2, V = V, struct= "AR", random =
list(~post_id|studyID, ~1|esID), data = dat)
[[alternative HTML version deleted]]
More information about the R-sig-meta-analysis
mailing list