Cisapride in Non-Ulcer Dispepsia
Meta-analysis on cisapride in non-ulcer dispepsia.
This meta-analysis is used as a data example in Hartung and Knapp (2001).
A data frame with the following columns:
| study | study label |
| event.cisa | number of events in cisapride group |
| n.cisa | number of observations in cisapride group |
| event.plac | number of events in placebo group |
| n.plac | number of observations in placebo group |
Hartung J & Knapp G (2001): A refined method for the meta-analysis of controlled clinical trials with binary outcome. Statistics in Medicine, 20, 3875–89
data(cisapride)
m.or <- metabin(event.cisa, n.cisa, event.plac, n.plac,
data = cisapride, sm = "OR", method = "Inverse",
studlab = study, addincr = TRUE)
s1 <- summary(update(m.or, sm = "RR"))
s2 <- summary(update(m.or, sm = "RR", hakn = TRUE))
#
s3 <- summary(m.or)
s4 <- summary(update(m.or, hakn = TRUE))
vars <- c("TE", "lower", "upper")
# Results for log risk ratio - see Table VII in Hartung and Knapp (2001)
#
res.rr <- rbind(data.frame(s1$fixed)[vars],
data.frame(s1$random)[vars],
data.frame(s2$random)[vars])
#
row.names(res.rr) <- c("FE", "RE", "RE (HaKn)")
names(res.rr) <- c("Log risk ratio", "CI lower", "CI upper")
#
res.rr
# Results for log odds ratio (Table VII in Hartung and Knapp 2001)
#
res.or <- rbind(data.frame(s3$fixed)[vars],
data.frame(s3$random)[vars],
data.frame(s4$random)[vars])
#
row.names(res.or) <- c("FE", "RE", "RE (HaKn)")
names(res.or) <- c("Log odds ratio", "CI lower", "CI upper")
#
res.orPlease choose more modern alternatives, such as Google Chrome or Mozilla Firefox.