最新卫生统计学实验课件秩和检验Nparway过程精品课件.ppt

目的要求掌握利用掌握利用univariate过程完成过程完成配对设计配对设计资资料的非参数检验；料的非参数检验；掌握利用掌握利用npar1way过程过程及选择项及选择项Wilcoxon的使用完成完全随机设计两组、的使用完成完全随机设计两组、多组资料的秩和检验；多组资料的秩和检验；熟悉利用熟悉利用freq过程过程完成完成随机区组设计随机区组设计资料资料的秩和检验且能对结果作分析。的秩和检验且能对结果作分析。 单样本秩和检验data a;input nf;cards;2.15 2.10 2.20 2.12 2.42 2.522.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.15 Test -Statistic- -p Value- Students t t 2.898449 Pr |t| 0.0145 Sign M 3.5 Pr = |M| 0.0654Signed Rank S 29.5 Pr = |S| 0.0059 Tests for NormalityTest -Statistic- -p Value-Shapiro-Wilk W 0.837966 Pr D 0.1500Cramer-von Mises W-Sq 0.092826 Pr W-Sq 0.1281Anderson-Darling A-Sq 0.651796 Pr A-Sq 0.0704 例例2 为比较两种刺激物的对人的影响力为比较两种刺激物的对人的影响力的持续时间，现随机抽取的持续时间，现随机抽取13个试验对象接个试验对象接受第一种刺激物，受第一种刺激物，6个试验对象接受第二个试验对象接受第二种刺激物。持续时间如下表种刺激物。持续时间如下表，问，问二种刺激二种刺激物的持续时间是否有差异？物的持续时间是否有差异？两独立样本秩和检验两独立样本秩和检验 持续时间表(minutes) 刺激物1组 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 刺激物2组 3.27 3.27 3.27 3.70 3.70 3.74 PROC NPAR1WAYPROC NPAR1WAY过程格式过程格式 PROC NPAR1WAY WilcoxonWilcoxon ； CLASS 变量； VAR 变量表； RUN；data sas1data sas1；do a=1 to 2do a=1 to 2；input ninput n；do j=1 to ndo j=1 to n；input xinput x；outputoutput； endend； endend；cardscards； 13131.94 1.94 2.92 2.92 2.92 2.92 3.27 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.743.27 3.27 3.27 3.70 3.70 3.746 63.27 3.27 3.27 3.70 3.70 3.743.27 3.27 3.27 3.70 3.70 3.74； proc sortproc sort；排序过程；排序过程by aby a；runrun；proc univariate proc univariate normalnormal；考虑各组的正态性；考虑各组的正态性var xvar x；by aby a；runrun；proc nparproc npar1 1way way wilcoxonwilcoxon；class aclass a；var xvar x；exactexact；输出精确的概率；输出精确的概率run； 结果 The NPAR1WAY Procedure Wilcoxon Scores (Rank Sums) for Variable x Classified by Variable a Sum of Expected Std Dev Meana N Scores Under H0 Under H0 Score-1 13 110.50 130.0 11.004784 8.5002 6 79.50 60.0 11.004784 13.250 每组得分和每组得分和 H0成立时的得分和成立时的得分和 H0成立时的标准差成立时的标准差 平均得分平均得分 Average scores were used for ties. Wilcoxon Two-Sample Test Statistic (S) 79.5000 Normal Approximation Z 1.7265 One-Sided Pr Z 0.0421 Two-Sided Pr |Z| 0.0843 t Approximation One-Sided Pr Z 0.0507 Two-Sided Pr |Z| 0.1014 T统计量统计量正态近似法正态近似法结果（超出结果（超出表范围时使表范围时使用）用）近似近似t检验结果检验结果 Exact Test One-Sided Pr = S 0.0527 Two-Sided Pr = |S - Mean| 0.1054小样本时可通过加小样本时可通过加exact选项要求输出精确的概率。选项要求输出精确的概率。大样本时可直接采用正态近似的结果。大样本时可直接采用正态近似的结果。例例3 多组独立样本资料的秩和检验多组独立样本资料的秩和检验data sasbookdata sasbook；do a=1 to 4do a=1 to 4；do j=1 to 10do j=1 to 10；input x input x ；outputoutput； endend； endend；cardscards；1.23 1.14 1.63 1.00 1.35 2.01 1.64 1.13 1.01 1.701.23 1.14 1.63 1.00 1.35 2.01 1.64 1.13 1.01 1.702.59 0.96 3.96 0.23 1.61 2.94 1.96 3.68 0.54 2.59 2.59 0.96 3.96 0.23 1.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.002.50 3.92 3.33 0.23 2.07 3.90 2.84 2.42 3.72 3.000.40 0.51 2.49 0.74 2.59 2.36 3.00 3.31 3.21 0.40 0.51 2.49 0.74 2.59 2.36 3.00 3.31 3.21 3.123.12；proc sortproc sort；by aby a；run； proc univariate proc univariate normalnormal；by aby a；var xvar x；runrun；proc nparproc npar1 1way way wilcoxonwilcoxon；class aclass a；var xvar x；run； Wilcoxon Scores (Rank Sums) for Variable x Classified by Variable a Sum of Expected Std Dev Mean a N Scores Under H0 Under H0 Score 1 10 130.00 205.0 32.006610 13.000 2 10 202.50 205.0 32.006610 20.250 3 10 275.00 205.0 32.006610 27.500 4 10 212.50 205.0 32.006610 21.250 Average scores were used for ties. Kruskal-Wallis Test Chi-Square 7.7513 DF 3 Pr Chi-Square 0.0514k3,且每组例数均大于且每组例数均大于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.711.6 12.7 11.8 12.0 9.4 9.1 10.4 9.89.8 8.7 9.9 12.08.3 8.0 8.6 8.68.6 9.8 9.6 10.6 8.9 9.0 10.6 11.48.3 8.2 8.5 10.8;proc freq;tables b*a*x/scores=rank cmh2; /*得到得到FRIEDMAN统统计量计量*/ run; Summary Statistics for a by x Controlling for bCochran-Mantel-Haenszel Statistics (Based on Rank Scores)Statistic Alternative Hypothesis DF Value Prob1 Nonzero Correlation 1 14.1304 0.00022 Row Mean Scores Differ 3 14.5443 0.0023 Total Sample Size = 32作业： page204，第4、6、7题