非参数统计的Python实现—— Mantel-Haenszel 检验

概念

很多研究都涉及分层数据结构，比如产品研究中需要根据城市和农村特点分别研究不同人群对产品的满意程度，这里城市和农村是问题的两个层，层中不同人群对产品的满意程度可以用列联表表示。回答处理与反应结果是否独立的问题时，需要首先按层计算差异，再将各层的差异进行综合比较，从而做出综合的判断。在本文中笔者仅考虑各层 2*2 的列联表处理，且试验组的阳性结果放在列联表第一行第一列位置。

实例 & 代码

对两家医院考察某治癌药的治癌效果，试验组 A 与对照组 B (安慰剂) 对比记录其疗效如下表所示。问治癌药是否有显著效果。
解答：

import pandas as pd
import numpy as np
import scipy.stats as stats
arr1=np.array([[50,15],[92,90]])
arr2