# NOT RUN {
data(amlodipine)
m <- metacont(n.amlo, mean.amlo, sqrt(var.amlo),
n.plac, mean.plac, sqrt(var.plac),
data = amlodipine, studlab = study)
s1 <- summary(m)
s2 <- summary(update(m, hakn = TRUE))
vars <- c("TE", "lower", "upper")
# Same results for mean difference as in Table III in Hartung and
# Knapp (2001)
#
res.md <- rbind(data.frame(s1$fixed)[vars],
data.frame(s1$random)[vars],
data.frame(s2$random)[vars])
#
res.md <- round(res.md, 5)
#
row.names(res.md) <- c("FE", "RE", "RE (HaKn)")
names(res.md) <- c("Absolute difference", "CI lower", "CI upper")
#
res.md
# }
Run the code above in your browser using DataLab