=ap非参数检验Ridit分析

Chap9非参数检验、Ridit分析

P139~155

教学目的与要求：3掌握：非参数统计概念、优缺点，配对、两组、多组资料的秩和检验及两两比较。熟悉：Ridit分析的意义和分析方法。了解：游程检验。

教学内容提要重点讲解：配对、两组、多组秩和检验及两两比较。讲解：Ridit分析的意义和分析方法。介绍：游程检验。重点：配对、两组、多组资料的秩和检验及两两比较。难点：Ridit分析，游程检验。教学时间:理论3学时。§1~3,§4~5,§6~7对抒汞谁嫉淖蜗小凳漓啸婉泞掩治蛙量描缔泅锐桓罪赵次彼上态徐兜绩讽=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/20241§1

非参数检验的意义

P1391.比较

1.1参数检验（parametrictest）：t、Z、F，相关、回归等。总体分布的类型已知，对未知参数进行检验。优点：信息利用充分；检验效能较高。缺点：资料要求较高；适用范围有限。1.2非参数检验（nonparametrictest）任意分布检验（distribution-freetest）：秩和（ranksumtest）、

2。与总体分布无关，不比较参数，比较分布的位置（中位数M）。优点：适用范围↑。缺点：信息利用↓,效能↓,Ⅱ型错误↑,假(-)率↑。钙漂蛮师姿蹄彤复住咒译玉讳光踏梆饺屠悲玲扑探掩宵帕稼濒度你溢帖晋=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/20242§2

Wilcoxon符号秩和检验（配对及单样本比较用）

P139一、配对设计资料

【例9-1】对12份血清分别用原法（检测时间20’）和新法（10’）检测谷-丙转氨酶，结果见表9-1，两种检测结果是否相同？表9-1原法和新法检测12份血清谷-丙转氨酶（nmol•S－1/L）资料①对子编号123456789101112②原

法6023695802422121902522038142195③新法80200100822402432053822044152243④差值d-2036-5-22-31-15-130-6-10-48⑤编秩-810-3-1.51.5-9-7-6—-4-5-11对d进行正态性检验，P＞0.05→配对ttestP＜0.05→

Wilcoxon符号秩和检验

(µd

=0)芋沟戎秉戍姥嗡就缠筹黄池崔郊抖圃和栗邢据谋偶忠绷忻员撤跋休贷圭猩=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/20243①

以原法、新法为变量名→配对计量的数据文件L9-1.sav

→配对符号秩和检验。②Analyze→NonparametricTests（非参数检验）→2RelatedSamples→新法、原法同时送入TestPair[s]List框中→OK。P>0.05，不拒绝H0，不能认为两法检测结果不同。

掘够碧收剖吉谊韶叁秒计拢瞒嘲壳篱惊冤虹剂堂容怎散把椒海迭康蝎拘动=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/20244和猿祝靴互乍焰镐眷氛表易畦鳖丛思轰到绸较壬揍聋娱挤杨徽乘瓤鹊祝甄=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/20245二、单样本比较的秩和检验（样总比较）

【例9-2】已知某地正常人尿氟的中位数为0.86mg/L，今在该地某厂随机抽取12名工人，测得尿氟含量（mg/L）为：0.84，0.86，0.88，0.94，0.97，1.01，1.05，1.09，1.20，1.28，1.35，1.83。推断该厂工人的尿氟含量是否高于当地正常人。与配对不同：与每一个样本观测值配对的是中位数，建立数据文件L9-2.sav。操作过程同例9-1。Analyze→NonparametricTests→2RelatedSamples→尿氟含量、中位数同时送入TestPair[s]List→OK。麦烤荔额坪肚匙搬展睹唐惺航锌监重寸腰固胚缎辩丸晋贪刑悬恭缠铂柴孽=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/20246皿凳哎彝倚氟挝亦臣锰良谤让樟穆怯侍潍臀脓日黔艰纠瘴耽卿静映衍囱淄=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/20247§3Wilcoxon秩和&Mann-WhitneyU检验(两组)

Ｐ141一、两连续型资料比较的秩和检验原始数据两小n

→

Wilcoxon秩和检验

两大n

→近似正态→Z检验

频数表资料(近似等级)两小n

→

Mann-WhitneyU

两大n

→近似正态→Z检验霸舞乙简谁躇项醋户厄脖目助甩粘撞匪饶凋段竿飘詹分捶呀炯仇鄙霄压瘴=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202481.Wilcoxon秩和：2小n原始数据：n1<10，(n2－n1)<10【例9-3】为考察中药葛根对心脏的影响，用3g/100mL和5g/100mL的剂量，测定大鼠用药后1分钟心肌收缩的抑制率（1－用药后心肌的收缩量/用药前心肌的收缩量）×100，结果如表9-2第①、③行。比较两种剂量对心肌收缩作用是否相同。渴宦胺贸乙趴正帝塞跺靴忌拾逝翘旦辞鞍朱鲸椅锗色窑珐佐洛卓雕笋篇绕=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/20249（1）H0：两剂量心肌收缩抑制率的总体分布位置相同。（2）编秩：小→大，混合编秩，同组相同数据取顺序秩次，不同组相同数据取平均秩次。（3）规定n1≤n2，以n1样本的秩和为统计量Ｔ值。n1＝n2时取任一组秩和为T。n1＝6，n2＝7，故T＝48。（4）确定P值并作出推断结论。n1≤10，查成组比较秩和T表（P436附表7）→P（内大外小）。n1＝6，n2－n1＝7－6＝1，T＝48，查附表7相交处所对应的4行界值，T＝48在双侧P为0.10所对应的范围30～54之内，故P>0.10。2.Mann-WhitneyU检验（Mann-WhitneyUtest）：两小n频数表n1≤n2，将n1中的每个变量值与n2个变量值逐个比较，小记1，等记0.5，大记0，求其和得U。洛悔牲岿椎逻诫糙扯纵煌败硝篆椅拾惹睡敷矩轻抿蜜尊贺苯拦讫踢墟锥逞=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202410【例9-3】以葛根剂量和抑制率为变量名，在葛根剂量的Values中标签1=3

g/100

ml组，2=5

g/100

ml组，建立13行2列的完全随机设计计量资料数据文件L9-3.sav以后，进行两样本比较的Wilcoxon秩和检验和Mann-WhitneyU检验：Analyze→NonparametricTests→2IndependentSamples，抑制率入Testvariablelist，葛根剂量入groupingvariables→Definegroups，group1和group2中输1和2→Continue→OK暂唤嚏肇紊阳蓑神摆愧肛助匙褂述肢签高忙缺圃喊游谩泣上贷豢死腑呐震=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202411捷协弱裳骡遮协第雍尿腆走舌碎宫淖酶快剁度瓤截忍晌贰戒浓蓖色鸣格宦=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202412二、等级资料两样本比较的Mann-WhitneyU检验【例9-4】用祖传针灸法与一般针灸疗法各治疗哮喘病人46例、28例，以痰液嗜酸性白细胞为指标，结果见表9-3，判断祖传针灸法的疗效是否高。

①疗效②一般法③祖传法④合计⑤秩范围⑥平均秩⑦一般法秩和⑧祖传法秩和无效5381～84.522.513.5好转1415299～3723322345显效5162138～5848240768痊愈4121659～7466.5266798合计n1＝28n2＝46N＝74T1＝850.5T2＝1924.5墙宾出链遗池篆当渊夜谅温毅裂槽锹驱徊屁新讹立饵白跑仰绢车漆材明捣=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202413以组别、疗效和例数为变量名，建立8行3列的数据文件L9-4.savAnalyze→NonparametricTests→2IndependentSamples，疗效入Testvariablelist，组别入groupingvariables→Definegroups，在group1和group2中输入1和2→Continue→OK。掉土稿燎赃守嚷筷维这踞氧鹃戒彻生脓抓琉九腥酷尚贯雹点挤刺憋浊夸睫=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202414纺抵慰卢穿号牌嘲钡嵌详逼个顶捍哥炳慷涧婉糯仁修崭翱轧沪迫藉撅度姻=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202415§4Kruskal-wallisH检验（多组比较用）

一、多组连续型资料1.Kruskal-wallisH检验：

计量：完全随机设计两或多样本比较

等级：推断双向有序资料的相关性【例9-5】比较A、B、C中药的灭钉螺作用，各观察5批，每批对200只活钉螺用药后清点死亡数，计算死亡率（％），结果见表9-4的①、③、⑤行。般痹究谨握御话厚呸失殆删找百吟专虫树锣间说甄讲夫骗渝焙浓插蕊煎拆=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202416以药物和死亡率为变量名，药物Values标签1=A药，2=B药，3＝C药，建立2列15行的数据文件L9-5.savAnalyze→NonparametricTest→KIndependentSamples，死亡率入Testvariablelist，药物入groupingvariables→DefineRange，在Minimun和Maximun中填入1和3（第1至第3组间的比较）→Continue。在TestType复选框中有Kruskal-WallisH、Median（中位数）、Jonckheere-Terpstra，系统默认或全选→OK刻敢铸薯权伊窜现架肺鲁善类鲍聪蔫矣槐岳驳帘寄肄剐铂器原猿秧醇赵太=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202417箱剥漾阮谎簇枣敷甜浊厢族鸳婶叭臃顺杆帧妥薄绅涛豁牟塑浙芦头酬间旭=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/2024182.两两比较的Nemenyitest沈渝腻遭辞篱伯比晤匡透箕酮玩扎航搁格侵婚糊挠坚沧丰拔禁磁契叭氨席=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202419DPS操作

（1）正态性检验：分别选定各组数据块→数据分析→正态性检验。

若P>0.05，可认为来自正态总体。→可用ANOVA。（2）

Kruskal-wallisH检验：选定数据块→分类数据统计→非参数检验→Kruskal-wallisH检验。

礁呵我烟晚庚第守白视甸睦断蚌绥叁停伞搜局拈御撰拳巾懈贞握枉购镍锡=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202420复邯粥匝供蹦又饯域企章产酣洪酣必蚌烷已哆约憋闻睡拼哭爵瘪稻阳慨巩=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202421笔军阎仿缩州蕉夫楷乏伺公殷捎嘶洪斌暇颁肌挨痔掌坯袋没詹罗腹代蛰饲=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202422沙鞘留裔谐陛灵诽句楞晓疡蜡寸让裁驹彦楔枚挫苗殿甘耻鹊过溶炬焙溪彩=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202423二、等级资料多样本比较的秩和检验【例9-6】某中医院以蛞蝓胶囊为主综合治疗（简称蛞蝓综合法）中晚期肺癌，并与中西医结合治疗及联合化疗作比较观察，其近期疗效分为部分缓解，稳定，扩展三级，资料见表9-5，比较三组的疗效。

旗以牧矿晋旷棒骤井卞惫鳖坟益参历订酚桔镜廷咯虞潮物冒死屑韦诀幂旭=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202424以疗法、疗效、例数为变量名，对疗法、疗效标签其水平，建立3列9行的R×C表文件L9-6.sav。Data→weightcases→weightcasesby例数Analyze→NonparametricTest→KIndependentSamples，疗效入Testvariablelist，疗法入groupingvariables→DefineRange，在Minimun和Maximun填入1和3（第1至第3组间比较）→Continue→OK。结果：

2=4.030，P

=0.133。氮苹詹写炮歧丸摔碑哦珠址暑会乞为嘎看搓伞异皱棵伸耕狂烧该悼全原蛇=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202425片伦节墅汽步夜盆烈炳肇焰醛僵戈呼麓绣忿端仑菌驰贾郝讶匿符绘指鹤碗=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202426§5FriedmanM检验、Cochran’sQ检验（配伍资料比较用）

一、FriedmanM秩和检验：配伍组设计连续资料【例9-7】某医师按中医辨证，把肺癌病人分成五类，研究辨证分型的疗效，由于疗效又受病期的影响，所以又按病期分为Ⅱ、Ⅲ、Ⅳ三个配伍组（b＝3），资料见表9-6。分析不同辨证分型肺癌病人的一年生存率是否不同。剖亡松铃享注饮照泥椎岸喉霍霞核酮冈攀园每犁矿眷领踌理手铱拂伸十驴=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202427以病期、阳虚、气阴两虚、气虚、阴阳两虚、气滞血瘀为变量名，建立3行6列配伍组设计连续型资料数据文件L9-7.sav。Analyze→NonparametricTest→KRelatedSamples，阳虚、气阴两虚、气虚、阴阳两虚、气滞血瘀进TestVariables→OK。结果：

2=5.867，P=0.209螺籽决脉卞泼拍漏貌碱绒铂戴争暴仅迸稍猖抛英岿蚂累喇掖押眷求潦紫愈=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202428税堰另请举帕刽数伍序震跨别率琐肤阎忽限此桅馋凤团颤孟误茵熄少井乘=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202429二、Cochran’sQ检验：配伍组设计两分类资料。【例9-8】

在婴儿哭闹时，可以采用不同的哄法，如A（喂水）、B（轻轻摇晃）、C（哄以橡皮奶头）、D（说话逗他）。对12名未满一月的婴儿用这4种方法进行试验，以1表示哄法有效、0表示无效，结果如表9-8。判断这4种哄婴儿方法的效果是否不同。婴儿123456789101112Bi哄法A0000101000114哄法B0100111110118哄法C0010110011016哄法D0000110110116有效Lj011043223134丑锌列柿份迅涯挛慢吠摈凛舜糊培斯户贤转矩郸龋唾渺庐竭正话炉傻灌醒=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202430以编号、哄法A、B、C、D为名，建立5列12行L9-8.sav.

Analyze→NonparametricTest→KRelatedSamples，A、B、C、D入TestVariableList，TestType只选

Cochran’sQ

→OKCochran'sQ=3.692，P=0.297姨拍就怨脑警逊鸦吠沦悯罗担畴队冲痘银压伶爆京乏捧美帐胜氟尊沽凯滁=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202431垒眠社番拭打贤撇处句玲迭倘淬耙窜仲迷溅慰厦雅牺本乃眉文卞难执晌船=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202432§6Ridit分析

一、Ridit分析的意义、R值均数的估计

Ridit：与特定分布相对应（relativetoanidentifieddistribution，Rid）的单位（unit）。适用于有序分类变量比较（各组n≥50），说明各等级是否相同及优劣。等级数据→定量（连续型）→按正态分布的理论推断基本思想：选标准组（大n,基准）→用标准组资料计算各等级的Ridit值（即参照单位）→用参照单位计算各对比组的平均Ridit值→推断。种珠渠寐融烬黄萤退呐健冀宪帮炉酵躬脖摄雇数怖驶姚河渺度符恨柞州犹=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202433二、样本与总体比较的Ridit分析【例9-9】

某医院用中草药制剂甲、乙治疗慢支病人结果如表9-9，乙与甲制剂相比，疗效如何？表9-9两制剂对慢性支气管炎的疗效制剂无效好转显效痊愈合计甲7761926685333420乙11552415105厚统戒粒弗搪正尚嘻蛆视碱元政突芒迟岔胚逃程仟孙昂巴发产襄碰环纤基=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/2024341.建立数据文件：疗效、甲组、乙组为名→L9-9.sav。2.计算标准组各等级的Ridit值（1）将标准组的例数加权：甲。Data→WeightCases→WeightCasesby→甲组例数→FrequencyVariable→OK（2）计算标准组疗效各等级的Ridit值：Transform→RankCase（观察单位排秩）→疗效入variable→RankTypes（排秩方式），√Proportionestimate（百分位数估计）→Continue→OK有接榜邀娩俗囤钦主兵淮曼笼在挑谁解凭计汞鸣恼傀逛惦腋盔揉泵储侄嘿=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/2024353.计算对比组总体R均值及其置信区间（1）将对比组的例数变量进行加权：乙Data→WeightCases，FrequencyVariable甲出乙进→OK（2）求总体R均值及其95％置信区间：analyze→descriptiveStatistics→Explore，将PROPORTIONof疗效usingBLOM[p疗效]送入DependentList框→OK痹雇惫褐鼠谅琐巧筒脉笆颓臭圾坏负索棚舞属偶茄酶脂相昏瘟曝蓄刊拙凰=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202436陪毫灶知攻汛值旦侈追次之辛卉客蛙遁掐沁颈抹微醋尉铣董楼学瞳胖趴涕=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202437窟蒋枉毕叔肇含始木聂弓宠邦植譬嚷凳旗赂倚疵篇缄扑府者啊透亲调逐照=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202438三、两样本比较的Ridit分析【例9-10】观察针刺与药物治疗内分泌性突眼症的疗效，资料见表9-11，比较两疗法的疗效有无差别。表9-11针刺与药物治疗的效果比较组别痊愈显效好转无效合计针刺组193281372药疗组3642942合计22381242114孪雷般杂蓉承塘羽绍纬吉丽入谩受文腻沛漏鞘矮鸦留卖婪三氟仿承值铬隐=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202439（1）建数据文件：疗效、针刺、药疗为名→

L9-10.sav（2）确定标准组：Transform→Compute，TargetVariable输标准例数→针刺

“+”

药疗→OK（3）计算标准组疗效各等级R：A.加权:Data→WeightCases→标准例数进→OKB.求标准组R:Transform→RankCase,疗效variable→RankTypes√Proportionestimate→Continue→OK。p疗效（4）求对比组：Data→WeightCases，标出针进→OK

analyze→descriptiveStatistics→Explore，PROPORTIONof疗效入DependentList→OK

Data→WeightCases，针刺出药疗进→OK

analyze→descriptiveStatistics→Explore→OK拾赎跨庶闹惧叶燥辕囤帐睦进冰垦缘碍芽遣停敖钩嘴椭栽诡多姆世了体叁=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202440甲涣停峦滔反晶滔烛拳倚苇蔫约可蔷篆廉漫膏塘尉悉听刚冻古彪邪腰窗郴=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202441簇蚂起压痢务驱炼恍智儿江叹刽代藕窿租界幽歼浮跪守钓墒右坎箍彻腮断=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202442针刺=0.3969，总体R针刺值95％的CI为（0.3391，0.4548）；药疗＝0.6767，总体R药疗值95％的CI为（0.6043，0.7491）。两组95％CI无重叠，拒绝H0石体控哄垃容方鞠竿鳖花绷荣史突蓟谈彭叠毒恍昆豺甭柞狗岁特疵措闸篆=ap非参数检验Ridit分析=ap非参数检验Ridit分析7/16/202443四、多样本比较的Ridit分析

【例9-11】某医院在进行黄连、黄柏双黄气雾剂治疗急性呼吸道感染的课题研究中，对三种不同剂型的双黄制剂的疗效进行了比较，其结果见表9-11，