分类数据统计分析与SAS编程.ppt

一 . 概 述 ;;二 . 2×2 表 ; 2.2.1 Mantel-Haenszel 卡方——Q;如上例江剪刀草组无效（263例）的期望值为： 方差 ： 本例 ： 样本足够大时 近似正太分布，由此得： 结论为两组疗效不同 。;2.2.2 Pearson 卡方——Qp;计算Pearson卡方先要计算出4个期望值 ，一个 已算出为239.87。其他则可由合计值减出。 Mantel-Haenszel 与Pearson 的关系如下： = 本例 ： = ;2.2.3 SAS 程序 [ 程序2-1 ] data bronchit; input treat $ outcome $ count; cards; jiang y 70 jiang n 263 dan y 102 dan n 180 ; proc freq order = data; weight count; tables treat * outcome/chisq; run; ; 检验结果 ： STATISTICS FOR TABLE OF TREAT BY OUTCOME 表中各列依次为统计量、自由度（DF）、卡方值及P值。 ; 2.7 配对资料四格表 ;2.7.3 程序 在SAS第6.10版以后，在proc freq 中tables语句中加入了agree选项作McNemar检验。本例程序如下： [ 程序2-2 ] data culture; input a $ b $ count; cards; y y 36 y n 34 n y 0 n n 135 ； proc freq order = data; weight count; tables a * b/agree; run; ;2.7.4 输出 除频数部分外，输出的结果如下： STATISTICS FOR TABLE OF A BY B McNemar’s Test 其中，Kappa系数是反映一致性的指标。 ;三. 多层 2×2 表;3.2 Mantel-Haenszel 检验 ; 3.2.2 符号 对于多层四格表，除四格表中数字以 代表外，以h表示层次，则第h层的四格表可用以下符号表示（表3-2）: 表3-2 多层四格表的符号表示 ;3.2.3 QMH 统计量 在无效假设为两组疗效相同时，第h层 的期望值 及方差 的计算公式为： = (3.1) = (3.2) 校正中心因素的两组之间疗效的差异可用Mantel-Haenszel 1959年提出的统计量表示。 = 式中q为层数，本例q＝3。 当q层的i行合计数 的合计数 大于30时近似服从自由度等于1的卡方分布。 ;Mantel-Haenszel方法，消除了层次因素的干扰而提高了检出关联的把握度。 当各层次的阳性率的方向一致时，如每层的治愈率都是试验组高于对照组，则QMH检验效果很好。如果各层的方向不一致则可能检验不出其关联。 ;3.2.4 程序 [ 程序3-1] data park2; input center treat $ response $ count @@; cards; 1 test y 12 1 test n 18 1 placebo y 15 1 placebo n 15 2 test y 31 2 test n 9 3 placebo y 34 2 placebo n 6