SPSS秩和检验--实例分析

SPSS统计——秩和检验第1页/共18页SPSS秩和检验应用条件①总体分布形式未知或分布类型不明；②偏态分布的资料：③等级资料：不能精确测定，只能以严重程度、优劣等级、次序先后等表示；④不满足参数检验条件的资料：各组方差明显不齐。⑤数据的一端或两端是不确定数值，如“50mg”等。一、配对资料的Wilcoxon符号秩和检验（Wilcoxonsigned-ranktest）例1对10名健康人分别用离子交换法与蒸馏法，测得尿汞值，如表9.1的第（2）、（3）栏，问两种方法的结果有无差别？表110名健康人用离子交换法与蒸馏法测定尿汞值(μg/l)样品号（1）离子交换法（2）蒸馏法（3）差值(4)=(2)(3)秩次（5）10.50.00.5222.21.11.1730.00.00.0—42.31.31.0656.23.42.8861.04.6-3.6-971.81.10.73.584.44.6-0.2-192.73.4-0.7-3.5101.32.1-0.8-5T+=+26.5T-=-18.5差值先进行正态性及方差齐性检验，看是否可以做参数检验，其检验效能高于非参数检验。（下同）H0：Md（差值的总体中位数）=0H1：Md≠0α=0.05T++T-=1+2+3+…n=n(n+1)/2①小样本（n≤50）--查T界值表基本思想：如果无效假设H0成立，则正负秩和的绝对值从理论上说应相等，都等于n(n+1)/4，既使有抽样误差的影响正负T值的绝对值相差也不应过大。反过来说，如果实际计算出的正负T值绝对值相差很大，我们只能认为H0成立的可能性很小。界值的判断标准若下限T上限，P值表中概率值若T≤下限或T≥上限，则P值≤表中概率值②大样本时（n50），正态近似法（u检验）基本思想：假定无效假设H0成立，则正负秩和的绝对值应相等，随着n增大T逐渐趋近于均数等于n(n+1)/4、方差为n(n+1)(2n+1)/24的正态分布。所以可用近似正态法计算u值。即：4824121413jjttnnnnnTuSPSS统计——秩和检验第2页/共18页*校正公式：当相同秩次个数较多时48)(24)12)(1(5.0|4/)1(|3jjttnnnnnTutj：第j个相同秩次的个数SPSS:建立变量名：录入数值：统计分析：SPSS统计——秩和检验第3页/共18页分析——非参数检验——两相关样本（配对样本）结果分析：SPSS统计——秩和检验第4页/共18页表一：第一行：b-a的负秩（NegativeRanks）有5个（右上角的a在表下方有注释），平均秩次为5.3，负秩和为26.5。第二行：正秩，正秩的个数，平均秩次，正秩和。表二：Z即为u值，可用正秩和18.5或负秩和26.5计算，习惯上用较小的秩和计算u值。p=0.635大于0.05，不拒绝H0，还不能认为两种方法有差别。二、两个独立样本比较的Wilcoxon秩和检验(Wilcoxonranksumtest)1．原始数据的两样本比较例2某实验室观察局部温热治疗小鼠移植肿瘤的疗效，以生存日数作为观察指标，试检验两组小鼠生存日数有无差别？实验组对照组生存日数秩次生存日数秩次109.5211212.5321515431516541617651718761819872020982321109.590以上221111SPSS统计——秩和检验第5页/共18页1212.51314n1=10T1=170n2=12T2=83时间资料不服从正态分布H0：两总体分布位置相同H1：两总体分布位置不同a=0.05记n较小组秩和为T，样本量n1。如果n1=n2，可取任秩和①查表法：查T界值表：n1≤10，n2n1≤10界值的判断标准：若下限T上限，P值表中概率值若T≤下限或T≥上限，则P值≤表中概率值②正态近似法当n1或n2-n1超出T界值表的范围时，随n增大，T的分布逐渐逼近均数为n(n+1)/4、方差为n(n+1)(2n+1)/24的正态分布，所以可用近似正态法计算u值。即：12/)1(5.0|2/)1(|2111NnnNnTu*校正公式（当相同秩次较多时）cuuC/)-/()-(-133NNttCjjSPSS建立变量名：录入数值：SPSS统计——秩和检验第6页/共18页统计分析：SPSS统计——秩和检验第7页/共18页结果分析：Z值为-3.630，p＜0.001，拒绝H02.频数表资料（或等级资料）的两样本比较例320名正常人和32名铅作业工人尿棕色素定性检查结果见下表。问铅作业工人尿棕色素是否高于正常人？人数秩和结果（1）正常人(2)铅作业工人(3)合计(4)秩次范围（5）平均秩次（6）正常人(7)=(2)(6)铅作业工人(8)=(3)(6)SPSS统计——秩和检验第8页/共18页-188261-2613.52431082101227-3832.565325++07739-4542.00294+++03346-4847.00141++++04449-5250.50202合计n1=20n2=3252--T1=308T2=1070取n较小组的秩和为T值，用校正公式计算。即：12/)1(5.0|2/)1(|2111NnnNnTucuuC/)-/()-(-133NNttCjjSPSS：建立变量名：录入数值：统计分析：SPSS统计——秩和检验第9页/共18页SPSS统计——秩和检验第10页/共18页结果分析：同两个独立样本比较的Wilcoxon秩和检验P＜0.001，拒绝H0三、多个样本比较的秩和检验(Kruskal-WallisHtest)1．原始数据法例4某研究者测定正常人、单纯性肥胖、皮质醇增多症者各10人的血浆总皮质SPSS统计——秩和检验第11页/共18页醇含量见下表，问这三组人的血浆总皮质醇含量有无差别？三组人的血浆总皮质醇含量测定值（μg/L）正常人单纯性肥胖皮质醇增多症测定值秩次测定值秩次测定值秩次0.410.629.8201.941.2310.2212.262.0510.6222.582.4713.0232.893.110.514.0253.110.54.11414.8263.7125.01615.6273.9135.91715.6284.6157.41921.6296.01813.62424.030Rini96.510117.51025110H0:：三组人的血浆总皮质醇含量总体分布位置相同H1：三组人的血浆总皮质醇含量总体分布位置不全相同a=0.05)1(3)1(122NnRNNHii1kSPSS建立变量名录入数值：SPSS统计——秩和检验第12页/共18页统计分析：SPSS统计——秩和检验第13页/共18页结果分析：SPSS统计——秩和检验第14页/共18页若g（组数）=3且最小样本列数大于5或g3时，H或HC近似服从自由度为g-1的卡方分布。H=18.130，自由度=2，P＜0.001，拒绝H0，三组总体分布位置不全相同，需做两两比较。2．频数表法：例5四种疾病患者痰液内嗜酸性粒细胞的检查结果见表6-5。问四种疾病患者痰液内嗜酸性粒细胞有无差别？四种疾病患者痰液内嗜酸性粒细胞比较白细胞（1）支气管扩张（2）肺水肿（3）肺癌（4）合计（6）秩次范围（7）平均秩次（8）-035111—116+2571912—3021++9532031—5040.5+++6221051—6055.5Ri739.5436.5409.5———ni17.0151760——iR43.5029.1024.09———)1(3)1(122NnRNNHiicHHcSPSS统计——秩和检验第15页/共18页)()(133NNttCjj1kSPSS建立变量名：录入数值：统计分析：SPSS统计——秩和检验第16页/共18页SPSS统计——秩和检验第17页/共18页结果分析：同多样本比较的秩和检验P＜0.001，拒绝H0，三组总体分布位置不全相同，需做两两比较。秩和检验方法要点和注意事项检验方法方法要点注意事项配对样本的符号秩检验1.依差值大小编秩，再冠以差值的符号，任取T+、T-作为T，1.编秩时若差值绝对值相同符号相反，取平均秩次。0差值SPSS统计——秩和检验第18页/共18页查附表9，T界值表。TT界值，Pα。2.n50，用z检验。省略。2.n5，不能得有意义结论。两独立样本的秩和检验（分布位置）1.按两组数据由小到大统一编秩，以n1较小者为T，查附表10T界值表。T在界值范围内，Pα。2.n110或n1-n210时，用z检验。1.编秩时若相同数据在不同组，取平均秩次。2.当相同秩次较多时，使用校正公式。成组设计多样本比较的秩和检验(K-W检验)1.将k组数据由小到大统一编秩，求各组秩和Ri。2.计算H值，用ν=k-1查χ2界值表，确定P值。3.拒绝H0时，应作多个样本两两比较的秩和检验。1.编秩时若相同数据在不同组，取平均秩次。2.当相同秩次较多时，使用校正公式。