Test if the proportions of two dichotomous variables are equal in the same population.
mcnemar_test(x, y = NULL)
mcnemar_test
returns an object of class "mcnemar_test"
.
An object of class "mcnemar_test"
is a list containing the
following components:
Sheskin, D. J. 2007. Handbook of Parametric and Nonparametric Statistical Procedures, 4th edition. : Chapman & Hall/CRC.
# test if the proportion of students in himath and hiread group is same himath <- ifelse(hsb$math > 60, 1, 0) hiread <- ifelse(hsb$read > 60, 1, 0) mcnemar_test(table(himath, hiread))#> Controls #> --------------------------------- #> Cases 0 1 Total #> --------------------------------- #> 0 135 21 156 #> 1 18 26 44 #> --------------------------------- #> Total 153 47 200 #> --------------------------------- #> #> McNemar's Test #> ---------------------------- #> McNemar's chi2 0.2308 #> DF 1 #> Pr > chi2 0.631 #> Exact Pr >= chi2 0.7493 #> ---------------------------- #> #> Kappa Coefficient #> -------------------------------- #> Kappa 0.4454 #> ASE 0.075 #> 95% Lower Conf Limit 0.2984 #> 95% Upper Conf Limit 0.5923 #> -------------------------------- #> #> Proportion With Factor #> ---------------------- #> cases 0.78 #> controls 0.765 #> ratio 1.0196 #> odds ratio 1.1667 #> ----------------------# using matrix mcnemar_test(matrix(c(135, 18, 21, 26), nrow = 2))#> Controls #> --------------------------------- #> Cases 0 1 Total #> --------------------------------- #> 0 135 21 156 #> 1 18 26 44 #> --------------------------------- #> Total 153 47 200 #> --------------------------------- #> #> McNemar's Test #> ---------------------------- #> McNemar's chi2 0.2308 #> DF 1 #> Pr > chi2 0.631 #> Exact Pr >= chi2 0.7493 #> ---------------------------- #> #> Kappa Coefficient #> -------------------------------- #> Kappa 0.4454 #> ASE 0.075 #> 95% Lower Conf Limit 0.2984 #> 95% Upper Conf Limit 0.5923 #> -------------------------------- #> #> Proportion With Factor #> ---------------------- #> cases 0.78 #> controls 0.765 #> ratio 1.0196 #> odds ratio 1.1667 #> ----------------------