Subgroup Analysis & Meta-Regression
Slides
Keywords: Fixed-Effects (Plural) Model, Ecological Fallacy
Code
library(meta)
# Copy and paste the entire code below into your Console.
# Then hit Enter.
thirdwave <-
structure(
list(
study = c("Call", "Cavanagh", "Danitz-Orsillo",
"de Vibe", "Frazier", "Frogeli",
"Gallego", "Hazlett-Stevens", "Hintz",
"Kang", "Kuhlmann", "Lever Taylor",
"Phang", "Rasanen", "Ratanasiripong",
"Shapiro", "Song", "Warnecke"),
es = c(0.70, 0.35, 1.79, 0.18, 0.42, 0.63,
0.72, 0.52, 0.28, 1.27, 0.10, 0.38,
0.54, 0.42, 0.51, 1.47, 0.61, 0.60),
se = c(0.260, 0.196, 0.345, 0.117, 0.144,
0.196, 0.224, 0.210, 0.168, 0.337,
0.194, 0.230, 0.244, 0.257, 0.351,
0.315, 0.226, 0.249),
rob = c("high", "low", "high", "low", "low",
"low", "high", "low", "low", "low",
"high", "low", "low", "low", "high",
"high", "high", "low"),
fu.weeks = c(4, 10, 2, 8, 5, 6, 6, 7,
12, 1, 3, 4, 5, 6, 3, 2, 5, 2)),
row.names = c(NA, -18L), class = "data.frame")
# Run REM meta-analysis
rem <- metagen(es, se, study, thirdwave)
# Subgroup analysis
update.meta(rem, by = rob)
# Run meta-regression on follow-up in weeks
metareg(rem, ~ fu.weeks)