需要注意的是卡方检验要求列联表中每格的数值或者期望值大于5,如果该条件不满足,那么R中就会使用Yate’s矫正公式进行计算:
A <- c(rep("male",15),rep("female",20)...,rep("male",15))#创建变量A
B <- c(rep("healthy",4),rep("sick",35),rep("healthy",11))# 创建变量B
C <- c(rep("smoker...",26), rep("nonsmoker",24)) # 创建变量C
mydata <- data.frame(A,B,C) # 利用以创建的变量构建数据框
mytable <- table(A,C)...mytable <- table(C, B, A) #将性别作为调整的因素
?
mantelhaen.test(mytable) #进行CMH检验,判断在矫正性别后,吸烟和健康之间是否有关
?