【文档说明】协和医学院课件-非参数检验.ppt，共(28)页，799.001 KB，由小橙橙上传

转载请保留链接：https://www.ichengzhen.cn/view-254605.html

以下为本文档部分文字说明：

2023/5/231非参数检验基础所流行病王艳红流行病学资料分析2一、相关理论知识的回顾非参数检验特点(定期复习十分重要)总体分布不作严格规定；不依赖于总体分布类型（如正态分布）——任意分布检验非参数检验的适用条件总体的分布不易确定偏态分布方差不齐等级资料（统计资料一般分为计量资料与计数

资料两大类，介于其间的一类资料称为等级资料.不同类型的资料应采用不同的分析方法。例如用某药治1.5指标的特异性指测量特定项目而不疗若干痢疾患者,其中治愈、好转、无效人受其它项目干扰的能力文献来源；化验结果分为:-、+、+

+、+++、++++5个有序等级）3一、相关理论知识的回顾基于秩次的非参数检验4一、相关理论知识的回顾（1）、Wilcoxon秩和检验应用•两组连续变量资料的秩和检验•两组有序变量资料的秩和检验基本思想：•假设两总体的分布相同（H0），两样本可认为是从同一

总体中抽取的随机样本；将二者混合后由小到大编秩，然后分别计算两组的平均秩和R1与R2，R1和R2如果近似相等，其差别可能是由随机抽样引起的。5一、相关理论知识的回顾6一、相关理论知识的回顾（1）、Wilcoxon秩和检验公式1——没有节点（即没有相同的观察值）公式2—

—有节点()()112112121122112nnnRTnnnn++−−=++()()112121121212121122(1)112()(1)giiinnnRTttnnnnnnnn=++−−=−++−++−7一、相关理论知识的回顾

（2）、Kruskal-Wallis秩和检验应用•多组连续变量资料的秩和检验•多组有序变量资料的秩和检验基本思想•同Wilcoxon秩和检验8一、相关理论知识的回顾（2）、Kruskal-Wallis秩和检验，当K=3且至少有一组例数大于5时，或者k>3时服从卡方分

布；若K=3且没有一组例数大于5时，服从H分布。公式1——没有结公式2——有结()21123(1)1kiiiRHNNNn==−++()()21313123(1)11kiiigjjiRNNNnHttNN==−++=−−−

9一、相关理论知识的回顾（3）、Friedman秩和检验应用•随机区组设计资料基本思想•各区组的观察值按由小到大的顺序进行编秩；如果各处理的效应相同，各区组内秩1，2，……，k应以相同的概率出现在各处理

组（列）中，各处理组的秩和应该大致相等，不可能出现较大的差别。10（3）、Friedman秩和检验公式11二、相关SAS过程步（1）、NPAR1WAY过程两组连续变量资料的秩和检验多组连续变量资料的秩和检验（2）、RANK过程

+ANOVA(或GLM)过程完全随机设计连续变量资料的秩和检验随机区组设计资料的秩和检验多重比较（3）、Freq过程等级资料两组连续变量资料的秩和检验完全随机设计连续变量资料的秩和检验随机区组设计资料的秩和检验单因素单因素、双因素、两

两比较单因素、双因素、等级资料12二、相关SAS过程步（1）、NPAR1WAY过程的基本结构procnpar1waywilcoxon;classgroup（分组变量）;varx（连续变量）;run;Wilcoxon是指资料的分析基于Wilcoxon得分两样本资料：Wilcoxo

n秩和检验多样本作Kruskal-Wallis秩和检验13二、相关SAS过程步（2）、RANK过程+ANOVA(或GLM)过程基本结构完全随机设计（Kruskal-Wallis秩和检验）procrankdata=sam

pleout=a;varx;ranksr;/*产生秩次*/run;procglmdata=a;classgroup;modelr=group;meansgroup/snk;run;procsortdata=a;byblock;run;pr

ocrankdata=aout=b;varx;byblock;ranksr;run;procglmdata=b;classtreatblock;modelr=treatblock;meanstreat/snk;run;随机区组设计（Friedman秩和检验）14二、相关SAS过程步（3）、

Freq过程的基本结构（CMH解法）procfreqdata=aa;tablesgroup*x/scores=rankcmh;run;15三、常见的非参数检验实例（一）、配对设计和单样本资料的非参数检验（二）、完全随机设计两组独立样本的秩和检验（三）、完全随机设计多组独立样本的秩

和检验（四）、随机区组设计资料的秩和检验16（一）、配对资料的非参数检验•拟对某种药物的降压效果进行评价，治疗前后病人的舒张压水平如下：•试用非参数检验的方法比较治疗前后舒张压有无差异？datadbp;inputbeforeafter@@;d=before-after;cards;1171

2312710814112010710711010011498115102138152127104122107;run;procunivariatedata=dbpnormal;vard;run;17（二）、完全随机设计两组独立样本的秩和检验实例

•某缺碘地区测定9名12岁男孩和8名女孩的尿碘含量（μg/l），结果如下：–男孩（9名）：102.00、45.70、103.00、105.00、98.30、113.00、64.00、113.00、54.00–女孩（8名）：61.00、49.70、45.

80、33.90、116.00、118.00、120.00、72.00•试用非参数检验的方法比较不同性别的儿童尿碘含量是否不同?18SAS数据集的建立dataaa;if_n_<10thensex=1;if_n_>=10thensex=2;inputx@@;cards;102.004

5.70103.00105.0098.30113.0064.00113.0054.0061.0049.7045.8033.90116.00118.00120.0072.00;run;SAS过程步procnpar1waywilcoxondata=aa;varx;classsex;run

;procfreqdata=aa;tablessex*x/scores=rankcmh;run;行均值得分的差值Rowmeanscoresdiffer19练习1•应用免疫酶法对鼻咽癌患者和非鼻咽癌者分别测定13人的血清病毒VCA-LOGA抗体滴度（倍），

结果如下–鼻咽癌患者：5204080808016016032032064012801280–非鼻咽癌患者：510102020204040808080160640•试用非参数检验法比较两组患者的血清病毒的VCA-LOGA抗体滴度是

否有差异？20（三）、完全随机设计多组独立样本的秩和检验实例•某研究者欲研究A、B两个菌种对小鼠巨噬细胞吞噬功能的激活作用，将60只小鼠随机分为三组，其中生理盐水为对照组，用常规巨噬细胞吞噬功能的监测方法，获得三组吞噬指数如下：–对照组：1.301.401.501.501.601.601.7

01.701.701.701.702.102.302.302.302.402.40–A菌组：1.801.802.202.202.202.202.302.302.302.662.662.682.682.702.802.803.003.103.103.103.103.204.304.30–B菌组：1

.501.801.802.002.302.302.402.402.402.402.502.502.602.602.602.702.703.103.10•试用非参数检验的方法比较三组吞噬指数有无差别？21SAS数据

集的建立dataaa;if_n_<18thengroup=1;if_n_>=18and_n_<42thengroup=2;if_n_>=42thengroup=3;inputx@@;cards;1.30

1.401.501.501.601.601.701.701.701.701.702.102.302.302.302.402.401.801.802.202.202.202.202.302.302.302.662.662.6

82.682.702.802.803.003.103.103.103.103.204.304.301.501.801.802.002.302.302.402.402.402.402.502.502.602.602.602.702.703.103.10;run;SAS过程步procnpar1wayw

ilcoxondata=aa;varx;classgroup;run;procrankdata=aaout=bb;varx;ranksr;run;procglmdata=bb;classgroup;modelr=group;meansgroup/dunnett("1");run;p

rocfreqdata=aa;tablesgroup*x/scores=rankcmh;run;Kruskal-Wallis检验结果：R-Square×Df总行均值得分的差值Rowmeanscoresdiffer22

练习2•正常人、单纯型肥胖和皮质醇增多症病人的血浆总皮质醇测定值如下（µg/L）：–正常人：1.41.92.22.52.83.13.73.94.67.0–单纯型肥胖：1.61.22.02.43.14.15.

07.410.813.6–皮质醇增多症：9.810.210.613.014.014.815.615.621.624.0•试用非参数的方法比较三组人群的血浆总皮质醇含量有无差别？23（四）、随机区组设计资料的秩和检验实例•用某新药治疗血吸虫病患者，

采用三天疗法，在治疗前及治疗后测定7名患者的血清谷丙转氨酶（SGPT）的变化，以观察该药对肝功能的影响，测定结果如下。•试用非参数检验的方法比较四个阶段SGPT有无差别？24SAS数据集的建立dataaa;dobloc

k=1to7;dotreat=1to4;inputx@@;output;end;end;cards;63188138549023822014454300839245140213100541751503672300163906420718587;run;SAS过程步pro

csortdata=aa;byblock;run;procrankdata=aaout=b;varx;byblock;ranksr;run;procglmdata=b;classtreatblock;modelr=treatblock;run;procfreqdata=aa;ta

blesblock*treat*x/scores=rankcmh;run;行均值得分的差值RowmeanscoresdifferFriedman检验结果：R-Square×（Df总－Df区组）25小结（1）

配对设计和单样本资料的非参数检验Wilcoxon符号检验和Wilcoxon符号秩和检验SAS中Univariate过程：“位置检验:Mu0=0”中的•符号M•符号秩S完全随机设计两组独立样本的秩和检验Wi

lcoxon秩和检验SAS中•NPAR1WAY过程（wilcoxon）•FREQ过程（tables…./scores＝rankcmh）“行均值得分的差值”26小结（1）完全随机设计多组独立样本的秩和检验Kruskal-Wall

is秩和检验SAS中•NPAR1WAY过程•Rank过程+Glm过程（或ANOVA过程）“R-Square(R2)×Df总”结果•Freq过程（tablesgroup*x/scores＝rankcmh）“行均值得分的差值”随机区组设计资料的秩和检验Friedman秩和检验

SAS中•Freq过程（tablesblock*group*x/scores＝rankcmh）“行均值得分差值”•Sort过程（block）+Rank过程+Glm过程“R-Square×（Df总－Df区组）”结果27小结（2）符合正态分布、方差齐性条件参数检验•t检验•

方差分析（单向、双向等）不符合正态分布、方差齐性条件数据变换使其满足条件，方差分析数据变换仍不满足条件，非参数检验•wilcoxon秩和检验•Kruskal-Wallis秩和检验•Friedman秩和检验28