# NOT RUN {
data(Fleiss1993cont)
#
# Generate additional variable with grouping information
#
Fleiss1993cont$group <- c(1, 2, 1, 1, 2)
#
# Do meta-analysis without grouping information
#
m1 <- metacont(n.psyc, mean.psyc, sd.psyc, n.cont, mean.cont, sd.cont,
data = Fleiss1993cont, sm = "SMD",
studlab = paste(study, year))
#
# Update meta-analysis object and do subgroup analyses
#
summary(update(m1, byvar = group))
# Same result using metacont function directly
#
m2 <- metacont(n.psyc, mean.psyc, sd.psyc, n.cont, mean.cont, sd.cont,
data = Fleiss1993cont, sm = "SMD",
studlab = paste(study, year), byvar = group)
summary(m2)
# Compare printout of the following two commands
#
as.data.frame(m1)
m1$data
# }
Run the code above in your browser using DataLab