适合性检验的SAS检验方法

例5.1　统计某羊场一年所产876 只羔羊，有公羔428 只，母羔448 只。试检验该羊场所产羔羊的性别是否符合1∶1的理论比例？

这是一个适合性检验的简单例子，按照1∶1的理论比例，该群羔羊的公羔数等于母羔数，即都为438 只。SAS 程序如图5.2 所示。

输出结果为：输出结果中分别给出了χ2值（Chi-square）、矫正的χ2值（Continuity adj，chisquare）、似然比χ2值（Likelihood ratio chi-square）、M-H 法χ2值（Mantel-Haenszel chi-square）、Fisher 精确概率值、Phi 系数和列联系数（Contingency Coefficient）等。在实际分析中，应根据具体资料选择不同的结果。

一般地，当自由度df>1，或理论次数E≥5 且样本数n≥30 时，选择χ2值；当自由度df=1，或1≤E<5，n≥30 时，选用矫正的χ2值；当E<1 或n<30 时，选用Fisher 精确概率结果；当对大样本条件的近似分布进行检验时，选用似然比χ2值；而M-H 法χ2值主要是指卡方检验中的分层分析，在畜牧试验分析中极少使用。

图5.2　例5.1 SAS 程序图（1）

本例χ2检验的结果，，表明该群羔羊的性别比例符合1∶1的理论比例。

其SAS 程序如图5.3、图5.4 所示。

(www.xing528.com)

图5.3　例5.1 SAS 程序图（2）

图5.4　例5.1 SAS 程序图（3）

例5.2　两对相对性状在F2的4 种表现型A_B_、A_bb、aaB_、aabb 的实际观察次数依次为152、39、53、6。检验这两对相对性状在F2的4 种表现型A_B_、A_bb、aaB_、aabb 的分离是否符合9∶3∶3∶1的理论比例？

这是一个4 个属性类别的适合性检验问题，k=4，df=k -1=3 >1，进行χ2检验不需要作连续性矫正，SAS 程序如图5.5 所示。

图5.5　例5.2 SAS 程序图

输出结果为：本例χ2检验的结果，χ2=5.824，P=0.120 >0.05，表明4 种基因型的观察次数符合9∶3∶3∶1的理论比例。