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

1
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.

Cox-reg

Hazard ratio estimated via Cox-regression.

Note

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

Author(s)

Kaspar Rufibach
kaspar.rufibach@gmail.com

References

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

Examples

1
2
3
## use Acute Myelogenous Leukemia survival data contained in package 'survival'
time <- leukemia[, 1]; status <- leukemia[, 2]; x <- as.factor(leukemia[, 3])
mantelCoxHR(time, status, x)

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.