Learn R Programming

meta (version 4.18-0)

cisapride: Cisapride in Non-Ulcer Dispepsia

Description

Meta-analysis on cisapride in non-ulcer dispepsia.

This meta-analysis is used as a data example in Hartung and Knapp (2001).

Arguments

Format

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

See Also

metabin

Examples

Run this code
# NOT RUN {
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.or
# }

Run the code above in your browser using DataLab