第二章 SAS 第三节卡方检验 - 2013年《SAS软件进行统计分析应用时的解析》课件.pptVIP

卡平方?2检验 一、 适合性检验 【例2-4】调查某水稻品种F2群体的米粒特性，有红色非糯米粒500粒，红色糯性米粒172粒，白色非糯米粒190粒，白色糯性米粒55粒，试测验这一资料中红色非糯、红色糯性、白色非糯、白色糯性四种米粒数量的实际观察比例是否符合于9:3:3:1的理论比值。 SAS编程 可以使用PROC FREQ过程进行适合性测验。适合性测验的SAS程序如下： data aa; input tr $ number @@; datalines; A 500 B 172 C 190 D 55 ; run; proc freq data=aa; weight number; tables tr / testF=(515.8125 171.9375 171.9375 57.3125); run; 在上述程序中，TESTF指令测验的假设分布为变量各水平的理论次数，括号中即为变量各水平的理论次数。也可以对变量各水平的理论比率进行测验，其指令为testP，括号中即为变量各水平的理论比率，理论比率总和为1。对于上程序，testF=(515.8125 171.9375 171.9375 57.3125)可替换为testP=（0.5625 0.1875 0.1875 0.0625）。 程序运行结果如图3-1。卡方值为1.5000，概率为0.4724，未达到显著，故接受H0，即红色非糯、红色糯性、白色非糯、白色糯性四种米粒数量的实际观察比例是否符合于9:3:3:1的理论比值。 二、 独立性检验 【例2-2】为研究追施氮肥与水稻倒伏性的关系，进行水稻灌浆期追施氮肥的试验，获得不同施肥处理下的倒伏株数的数据，列于表3-2，试分析水稻植株倒伏与追施氮肥是否有关。 1、SAS编程 SAS程序如下： data tem; do tr=1 to 2; do df=1 to 2; input x @@; output; end; end; cards; 132 106 49 182 ; proc freq; tables tr*df / chisq; weight x;/* 语句指定存放表格单元观测数的变量 */ run; 其结果包括两向表、统计量和Fisher精确检验，其中统计量结果如图3-3，其中统计量有卡方和连续校正卡方，对于自由度为1时，统计量结果选卡方，对于自由度大于1时，统计量结果选卡方值。对于本例，自由度为1，因此，选连续校正卡方值，即卡方值为56.5885，无效假设H0为正确卡平方值的概率为0.0001，达到极显著，否定H0，推断水稻植株倒伏与追施氮肥有关。 图3-3 独立性测验的SAS程序分析结果 2、SAS/Analyst法 1. 启动软件，点击解决方案 → 分析 → 分析家。 2. 输入数据。或打开已有的数据集，或由其它格式文件导入。若用excel文件导入，事先将数据表整理成“Danfei”（氮肥）、“Daofu”（倒伏）、“Zhushu”（株数）形式 。 3. 点击统计→假设检验→表分析。 4. 弹出“Table Analyst：imported”（表分析：导入）对话框，从左侧候选变量列表中，点选变量“Danfei” → 点击按钮“Row”，进入“Row”（行）变量列表框内，再点选变量“Daofu”→ 点击按钮“Column”（列）变量列表框内，再点选变量“zhushu”→ 点击按钮“Cell Counts”（单元格计数）列表框内，见图3-4； 图3-4 Table Analyst：imported （表分析：导入）对话框 5. 点击“Statistics”（统计）按钮，弹出“Table Analyst：Statistics”（表分析：统计）对话框，点选“Chi-square statistics”（图3-5），其他选项保留默认设置 →“OK”。 6. 点击“Tables”（表）按钮，弹出“Table Analyst：Tables”（表分析：表）对话框，在 图3-6 Table Analyst：Tables（表分析：表）对话框 “Frequencies”（频数）选项栏，点选“Observed”（观察的）；在“Percentage”（百分比）选项栏，点选“Cell”、“Row”、“Column”（图3-6），其他选项保留默认设置 → “OK” → “OK”。结果输出于新生成的Analysis窗口。同编程输出结果。 图3-6 Table Analyst：Tables （表分析：表）对话框 (二) r×c表的独立性检验 【例3-3】多小叶大豆突变系牡5796-3的叶片均有5片复叶率，以其作为父本分别与小叶正常的栽培大豆品种垦丰16、合丰50、垦鉴豆43、绥农28杂