wilcoxon符号秩检验学习.pptxVIP

第1页/共33页例 2. 4 下面是10个欧洲城镇每人每年平均消费的酒量（相当于纯酒精数）（单位：升）。数据已经按升幂排列。 4.12 5.18 7.63 9.74 10.39 人们普遍认为欧洲各国人均年消费酒量的中 位数相当于纯酒精8升，也就是me0=8。由数据 算得的中位数为。因此，我们的检验设为： H0：me＝8，H1：me 8第2页/共33页 先计算每个样本值和原假设中me0的值之差，即Xi－8。 考虑这些差的绝对值并将绝对值从小到大排序，从而求出这些绝对值的秩。 再计算比8大的样本对应的绝对值的秩之和，如果这个和比较大，我们就拒绝原假设，接受备择假设。第3页/共33页问题一般提法： 假定样本X1, … , X n来自分布连续对称的总体X，在此假定下总体X的中位数等于均值。 问题主要是检验中位数，即原检验为H0：me=me0，相对于各种单双边的备择假设。第4页/共33页注： （1）与符号检验不同： Wilcoxon符号秩检验假设总体分布是对称的。（2）在总体分布对称的假设下，即设总体X的分布关于点θ对称，则X的均值和中位数相同，且均为θ。所以检验总体中位数可等价于检验总体对称中心。即检验的原假设 H0：M=M0 等价于 H0：θ=θ0（相对于各种单双边的备择假设）。第5页/共33页检验步骤： H0： θ＝θ0 （对应于各单双边备择假设） Step 1. 计算i=1, 2, … , n。记差为z i. Step 2. 将差z i.的绝对值，即… ,按从小到大的顺序排列。由于总体服从连续型分布，不妨假定样本互不相等，都不等于0，且样本差的绝对值也互不相等。所以可得到样本z i.的绝对值的秩，不妨记 的秩为R i。 第6页/共33页Step 3. 符号秩和检验统计量为 其中 或者取检验统计量为 其中主要取W＋为检验统计量。第7页/共33页Step 4 设w＋表示由样本算出的W＋的值。 （1） H0： θ＝θ0 , H1： θθ0 p值＝P( W＋≥ w＋ )； （2） H0： θ＝θ0 , H1： θθ0 p值＝P( W＋≤ w＋ )； （3） H0： θ＝θ0 , H1： θ≠θ0 p值＝2min{P( W＋≥ w＋ )，P(W＋≤ w＋)}第8页/共33页 对Step 4的注解： 对于对称中心不为0的总体分布，可以转化为中心为0的情况进行检验！ 现不妨假设θ0＝0，则原假设变为 H0：θ＝0 对于这种检验，通过严格的证明来说明p值的选取。第9页/共33页 （1）H0： θ＝0 , H1： θ0。 若H1成立，则总体X的分布关于点θ对称。从而有， P( X0 ) P( X0 ) ，且对任意正数a， P( Xa ) P( X－a )。 所以当H1成立，不仅观察到的取正值的样本数据的个数比较多，且取正值的样本数据的拒绝值也比较大。由此，H1成立时，W＋的值较大 。所以p值＝P( W＋≥ w＋)。 第10页/共33页例 2. 2中我们的检验设为： H0：M＝8，H1：M 8下面来用Wilcoxon符号秩检验，等价于检验 H0：θ＝8，H1： θ 8第11页/共33页检验步骤 Step 1. 对于 i=1, 2, … , n，计算得到新的样本zi和它们对应的秩如下：样本 xi4.125.187.639.7410.3911.9212.3212.913.5414.45 zi的符号－－－ ＋ ＋ ＋ ＋ ＋ ＋ ＋zi的绝对值3.882.190.371.742.393.924.324.895.546.45 秩 5 3 1 2 4 6 7 8 9 10第12页/共33页 Step 2. 计算W＋。 W+=2+4+6+7+8+9＋10＝46 利用W＋的分布，辅以统计软件，可计算出p值＝ 。Step 3. 所以给定α＝时，此时可拒绝原假设，认为欧洲人均酒精年消费多于8升。 第13页/共33页W＋的分布性质 设独立同分布样本x1,…,xn来自连续对称总体X,X分布的对称中心为θ。为方便讨论，不妨设原假设为 H0：θ＝0，即总体分布关于原点0对称的条件下，讨论W＋的性质。 注：W＋与W－有下列关系： W++ W- = n(n＋1)/2 第14页/共33页（关键）性质 令 则在总体的分布关于原点0对称时，W＋与S同分布。注： S是W＋当Ri＝i时的特殊情况。研究W＋的分布可转为研究S的分布。第15页/共33页概率分布性质 在总体的分布关于原点0对称时，W＋的概率分布为 P ( W+ = d )=P ( S＝d ) =t n(d)/2n,其中，d＝0, 1, 2, … , n(n+1)/2，tn (d)表示从1, 2, … , n这n个数中任取若干