第十 秩和检验.pptVIP

第十章 秩 和 检 验 （一）方法步骤 data li10_1; input x1 x2@@; d=x1-x2; cards; 39 55 42 47 51 53 43 41 55 54 45 63 22 42 48 46 40 43 45 49 40 37 49 52 ; proc univariate normal; var d ; run; 0.01 0.02 0.05 0.10 0.10 0.05 0.02 0.01 76 79 84 89 141 146 151 154 data li10_3; do c=1 to 2; input n@@; do i=1 to n; input x@@; output; end; end; cards; 10 10 12 15 15 16 17 18 20 23 90 12 2 3 4 5 6 7 8 9 10 11 12 13 ; proc univariate normal; class c; var x; proc npar1way wilcoxon; class c; var x; run; 4、确定P值和得出推断结论 查附表2 t界值表，?＝?得单侧P0.01，按?＝0.05水准拒绝H0，接受H1，故可以认为家庭功能良好的糖尿病患者的血糖控制情况优于家庭功能障碍的糖尿病患者。 Data li0_4; input c x f@@; cards; 1 1 6 1 2 12 1 3 20 2 1 30 2 2 28 2 3 4 ; proc univariate normal; class c; var x; freq f; proc npar1way; class c; var x; freq f; run; 第三节 完全随机设计多样本比较的秩和检验 一、多组数值变量资料的秩和检验 例10-5 为研究精氨酸对小鼠截肢后淋巴细胞转化功能的影响，将21只小鼠分成3组：A组为对照组，B组为截肢组，C组为截肢加精氨酸治疗组。观察脾淋巴细胞对肝素酶（HPA）刺激的增值反应，测量指标是3H吸收量（cpm），数据如表10-5所示，试分析各组测量值是否不同。（经检验这三组来自的总体的方差不齐） 1、建立假设 H0:三组3H吸收量的总体分布相同 H1:三组3H吸收量的总体分布不同或不完全 相同 ?＝0.05 2、编秩 先将两组数据放在一起，从小到大统一编秩。编秩方法同两样本秩和检验。见表10-5（2）、（4）、（6）栏。 3、求秩和 分别将各组秩次相加，分别求得R1、 R2和 R3。 4、计算统计量 （10－5） 式中Ri为各组的秩和，ni为各组对应的例数，N：总例数，本例N＝21。 5、确定P值和作出推断结论 （1）当组数=3，每组例数ni?5，可查附表12H界值表得到P值。 （2）当不满足条件（1）时，H近似地服从自由度k-1的?2分布，可查?2界值表得P值。 本例ni＝75,故查?2界值表，?20.05，2＝5.99，H＝9.848 ?20.05，2 ＝5.99,所以，P0.05,按?＝0.05水准拒绝H0，接受H1，故可以认为三组3H吸收量不同或不完全相同。 当相同的数值变量较多时，应按公式（10－6）计算校正值Hc data li10_5; do c=1 to 3;do i=1 to 7; input x@@;output;end;end; cards; 3012 9458 8419 9580 13590 12787 6600 2532 4682 2025 2268 2775 2884 1717 8138 2073 1867 885 6490 9003 690 ; proc univariate normal;class c;var x; proc anova;class c;model x=c; means c/hovtest; proc npar1way wilcoxcon;class c;var x; proc freq ; tables c*x/scores