# 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)