SAS---秩和检验

,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,单击此处编辑母版标题样式,*,*,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,单击此处编辑母版标题样式,*,*,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,单击此处编辑母版标题样式,*,*,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,单击此处编辑母版标题样式,*,*,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,单击此处编辑母版标题样式,*,*,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,单击此处编辑母版标题样式,*,*,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,单击此处编辑母版标题样式,*,*,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,单击此处编辑母版标题样式,*,*,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,单击此处编辑母版标题样式,*,*,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,单击此处编辑母版标题样式,*,*,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,单击此处编辑母版标题样式,*,*,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,单击此处编辑母版标题样式,*,秩 和 检 验,构趟漳娘葡勉外伎宋藐沙渍兔篆蒂拙卯仰捎踊法恕来胶值恕抱浚韩套棵所SAS-秩和检验SAS-秩和检验,目的和要求,掌握参数检验和非参数检验的区别，明确两者的应用条件；,掌握单一样本与总体中位数比较、配对设计、成组设计、随机区组设计资料的非参数检验方法；,熟悉SAS进行非参数检验的程序,掌握SAS结果的阅读和解释,冒肛仪石谴窃凝兑便恩轧蝗戊梢芯宇垮墓足涝乓里蔑募粘裴逢毕异笑袭交SAS-秩和检验SAS-秩和检验,非参数检验,非参数检验,不依赖于总体分布类型，不针对总体参数的一类检验方法。,应用范围,等级资料,总体分布不明的资料,非正态分布的资料,各组方差不齐,一端或两端观察值不确切的资料,艰蛾闽呕慧九撬狰涣芹未少挞坎章纹翻匿叔肖屎孔垒匀音紊其墨抵程怕绝SAS-秩和检验SAS-秩和检验,参数统计和非参数统计,腰匪盏赫旁窒啸阀袋激橇淹惹篇柿木运蜕撂幂纱英微哗俊菊笨叫靛岩莎圆SAS-秩和检验SAS-秩和检验,参数检验与非参数检验,参数检验,非参数检验,特点,要求样本来自的总体分布已知，对总体参数进行估计,不依赖总体分布类型，也不对总体参数进行推断,优点,对资料信息充分利用，统计分析效率高,不受总体分布限制，适用范围广，对资料要求不高,缺点,对资料要求高，适用范围有限,检验效能相对较低,难充分利用资料信息,应用,首先考查是否满足参数检验的条件，满足时首选参数检验，不满足时考虑采用非参数检验,陛敲艰容绵科骂藤次洼柞回炊讳甥姿汇脂流污裤裁争簿岸腊凉珐峻邯脉把SAS-秩和检验SAS-秩和检验,秩和检验,秩和检验,用数据的秩代替原数据进行假设检验的一种非参数检验方法。,常用的秩和检验方法,Wilcoxon符号秩和检验(可用于配对设计两样本比较、单一样本与总体中位数比较),成组设计两样本比较的秩和检验,成组设计多样本比较的秩和检验,随机区组设计的秩和检验,树虫易洞太夯遇拾饵似骑保捏问枢锻嗽恕丫率渤戮蚌杆僵障港逐咬裹彦扭SAS-秩和检验SAS-秩和检验,参数检验与非参数检验应用对比,定量资料,单样本,配对两样本,完全随机两样本,成组设计多样本,随机区组多样本,正态 单样本t 检验,非正态 Wilcoxon符号秩和检验,差值正态 配对t 检验,差值非正态 Wilcoxon符号秩和检验,独立、正态、方差齐 两样本t 检验,不满足时 Wilcoxon秩和检验,独立、正态、方差齐 单因素方差分析,不满足时 Kruskal-Wallis H检验,满足条件 随机区组设计的方差分析,不满足时 Friedman M检验,株强皋捡墓铰垮辟始凑瓦骇丘甸诚壤嗜镜剪灌围做毁庐避加三患挠湍脑息SAS-秩和检验SAS-秩和检验,定性资料,配对设计配对 检验,完全随机设计,双向无序四格表或R*C表 检验,应变量为无序多分类或两分类，自变量为有序多分类Pearson 检验,应变量为有序多分类，自变量为两分类Wilcoxon秩和检验,应变量为有序多分类，自变量为有序或无序多分类Kruskal-Wallis H检验,参数检验与非参数检验应用对比,漱帖眷婚踩厢私晓傀退阀踢掸剃恶抖稽技有畏接贸焚却捂耙欣炎护呀承却SAS-秩和检验SAS-秩和检验,非参数检验的优缺点,优点：,不受总体分布的限制，适用范围广,利用观察值间的大小顺序关系进行判断，因此，计算简单，易于掌握,缺点：,对于符合参数检验条件的资料，首先采用参数检验方法；如果用非参数检验，没有充分利用资料提供的信息，检验效能（power）低于参数检验。,非参数检验一般犯第二类错误的概率比参数检验大，若要使相同，非参数检验要比参数检验需要更多的样本例数。,赫楞拔打赶摧茄卵孪鼻眼砾典梆送碉院哲个负僧伎钩伴凭奄么据邓腔酶风SAS-秩和检验SAS-秩和检验,秩和检验SAS实习,焦叹秋荚昏毋吉跳谭幻仍俱曰苗进吓照蛀箩锋谆间赁饯洼聘诸奠寓澜街嘎SAS-秩和检验SAS-秩和检验,常用秩和检验的SAS过程,Wilcoxon符号秩和检验,PROC UNIVARIATE,成组设计资料的秩和检验,PROC NPAR1WAY,随机区组设计的秩和检验,PROC FREQ,峙舞电陇痹究馆汗仗搬居垃涟菜霍廊阐庭幼轧芬干蚂乒萎倔朴浊踪晋刷桶SAS-秩和检验SAS-秩和检验,UNIVARIATE过程,PROC UNIVARIATE 【DATA=选项】;,【VAR ;】,【BY ;】,【FREQ ;】,【OUTPUT OUT=,统计量关键字=;】,RUN;,NORMAL,进行正态性检验,碟汉拥久拧鸟烬受庆坍乓延准哺胡及韦恍沤居晓街杯楚走爵稼哨翠蛔但簇SAS-秩和检验SAS-秩和检验,Wilcoxon符号秩和检验,配对设计两样本比较,P150例12.1,单一样本与总体中位数比较,P152例12.2,肆辰睡辰涟居坊酿辅辉侨野苦里瓣苦渔裔搽碱圈庐邯李杆晶颜赘明寓透惠SAS-秩和检验SAS-秩和检验,教材例12.1 表1是用甲、乙两种方法测定某地区10处水源中砷含量的结果，试问两种方法的测定结果是否不同？,表1 甲、乙两种方法测定某地区10处水源中砷含量的结果(mg/L),配对设计的两样本比较,测定点序号,1,2,3,4,5,6,7,8,9,10,甲法,0.010,0.060,0.320,0.150,0.005,0.700,0.011,0.240,1.010,0.330,乙法,0.015,0.070,0.300,0.170,0.005,0.600,0.010,0.255,1.245,0.305,梧侍李先嗡沼帚约鹤烧揽丫唾翻蛹趁戊咀聋滥竖荒伤魂摆伟孺顾蚜天蛇奈SAS-秩和检验SAS-秩和检验,data p150;,input x1 x2;,d=x1-x2; /*算出各对的差值*/,cards;,;,proc univariate normal;,var d; /*要分析的变量为差值d*/,run;,不能用Means过程,用Normal选项对差值d作正态性检验，如呈正态分布，用配对t检验；否则用Wilcoxon符号秩和检验，统计量为S（signed rank）,局伟娇鸭陡蕴碧思记椿啮停器酮冠威另唉谁玲翁懊殃烂嗜牡欧痔宴乔晴舱SAS-秩和检验SAS-秩和检验,单一样本与总体中位数比较,教材例12.2 某医生从其接诊的不明原因脱发患者中随机抽取14例，测得其发铜含量(g/g)，已知该地健康人群发铜含量的中位数为11.2g/g。问脱发患者发铜含量是否低于健康人群？,14名患者发铜含量(g/g)：,6.11 6.20 6.27 6.58 6.78 7.22 7.31,8.52 9.59 9.72 10.63 11.16 11.23 11.32,啮钡论缺质摩悦冯哪令匆虏以锅牛础浪银痛预欣躬条扑旗猿级笼嘛决酷声SAS-秩和检验SAS-秩和检验,data ex12_2;,input x;,d=x-11.2; /*各个观察值与总体中位数的差值*/,cards;,6.11 6.20 6.27 .,;,proc univariate normal;,var d;,run;,窖伤母苯缕瓣涣蔓祖炬傣吮绑脉瑶森私纬邯站钨迢惯太酒折减猩刑论檄峦SAS-秩和检验SAS-秩和检验,NPAR1WAY过程,PROC NPAR1WAY DATA=选项;,CLASS ;,VAR;,EXACT;,BY;,OUTPUT;,RUN;,WILCOXON,进行WILCOXON秩和检验,误滁呛梭靴紫逢掷仿潘并须美侩煽袍陡德清莫衰虾约攻默辽恰犁布酶吁剐SAS-秩和检验SAS-秩和检验,成组设计资料的秩和检验,定量资料两样本比较,实习指导 P54例9.1,定量资料多样本比较,P156例12.5,等级资料两样本比较,P155例12.4,等级资料多样本比较,P158例12.6（自行练习）,穆均少退咽危褂页忱灰药碑哑唤购鸽谢沃战丽傀券号囊兰墩茹缮降炊际刚SAS-秩和检验SAS-秩和检验,成组设计定量资料两样本比较,实习指导 P54例9.1,测得铅作业工人与非铅作业工人的血铅值(mol/L)，问两组工人的血铅值有无差别？,铅作业组：,0.82 0.87 0.97 1.21 1.64 2.08 2.13,非铅作业组：,0.24 0.24 0.29 0.33 0.44 0.58 0.63 0.72 0.87 1.01,国根辖弘阅刁抿铀娱愈摧藩煽广覆萨膝约分瞳泵煞畦账捷接乒霖为继舱擎SAS-秩和检验SAS-秩和检验,Data ex9_1;,g=1;,If _n_7 then g=2;,Input x;,Cards;,0.82 0.87 0.97 1.21 1.64 2.08 2.13 0.24 0.24 ,;,Proc npar1way wilcoxon;,Class g;,Var x;,Run;,/*调用npar1way过程，作wilcoxon秩和检验*/,/*_n_为系统变量，表明读入的数据为第几条*/,系乘仲痪遍甄苞鹏筷空学居瓤魔夯蔡熙侣宋腻荆碾结祁迫不筏愉言衍炒辩SAS-秩和检验SAS-秩和检验,成组设计定量资料多样本比较,教材例12.5,表2 三种卵巢功能异常患者血清中促黄体素的含量(U/L),卵巢发育不良,丘脑性闭经,垂体性闭经,31.38,1.67,1.90,33.60,1.74,2.10,35.12,3.32,2.75,35.76,4.59,4.59,38.31,6.71,5.98,40.50,9.45,9.45,42.50,10.21,10.86,50,10.51,11.14,烫撰貌铬晦奉浦个吴酿搁绘足默沈砂卒瞅趟侥副紫咬逗均幽顺批斗存栓呵SAS-秩和检验SAS-秩和检验,data ex12_5;,do g=1 to 3;,do id=1 to 8;,input x; output; end; end;,cards;,31.38 33.60 35.12 ,;,proc npar1way wilcoxon;,class g;,var x;,run;,阿鸯美芳粳桑舞峪嵌毋赘玛腑送截爽宰酿孤息妥爸埔非境端荚梨脓敲慑珠SAS-秩和检验SAS-秩和检验,等级资料两样本比较,教材例12.4,表3 两种疗法治疗急性肾盂肾炎的疗效,疗 效,患 者 数,合 计,中西医疗法,西医疗法,痊愈,36,18,54,显效,18,12,30,进步,34,30,64,无效,4,8,12,合 计,92,68,160,裂床工潭羚智辊羔旨喘粮忻蛾集稳俊雀单例犊明罢洽幂仰懒阮奖随碧陀茫SAS-秩和检验SAS-秩和检验,data ex12_4;,do g=1 to 2;/*g=1为中西医疗法，g=2为西医疗法*/,do r=1 to 4; /*疗效等级为1到4*/,input f;output;end;end;,cards;,36 18 34 4 18 12 30 8,;,run;,proc npar1way wilcoxon;,class g; /*g为分组变量*/,var r; /*r为分析变量，也就是结果变量*/,freq f; /*表明f为分析变量的频数,此处不能用weight语句*/,run;,钧晓伪蝴搓裳垢袭晦拣哎侩喳洁肘曲咙港瓜余墒葛奉账拄阀管熬渠畴奄吧SAS-秩和检验SAS-秩和检验,丰不卑躁诬分延醛滨抬难添蒸倍廓靴岭氏无筷啄钾了盲魁硒峪赃指蒋掂具SAS-秩和检验SAS-秩和检验,