data(amlodipine)
m <- metacont(n.amlo, mean.amlo, sqrt(var.amlo),
n.plac, mean.plac, sqrt(var.plac),
data = amlodipine, studlab = study,
method.tau = "DL")
m.hk <- update(m, method.random.ci = "HK")
# Same results for mean difference as in Table III in Hartung and
# Knapp (2001)
#
vars.common <- c("TE.common", "lower.common", "upper.common")
vars.random <- c("TE.random", "lower.random", "upper.random")
#
res.common <- as.data.frame(m[vars.common])
names(res.common) <- vars.random
#
res.md <- rbind(res.common,
as.data.frame(m[vars.random]),
as.data.frame(m.hk[vars.random]))
#
res.md <- round(res.md, 5)
#
row.names(res.md) <- c("CE", "RE", "RE (HaKn)")
names(res.md) <- c("Absolute difference", "CI lower", "CI upper")
#
res.md
Run the code above in your browser using DataLab