mantelCoxHR: Mantel-Cox estimator of the hazard ratio

In felix-hof/biostatUZH: Misc Tools of the Department of Biostatistics, EBPI, University of Zurich

Mantel-Cox estimator of the hazard ratio

Description

Estimate the hazard ratio and compute a confidence interval for it via a special application of the Mantel-Haenszel procedure. A separate 2 x 2 table is constructed for each event time. The underlying assumption is that the hazard ratio is constant over the follow-up period.

Usage

mantelCoxHR(time, event, group, conf.level = 0.95)

Arguments

time	Event times, censored or observed.
event	Censoring indicator, 1 for event, 0 for censored.
group	Factor with two levels, e.g. treatment group.
conf.level	Significance level for confidence interval for hazard ratio.

Value

mantelCox.hr	Hazard ratio estimate.
ci.hr	Wald confidence interval at the level specified by conf.level.
p.val.logrank	p-value of logrank test for a comparison of the survival curves between the two groups.
coxph.hr	Hazard ratio estimated via Cox-regression.

Note

In general the Mantel-Cox estimate and the hazard ratio estimate received via Cox-regression do not coincide.

Author(s)

Kaspar Rufibach
[e-mail]

References

Kirkwood, B.R. and Sterne, J.A.C. (2003). Essential Medical Statistics. Blackwell Science. See p. 283 ff.

Examples

## load Acute Myelogenous Leukemia survival data
data("leukemia", package = "survival")

mantelCoxHR(time = leukemia[, 1],
            event = leukemia[, 2],
            group = as.factor(leukemia[, 3]))

felix-hof/biostatUZH documentation built on Sept. 27, 2024, 1:48 p.m.