1、spss非参数检验K多个独立样本检验Kspss-非参数检验-K多个独立样本检验（ Kruskal-Wallis检验）案例解析 最近经常失眠，好痛苦啊！ 大家有什么好的解决失眠的方法吗？希望知道的能够告诉我，谢谢啦，今天和大家一起探讨和分下一下SPSS-非参数检验-K个独立样本检验 （ Kruskal-Wallis检验）。 还是以SPSS教程为例：假设：HO: 不同地区的儿童，身高分布是相同的 H1： 不同地区的儿童，身高分布是不同的不同地区儿童身高样本数据如下所示：提示：此样本数为4个（北京，上海，成都 ，广州）每个样本的样本量（观察数）都为5个即：K=43 n=5, 此时如果样本逐渐增大，呈
2、现出自由度为K-1的平方的分布，（即指：卡方检验）点击“分析”非参数检验旧对话框K个独立样本检验，进入如下界面：将“周岁儿童身高”变量拖入右侧“检验变量列表”内， 将“城市（CS)变量” 拖入“分组变量”内，点击“定义范围” 输入“最小值”和“最大值”（这里的变量类型必须为“数字型”）如果不是数字型，必须要先定义或者重新编码。在“检验类型”下面选择“秩和检验”（ Kruskal-Wallis检验）点击确定运行结果如下所示：对结果进行分析如下：1：从“检验统计量a,b”表中可以看出：秩和统计量为：13.900 自由度为：3=k-1=4-1下面来看看“秩和统计量”的计算过程，如下所示：假设“秩和统
3、计量”为 kw 那么:其中：n+1/2 为全体样本的“秩平均” Ri./ni 为第i个样本的秩平均 Ri.代表第i个样本的秩和， ni代表第i个样本的观察数）最后得到的公式为：北京地区的“秩和”为： 秩平均*观察数（N) = 14.4*5=72上海地区的“秩和”为：8.2*5=41成都地区的“秩和”为：15.8*5=79广州地区的“秩和”为：3.6*5=18接近13.90 （由于中间的计算，我采用四舍五入，丢弃了部分数值，所以，会有部分误差）2：“检验统计量a,b”表中可以看出：“渐进显著性为0.003， 由于0.0030.01 所以得出结论：H1： 不同地区的儿童，身高分布是不同的秩和检验
4、前面介绍的均数的区间估计及假设检验，都是要求个体变量值服从正态分布，或根据中心极限定理，当样本较大时，样本均数服从正态分布。这种要求样本来自总体分布型是已知的，在此基础上对总体参数进行估计或检验，称为参数统计（parametric statistics）。但在医学研究中，许多数据不符合参数统计的要求，这时有两种处理的方法。一是，进行数据转换，使其符合参数统计方法的要求。二是，选择非参数检验方法，非参数检验(non-parametric test)方法是对样本来自的总体分布不作要求（如不要求样本来自正态分布）的一类假设检验方法。非参数检验的主要优点是对样本的总体分布不作要求，适用的范围广，尤其是
5、当变量中有不确定数值时，如50mg”等。 一、配对资料的Wilcoxon符号秩和检验（Wilcoxon signed-rank test） 例1 对10名健康人分别用离子交换法与蒸馏法，测得尿汞值，如表9.1的第（2）、（3）栏，问两种方法的结果有无差别？表1 10名健康人用离子交换法与蒸馏法测定尿汞值(g/l)样品号（1）离子交换法（2）蒸馏法（3）差值(4)=(2) (3)秩次（5）10.50.0 0.5 222.21.1 1.1 730.00.0 0.042.31.3 1.0 656.23.4 2.8 861.04.6-3.6-971.81.1 0.7 3.584.44.6-0.2-19
6、2.73.4-0.7 -3.5101.32.1-0.8-5T+=+26.5 T-=-18.5差值先进行正态性及方差齐性检验，看是否可以做参数检验，其检验效能高于非参数检验。（下同）H0：Md（差值的总体中位数）=0 H1：Md0 =0.05T+T-=1+2+3+n=n(n+1)/2 1 小样本（n50）-查T界值表基本思想：如果无效假设H0成立，则正负秩和的绝对值从理论上说应相等，都等于n(n+1)/4，既使有抽样误差的影响正负T值的绝对值相差也不应过大。反过来说，如果实际计算出的正负T值绝对值相差很大，我们只能认为H0成立的可能性很小。界值的判断标准 若下限T表中概率值 若T下限或T上限，则
7、P值表中概率值 2 大样本时（n50），正态近似法（Z检验）基本思想：假定无效假设H0成立，则正负秩和的绝对值应相等，随着n增大T逐渐趋近于均数等于n(n+1)/4、方差为n(n+1)(2n+1)/24的正态分布。所以可用近似正态法计算Z值。即：*校正公式：当相同秩次个数较多时 tj：第j个相同秩次的个数SPSS:建立变量名：录入数值：统计分析：分析非参数检验两相关样本（配对样本）结果分析：表一：第一行：b-a的负秩（Negative Ranks）有5个（右上角的a在表下方有注释），平均秩次为5.3，负秩和为26.5。第二行：正秩，正秩的个数，平均秩次，正秩和。表二：可用正秩和18.5或负秩和
8、26.5计算，习惯上用较小的秩和计算Z值。p=0.635大于0.05，不拒绝H0，还不能认为两种方法有差别。二、两个独立样本比较的Wilcoxon秩和检验(Wilcoxon rank sum test)1原始数据的两样本比较例2 某实验室观察局部温热治疗小鼠移植肿瘤的疗效，以生存日数作为观察指标，试检验两组小鼠生存日数有无差别？ 实验组对照组生存日数秩次生存日数秩次10 9.5 2 112 12.5 3 21515 4 31516 5 41617 6 51718 7 61819 8 72020 9 8232110 9.590以上2211111212.51314n1=10T1=170n2=12T
9、2=83时间资料不服从正态分布H0：两总体分布位置相同 H1：两总体分布位置不同 a=0.05记n较小组秩和为T，样本量n1。如果n1=n2，可取任秩和1 查表法：查T界值表：n110，n2n110界值的判断标准： 若下限T表中概率值 若T下限或T上限，则P值表中概率值2 正态近似法当n1或n2-n1超出T界值表的范围时，随n增大，T的分布逐渐逼近均数为n(n+1)/4、方差为n(n+1)(2n+1)/24的正态分布，所以可用近似正态法计算Z值。即： *校正公式（当相同秩次较多时）SPSS建立变量名：录入数值：统计分析：结果分析：Z值为-3.630，p0.001，拒绝H02. 频数表资料（或等
10、级资料）的两样本比较 例3 20名正常人和32名铅作业工人尿棕色素定性检查结果见下表。问铅作业工人尿棕色素是否高于正常人？ 结果（1）人数秩次范围（5）平均秩次（6）秩和正常人(2)铅作业工人(3)合计(4)正常人(7)=(2)(6)铅作业工人(8)=(3)(6)-18 8261-2613.5243108 2101227-3832.5 65325+0 7 739-4542.00294+0 3 346-4847.00141+0 4 449-5250.50202合计n1=20n2=3252-T1=308T2=1070取n较小组的秩和为T值，用校正公式计算。即： SPSS：建立变量名：录入数值：统计
11、分析：结果分析：同两个独立样本比较的Wilcoxon秩和检验Mann-Whitney TestP0.001，拒绝H0三、多个样本比较的秩和检验(Kruskal-Wallis H test)1原始数据法例4 某研究者测定正常人、单纯性肥胖、皮质醇增多症者各10人的血浆总皮质醇含量见下表，问这三组人的血浆总皮质醇含量有无差别？ 三组人的血浆总皮质醇含量测定值（g/L）正常人单纯性肥胖皮质醇增多症测定值秩次测定值秩次测定值秩次0.41 0.62 9.8201.94 1.23 10.221 2.26 2.05 10.6222.58 2.47 13.0232.89 3.1 10.5 14.0253.1
12、10.5 4.1 14 14.8263.7 12 5.0 16 15.6273.9 13 5.9 17 15.6284.6 15 7.4 19 21.6296.0 18 13.6 24 24.030Ri ni 96.5 10 117.5 10 251 10H0:：三组人的血浆总皮质醇含量总体分布位置相同H1：三组人的血浆总皮质醇含量总体分布位置不全相同a=0.05 SPSS建立变量名录入数值：统计分析：结果分析：若g（组数）=3且最小样本例数大于5或g3时，H或HC近似服从自由度为g-1的卡方分布。H=18.130，自由度=2， P0.001，拒绝H0，三组总体分布位置不全相同，需做两两比较。
13、2频数表法：例5 （单向有序分类变量的多个样本比较）用A、B、C三种药物治疗单纯性慢性支气管炎，结果见表第（1）（5）栏, 问三种药物的总体疗效是否不同。表 三种药物疗效比较的秩和检验计算过程疗效（1）药物合计（5）秩次范围（6）平均秩次（7）秩和R1(8)=(2) (7)秩和R2(9)=(3) (7)秩和R3(10)=(4) (7)A（2）B（3）C（4）治愈175325 12513.0221.065.039.0显效511117792610465.03315.0715.01105.0好转335247132105236170.55626.58866.08013.5无效7242657237293
14、265.01855.06360.06890.0合计1089293293 11017.516006.016047.5检验步骤如下：(1) 建立检验假设H0: 三种药物疗效的总体分布相同H1: 三种药物疗效的总体分布不同或不全相同 =0.05(2) 编秩用各疗效等级的合计值排序确定秩次范围，如表第（6）栏所示，A、B、C三种药物总的治愈人数是25，他们的秩次范围是125。同理疗效为“显效”组的秩次范围是26104，以此类推。再对第（6）栏秩次范围的上下限求和取平均值得各等级的平均秩次，如第（7）栏所示，疗效为“治愈”组的平均秩次为=13。(3) 求秩和分别用第（2）（4）栏各等级的频数与（7）栏平
15、均秩次相乘再求和，如第（8）（10）栏所示。(4) 计算统计量H值将第（8）（10）栏的总秩和T1、T2、T3代入公式(11.15)计算H值。若各样本相同秩次较多时（如超过25%），由公式（11.15)计算所得H值偏小，应按公式（11.16）和公式（11.17）对H值作校正计算HcH =() -3(N+1) H = (+) -3(293+1) = 48.23C =1- (tj3-tj) / (N3-N) C = 1- (253-25) (793-79) (1323-132) (573-57) /(2933-293) = 0.9086 Hc=H/C Hc = 48.23/0.9086=53.08
16、 (5) 确定P值和作出统计结论若组数k=3, 每组例数均小于或等于5，可查附表H界值表，得出P值。本例各组例数均大于5，已超出附表的范围，则H值近似服从=k-1的2分布，可查附表的2界值表。故按=3-1=2，查2界值表得20。05（2）=5.99，因为53.085.99，故PT界值，P。2. n50，用z检验。编秩时若差值绝对值相同符号相反，取平均秩次。0差值省略。两独立样本的秩和检验（分布位置）1. 按两组数据由小到大统一编秩，以n1较小者为T，查附表10 T界值表。T在界值范围内，P。2. n110或n1- n210时，用z检验。1. 编秩时若相同数据在不同组，取平均秩次。2. 当相同秩次较多时，使用校正公式。成组设计多样本比较的秩和检验(K-W检验)1. 将k组数据由小到大统一编秩，求各组秩和Ri。2. 计算H值，用=k-1查2界值表，确定P值。3. 拒绝H0时，应作多个样本两两比较的秩和检验。1. 编秩时若相同数据在不同组，取平均秩次。2. 当相同秩次较多时，使用校正公式。
