首先把四个表的数据转化为一个四列的矩阵数据，从第一列开始到第四列分别表示四格表中的abcd

for i=1:7
m(i)=a(i,1)/(a(i,1)+a(i,3));
n(i)=a(i,2)/(a(i,2)+a(i,4));
b(i)=[m(i)/(1-m(i))]/[n(i)/(1-n(i))];（or值）
s(i)=exp(1)^[log(b(i))+1.96*(1/a(i,1)+1/a(i,2)+1/a(i,3)+1/a(i,4))^0.5];（上限）
x(i)=exp(1)^[log(b(i))-1.96*(1/a(i,1)+1/a(i,2)+1/a(i,3)+1/a(i,4))^0.5];（下限）
c(i)=[a(i,1)+a(i,2)+a(i,3)+a(i,4)]*[a(i,1)*a(i,4)-a(i,2)*a(i,3)]^2/[(a(i,1)+a(i,2))*(a(i,1)+a(i,3))*(a(i,3)+a(i,4))*(a(i,2)+a(i,4))];（关联卡方值）
d(i)=chi2cdf(c(i),1);
e(i)=1-chi2cdf(c(i),1);（p值）
end

可以用以下例子检验：

六、有300例肝癌病人，某医师对其中105例用甲疗法，其余178例用乙疗法。观察一年后，甲疗法组死亡15例，乙疗法组死亡10例。试分析两组病人的死亡率是否有差异。（列出四格表）

疗法                死亡人数           存活人数             合计

甲疗法

乙疗法

合计
15                 90

10                 168
105

178

25                 258                   283

1.         建立检验假设并确定检验水准

H0两组病人死亡率相同π1=π2

H1两组病人死亡率不同π1≠π2

α=0.05

2.         计算检验统计量

T12=(25×105)／283=9.28

即所有理论频数T＞5,且n＝283＞40

故用专用公式:

X2=(ad－bc)2n／(a+c)(b+d)(a+b)(c+d)＝(15×168－90×10)×283／25×258×105×178＝6.16