# Mentel-Haneszel Method with Odds Ratio Estimation
library("lawstat")

boymatrix <- matrix(c(20,100,15,150),nrow=2,byrow=TRUE,   # enter data counts
            dimnames = list("Sleep"  = c("Low","High"),   # label column variable categories
                            "Result" = c("Fail","Pass"))) # label row variable categories
print(boymatrix)


girlmatrix <- matrix(c(30,100,25,200),nrow=2,byrow=TRUE,  # enter data counts
            dimnames = list("Sleep"  = c("Low","High"),   # label column variable categories
                            "Result" = c("Fail","Pass"))) # label row variable categories
print(girlmatrix)


myarray <- array(c(boymatrix,girlmatrix),dim=c(2,2,2))    # prepare data for M-H method
cmh.test(myarray)                                         # run Mantel-Haneszel Method

library("DescTools")

BreslowDayTest(myarray, OR = NA, correct = FALSE)