卡方上机实习程序.docVIP

两组率的比较 实习指导 P46 例8.1 data ex8_1; /*定义数据集名为ex8_1*/ input r c number; /*定义输入的变量为r c number*/ cards; /*直接输入数据，数据块开始*/ 1 1 54 1 2 8 2 1 44 2 2 20 ; run; proc freq data=ex8_1; /*调用freq过程，指定数据集*/ tables r*c/chisq expected;/*输出r*c表，注意顺序，行变量*列变量；同时利用Chisq选项作卡方检验，要求输出理论频数*/ weight number; /*指定表示频数的变量*/ run; The SAS System The FREQ Procedure Frequency Expected Percent Row Pct Col Pct Table of r by c r c Total 1 2 1 5448.22242.8687.1055.10 813.7786.3512.9028.57 62?49.21?? 2 4449.77834.9268.7544.90 2014.22215.8731.2571.43 64?50.79?? Total 9877.78 2822.22 126100.00 Statistics for Table of r by c Statistic DF Value Prob Chi-Square 1 6.1331 0.0133 Likelihood Ratio Chi-Square 1 6.3036 0.0120 Continuity Adj. Chi-Square 1 5.1175 0.0237 Mantel-Haenszel Chi-Square 1 6.0844 0.0136 Phi Coefficient ? 0.2206 ? Contingency Coefficient ? 0.2154 ? Cramers V ? 0.2206 ? ? Fishers Exact Test Cell (1,1) Frequency (F) 54 Left-sided Pr = F 0.9968 Right-sided Pr = F 0.0112 ? ? Table Probability (P) 0.0080 Two-sided Pr = P 0.0179 Sample Size = 126 配对四格表的卡方检验 实习指导 P50 例8.3 data ex8_3; /*定义数据集名为ex8_3*/ input r c number; /*定义输入的变量为r c number*/ cards; /*直接输入数据，数据块开始*/ 1 1 12 1 2 16 2 1 4 2 2 13 ; proc freq; /*调用freq过程，指定数据集*/ tables r*c/agree; /*注意配对四格表要用选项/agree，SAS中不对b+c40的情况做校正，应用时请自己做校正*/ weight number; /*指定表示频数的变量*/ run; The SAS System The FREQ Procedure Frequency Percent Row Pct Col Pct Table of r by c r c Total 1 2 1 1226.6742.8675.00 1635.5657.1455.17 2862.22?? 2 48.8923.5325.00 1328.8976.4744.83 1737.78?? Total 1635.56 2964.44 45100.00 Statistics for Table of r by c McNemars Test Statistic (S) 7.2000 DF 1 Pr S 0.0073 ? Simple Kappa Coefficient Kappa 0.1697 ASE 0.1244 95% Lower Conf Limit -0.0740 95% Upper Conf Limit 0.4135 Sample Size = 45 行列表资料的卡方检验 实习指导 P49 例8.2 data ex8_2; /*定义数据集名为ex8_2*/ do a=1 to 4; /*循环开始，循环控制变量为a,取值从1到4*/ do b=1 to 3; /*循