### load los data
data(dat.los)
### calculate mean differences and corresponding sampling variances
dat <- escalc(m1i=m1i, sd1i=sd1i, n1i=n1i, m2i=m2i, sd2i=sd2i, n2i=n2i,
measure="MD", data=dat.los, append=TRUE)
dat
### meta-analysis of mean differences using a random-effects model
res <- rma(yi, vi, data=dat)
res
### meta-analysis of standardized mean differences using a random-effects model
res <- rma(m1i=m1i, sd1i=sd1i, n1i=n1i, m2i=m2i, sd2i=sd2i, n2i=n2i,
measure="SMD", data=dat.los, slab=source)
res
### draw forest plot
forest(res, xlim=c(-7,5), alim=c(-3,1))
text(-7, 11, "Study/Source", pos=4)
text( 5, 11, "Observed SMD [95% CI]", pos=2)
Run the code above in your browser using DataLab