1、秩和检验（二）,等级资料的秩和检验 多组独立样本的秩和检验,PROC NPAR1WAY Wilcoxon ；CLASS 变量；VAR 变量表；RUN；,SAS中用于秩和检验的过程,某医科大学营养教研室为了解居民体内核黄素营养状况，于某年夏冬两个季节收集成年居民口服5mg核黄素后4小时的负荷尿，测定体内核黄素含量，结果见表10-4第(1)、(2)、(3)栏，试比较该地居民夏冬两个季节体内核黄素含量有无差别？,例1 等级资料,研究指标“核黄素的营养状况”是一个等级指标，不适宜用卡方检验。,data ex3; do RIB=1 to 3; *rib代表核黄素的状况; do season=1 to 2
2、; *season达标不同季节; input count; output;end;end; cards; 10 22 14 18 16 4 proc npar1way wilcoxon; class season; *按季节分组; freq count; *freq语句指定有序资料中各个等级的频数; var RIB; *分析的指标就是核黄素的状况rib; run;,程序,结果,秩和,结果,T=2036，P=0.0019，按照0.05水准，拒绝H0， 故可认为夏冬两个季节居民体内核黄素含量有差别。,例2 多组独立样本资料的秩和检验 比较四组的降胆固醇效果是否相同。,data ex2； do a=
3、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.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 npar1way wilcoxon； class a； var x； ru
4、n；,程序,Wilcoxon Scores (Rank Sums) for Variable xClassified by Variable 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,k3,且每组例数均大于5，则此时不能查表，且H值近似服从自由度为k-1的卡方分布。故本例2=7.75，P=0.05140.05，不拒绝H0，尚不能认为四组的降血脂效果有差异。,为研究西布曲明片减肥效果，将BMI大于等于25的193名成年女性分成两组，实验组采用服用西布曲明片和饮食疗治疗，对照组仅采用饮食治疗，结果见表9-16，请问西布曲明片有无减肥作用？,作业,
