《卫生统计学 实验课件 秩和检验(npar1way过程)》由会员分享，可在线阅读，更多相关《卫生统计学 实验课件 秩和检验(npar1way过程)（30页珍藏版）》请在金锄头文库上搜索。
1、秩和检验（非参数检验）,实习十二,目的要求,掌握利用univariate过程完成配对设计资料的非参数检验； 掌握利用npar1way过程及选择项Wilcoxon的使用完成完全随机设计两组、多组资料的秩和检验； 熟悉利用freq过程完成随机区组设计资料的秩和检验且能对结果作分析。,非参数统计的适用情况,（1）等级资料； （2）偏态分布； （3）分布不明； （4）个别数据偏离过大； （5）各组方差明显不齐。,非参数统计的优缺点,优点：不受总体分布的限制，适用范围广。 缺点：没有充分利用资料提供的信息。检验效能低于参数检验。适用参数检验条件的资料，应首选参数检验。,例1某研究者欲研究保健食品对小鼠抗
2、疲劳作用，将同种属的小鼠按性别和年龄相同、体重相近配成对子，共10对，并将每对中的两只小鼠随机分到保健食品两个不同的剂量组，过一定时期将小鼠杀死，测得其肝糖原含量（mg/100g），结果见下表， 问不同剂量组的小鼠肝糖原含量有无差别？,配对资料的秩和检验,data SAS1; input x1 x2 ; d=x1-x2; Cards; 620.16 958.47 866.50 838.42 641.22 788.90 812.91 815.20 738.96 783.17 899.38 910.92 760.78 758.49 694.95 870.80 749.92 862.26 793.9
3、4 805.48 ;,proc univariate normal; var d; run;,Tests for Location: Mu0=0Test -Statistic- -p Value-Students t t -2.25916 Pr |t| 0.0502Sign M -3 Pr = |M| 0.1094Signed Rank S -21 Pr = |S| 0.0313Tests for NormalityTest - -Statistic- -p Value-Shapiro-Wilk W 0.843269 Pr D 0.1362 Cramer-von Mises W-Sq 0.10
4、6441 Pr W-Sq 0.0824 Anderson-Darling A-Sq 0.657567 Pr A-Sq 0.0625,符号秩和,单样本秩和检验,data a; input nf; cards; 2.15 2.10 2.20 2.12 2.42 2.52 2.62 2.72 2.99 3.19 3.37 4.57 ; proc univariate normal mu0=2.15; var nf; run;,Tests for Location: Mu0=2.15Test -Statistic- -p Value-Students t t 2.898449 Pr |t| 0.014
5、5Sign M 3.5 Pr = |M| 0.0654 Signed Rank S 29.5 Pr = |S| 0.0059Tests for NormalityTest -Statistic- -p Value- Shapiro-Wilk W 0.837966 Pr D 0.1500 Cramer-von Mises W-Sq 0.092826 Pr W-Sq 0.1281 Anderson-Darling A-Sq 0.651796 Pr A-Sq 0.0704,例2 为比较两种刺激物的对人的影响力的持续时间，现随机抽取13个试验对象接受第一种刺激物，6个试验对象接受第二种刺激物。持续时间
6、如下表，问二种刺激物的持续时间是否有差异？,两独立样本秩和检验,持续时间表(minutes),PROC NPAR1WAY过程格式,PROC NPAR1WAY Wilcoxon ；CLASS 变量；VAR 变量表；RUN；,data sas1； do a=1 to 2； input n； do j=1 to n； input x； output； end； end； cards；,13 1.94 1.94 2.92 2.92 2.92 2.92 3.27 3.27 3.27 3.27 3.70 3.70 3.74 6 3.27 3.27 3.27 3.70 3.70 3.74 ；,proc so
7、rt；排序过程 by a； run； proc univariate normal；考虑各组的正态性 var x；by a； run； proc npar1way wilcoxon； class a； var x；exact；输出精确的概率 run；,结果,The NPAR1WAY ProcedureWilcoxon Scores (Rank Sums) for Variable xClassified by Variable aSum of Expected Std Dev Mean a N Scores Under H0 Under H0 Score - 1 13 110.50 130.0
8、 11.004784 8.500 2 6 79.50 60.0 11.004784 13.250,每组得分和 H0成立时的得分和 H0成立时的标准差 平均得分,Average scores were used for ties.Wilcoxon Two-Sample TestStatistic (S) 79.5000Normal ApproximationZ 1.7265One-Sided Pr Z 0.0421Two-Sided Pr |Z| 0.0843t ApproximationOne-Sided Pr Z 0.0507Two-Sided Pr |Z| 0.1014,T统计量,正态近似
9、法 结果（超出 表范围时使 用）,近似t检验结果,Exact TestOne-Sided Pr = S 0.0527Two-Sided Pr = |S - Mean| 0.1054,小样本时可通过加exact选项要求输出精确的概率。 大样本时可直接采用正态近似的结果。,例3 多组独立样本资料的秩和检验,data sasbook； do a=1 to 4； do j=1 to 10； input x ； output； end； end； cards；,1.23 1.14 1.63 1.00 1.35 2.01 1.64 1.13 1.01 1.70 2.59 0.96 3.96 0.23 1.
10、61 2.94 1.96 3.68 0.54 2.59 2.50 3.92 3.33 0.23 2.07 3.90 2.84 2.42 3.72 3.00 0.40 0.51 2.49 0.74 2.59 2.36 3.00 3.31 3.21 3.12 ； proc sort； by a； run；,proc univariate normal； by a； var x； run； proc npar1way wilcoxon； class a； var x； run；,Wilcoxon Scores (Rank Sums) for Variable xClassified by Varia
11、ble aSum of Expected Std Dev Meana N Scores Under H0 Under H0 Score1 10 130.00 205.0 32.006610 13.0002 10 202.50 205.0 32.006610 20.2503 10 275.00 205.0 32.006610 27.5004 10 212.50 205.0 32.006610 21.250,Average scores were used for ties.Kruskal-Wallis TestChi-Square 7.7513DF 3Pr Chi-Square 0.0514,k
12、3,且每组例数均大于5，则此时不能查表， 且H值近似服从自由度为k-1的卡方分布。,例4：随机区组设计的秩和检验 data ex107; do b=1 to 8; do a=1 to 4; input x; output; end;end; cards; 8.4 9.6 9.8 11.7 11.6 12.7 11.8 12.0 9.4 9.1 10.4 9.8 9.8 8.7 9.9 12.0 8.3 8.0 8.6 8.6,8.6 9.8 9.6 10.6 8.9 9.0 10.6 11.4 8.3 8.2 8.5 10.8 ; proc freq; tables b*a*x/scores=
13、rank cmh2; /*得到FRIEDMAN统计量*/ run;,Summary Statistics for a by xControlling for bCochran-Mantel-Haenszel Statistics (Based on Rank Scores)Statistic Alternative Hypothesis DF Value Prob 1 Nonzero Correlation 1 14.1304 0.0002 2 Row Mean Scores Differ 3 14.5443 0.0023Total Sample Size = 32,作业： page204，第4、6、7题,
