第八章、非参数检验课件.pptVIP

随机区组设计的秩和检验* 上表中rij表示第i处理组(共k个处理组)第j个区组(共b个区组)的观察值对应的相应秩次 任意区组的秩和为： 区组号 处理组1 处理组2 …… 处理组k 区组1 r11 r21 …… rk1 区组2 r12 r22 …… rk2 …… …… …… …… …… 区组b r1b r2b …… rkb 随机区组设计的秩和检验* 则k个处理组，b个区组所有观测对象的秩次之和为： 如果H0成立，则认为各处理组的秩次之和应该差不多，理论上都等于T/k(每组的例数均为b例，共有k组，每组的秩和Ti等于b(k+1)/2) 故任意组的秩和的期望值与方差为： 随机区组设计的秩和检验* 但是与前面H统计量不同的是，所有各组的秩和经过标准化后求得的平方和并不满足卡方分布 随机区组设计的秩和检验* 六、多个样本间的多重比较* 无论是对完全随机设计多个样本比较用Kruskal-Wallis秩和检验，还是对随机化区组设计用Friedman秩和检验，当推断结论为拒绝H0，接受H1时，与方差分析类似，只能得出各总体分布不同或不全相同的结论，但不能说明任两个总体分布不同 若要对每两个总体分布做出有无不同的推断，需要作组间的多重比较 秩和检验注意事项 秩和检验功效问题： 若满足t检验的前提条件，在样本含量较小时t检验和秩和检验的功效均很低，t检验的功效略高一些；而当样本含量增加时，两种检验的功效随之增加，且差别不大 若样本来自非正态分布，在样本含量小时， t检验的功效很低，Wilcoxon秩和检验的功效较高；随着样本含量的增大，Wilcoxon秩和检验的功效更高于t检验的功效，有时甚至高出很多 如果不清楚资料是否符合参数检验的检验条件，不要贸然应用参数检验方法 阶段总结 检验方法的选择 ： 一组样本资料：若来自正态总体，可用t检验；若来自非正态总体或总体分布无法确定，可用Wilcoxon符号秩和检验方法 配对设计资料：二分类变量，可用McNemar检验；连续型变量，若来自正态总体，可用配对t检验，否则可用Wilcoxon符号秩和检验 两组独立样本：连续型变量，若来自正态总体可用t检验，否则可用Wilcoxon秩和检验；二分类变量，可用卡方检验；无序多分类变量，可用卡方检验；有序多分类变量，宜用Wilcoxon秩和检验 多组独立样本：连续型变量值，来自正态总体且方差相等，可用方差分析；否则，进行数据变换使其满足正态性或方差齐的要求后，采用方差分析；数据变换仍不能满足条件时，可用Kruskal-Wallis秩和检验；二分类变量或无序多分类变量，可用卡方检验；有序多分类变量宜用Kruskal-Wallis秩和检验。 随机区组设计：连续型变量，来自正态总体且方差相等，可用随机区组设计方差分析；否则进行数据变换使其满足正态性或方差齐的要求后，用方差分析进行分析；数据变换仍不能满足条件时，可用Friedman秩和检验 秩和分布的正态近似性 若秩次相同的现象出现较多（如超过不等0的总对子数n的25%），用上式求得的u值偏小，应计算校正的统计量值uc 式中tj为第j(j=1,2…)次相持时，所含相同秩次的个数 ；在本例中，共有2例秩次相持，故S(tj3-tj)= (23-2)=6 单样本资料的符号秩和检验 若单组随机样本来自正态总体，比较其总体均数与某常数是否不同，可用t检验 若样本来自非正态总体或总体分布无法确定，也可用Wilcoxon符号秩和检验，检验总体中位数是否等于某已知数值 所不同的只是差值为各观察值与已知总体中位数之差，其他符号的意义同配对设计资料 单样本资料的符号秩和检验 已知某地正常人的尿氟含量的中位数为45.30μmol/L；在该地的某厂随机抽取12名工人，测得尿氟含量见下表，问该厂工人的尿氟含量是否高于正常人？ 尿氟含量 差值d 秩次 尿氟含量 差值d 秩次 44.21 -1.09 53.26 7.96 45.30 0 — 54.37 9.07 46.39 1.09 57.16 11.86 49.47 4.17 67.37 22.07 51.05 5.75 71.05 25.75 53.16 7.86 87.37 42.07 -1.5 1.5 3 4 5 6 7 8 9 10 11 单样本资料的符号秩和检验 计算步骤如下： 求各数据与已知总体中位数的差值d H0：差值的总体中位数等于零，即：Md=0 H1：差值的总体中位数不等于零，即Md≠0 a=0.05(单侧) 按差值的绝对值由小到大编秩，并按差值的正负给秩次加上正负号；编秩时，若差值为0，舍去不计；若差值的绝对值相等，称为相持（tie），这时取平均秩次 合计T+=6