公式与数值方法:χ² = Σ(O−E)²/E(E = 行合计×列合计/N,df=(R−1)(C−1));Yates 修正仅用于 2×2:Σ(|O−E|−0.5)²/E;G = 2ΣO·ln(O/E)(O=0 项计 0)。p 值由正则化不完全伽马函数 Q(k/2, χ²/2) 计算(级数展开与连分式双路径,精度 ~1e-7;df=1 时与 erfc(√(χ²/2)) 恒等、df=2 时与 e^(−χ²/2) 恒等,均作自检)。Fisher 2×2 精确检验用 lgamma(Lanczos g=7)计算超几何概率,双尾采用 SciPy 约定(概率 ≤ 观测表的所有表求和);参考例:品茶实验 [[3,1],[1,3]] 双尾 p = 34/70 = 0.4857。效应量:V = √(χ²/(N·min(R−1,C−1))),phi = √(χ²/N)(2×2),C = √(χ²/(χ²+N));V 的经验解读(0.1/0.3/0.5)来自 Cohen。调整残差 d = (O−E)/√(E(1−p行)(1−p列)),|d|>1.96(α=.05)或 >2.58(α=.01)的单元格标注显著并附 BH-FDR 多重比较提示。Cochran 规则:所有 E ≥ 5(或 E<5 的单元格 ≤20% 且无 E<1),否则建议 Fisher/模拟。Chapman 估计(可选,仅 2×2):N̂ = (a+1)(b+1)/(c+1) − 1(Lincoln-Petersen 修正)。检验向量:[[30,20],[20,30]] → χ²=4、p=0.0455、Yates χ²=3.24、G=4.0271、V=0.2、调整残差全部 ±2。