data(woodyplants)
# Meta-analysis of response ratios (Hedges et al., 1999)
#
m1 <- metacont(n.elev, mean.elev, sd.elev, n.amb, mean.amb, sd.amb,
data = woodyplants, sm = "ROM", studlab = paste(obsno, papno, sep = " / "))
print(m1, prediction = TRUE)
# Meta-analysis for plants grown with low soil fertility treatment
#
m2 <- update(m1, subset = (treat == "fert" & level == "low"))
print(m2, prediction = TRUE)
# Meta-analysis for plants grown under low light conditions
#
m3 <- update(m1, subset = (treat == "light" & level == "low"))
print(m3, prediction = TRUE)
Run the code above in your browser using DataLab