mantelHaenszel {biostatUZH} | R Documentation |
Compute Cochran-Mantel-Haenszel chi-squared test of the null that two nominal variables are conditionally independent in each stratum, assuming that there is no three-way interaction.
mantelHaenszel(exposure, outcome, stratum)
exposure |
Binary variable coding whether patient was exposed (1) or not (0). |
outcome |
Binary variable coding outcome of patient. |
stratum |
Nominal variable containing information about matching, i.e. strata. |
tab |
Table that counts numbers of strata for each case-control combination. |
test.stat |
Test statistic for chi^2 test. |
p.val |
p-value of chi^2 test. |
Kaspar Rufibach
kaspar.rufibach@gmail.com
Agresti, A. (2002). Categorical data analysis. New York: Wiley.
Similar functionality is provided in mantelhaen.test
and clogit
in survival.
See the examples below for a comparison.
# generate data set.seed(1977) exposure <- rep(c(1, 0, 0, 0, 0), 41) outcome <- sample(c(rep(1, 62), rep(0, 5 * 41 - 62))) strata <- rep(1:41, each = 5) # via conditional logistic regression logreg <- clogit(outcome ~ exposure + strata(strata), method = "approximate") summary(logreg) # R function in library 'stats' mh <- mantelhaen.test(x = outcome, y = exposure, z = strata) # this function mH <- mantelHaenszel(exposure, outcome, strata) # compare p-values summary(logreg)$coef[5] mh$p.value mH$p.val