SAS讲义_第二十七课符号检验和Wilcoxon符号秩检验

合集下载

第二十七课符号检验和Wilcoxon符号秩检验

第二十七课符号检验和Wilcoxon符号秩检验

0.0898 0.2120 0.3953 0.6047 0.7880 0.9102 0.9713 0.9935 0.9991 0.9999 1.0000
从表 27.2 的累计概率列中我们看到, 正号出现的次数大于 10 的概率为 1-0.9713=0.0287, 或者换一种方法计算为=0.0001+0.0009+0.0056+0.0222=0.0287, 二者的微小差异是因为小数点 后舍入问题造成的。而试验的结果:正号出现的次数为 11,大于 10,出现的概率不会超过 0.0287,我们开始设定的显著性水平为 0.1,由于 0.0287<0.1,所以我们拒绝原假设,接受备 选假设。如果我们的原假设为 p =0.5,既训练前后学生素质相等,那么就是双侧检验,应该 加上正号出现的次数小于 4 的概率 0.0287,即 2×0.0287=0.0574<0.1,同样是拒绝原假设,接 受区间为 4 次到 10 次,而拒绝区间为小于等于 3 次(小于 4 次)或大于等于 11 次(大于 10 次) 。 2. 大样本时的正态近似概率计算 当 n 20 时,样本可以认为是大样本。我们可以利用二项分布的正态近似,即对于
由于试验的结果只有两种可能,正号或负号,对每一个学生试验出现正号的假定概率为 p =0.5,负号为 1— p =0.5,这样整个试验的概率是相同的,并且每一个试验是相互独立的。 因此在 n =14 次独立的试验中,正号出现的次数服从二项分布 B(14,0.5) ,见表 27.2 所示。 表 27.2 二项分布的概率和累计概率 n=14,p=0.5 正号出现的次数 0 1 2 3
上海财经大学经济信息管理系IS/SHUFE Page 4 of 8

S T

Wilcoxon符号秩检验

Wilcoxon符号秩检验

04 Wilcoxon符号秩检验 的优缺点
优点
无需假设数据分布
Wilcoxon符号秩检验是一种非参 数检验方法,不需要假设数据服 从特定的分布,因此对于不符合 正态分布的数据也能得到较为准 确的结果。
对异常值不敏感
由于Wilcoxon符号秩检验是基于 秩次的检验方法,因此对于异常 值的存在并不敏感,能够得到较 为稳健的结果。
适用于配对样本
Wilcoxon符号秩检验适用于配对 样本的比较,能够充分利用样本 信息,提高检验的效能。
缺点
检验效能较低
相比于参数检验方法,如t检验,Wilcoxon符号秩检验的检 验效能较低,即当存在真实的差异时,该方法可能无法准 确地检测出差异。
对样本量要求较高
为了得到较为准确的检验结果,Wilcoxon符号秩检验需要 较大的样本量。当样本量较小时,该方法的准确性可能会 受到影响。
正态分布,因此适用于更广泛的数据
类型。
符号秩检验的定义
符号
在Wilcoxon符号秩检验中,首先计算每对观测值 之间的差值,并根据差值的正负赋予相应的符号 (+或-)。
检验统计量
根据符号和秩次计算检验统计量,通常使用 Wilcoxon符号秩统计量(W)或标准化后的z统 计量。这些统计量用于衡量两组观测值之间的差 异显著性。
非参数统计方法
Wilcoxon符号秩检验是一种非参数统 计方法,用于比较两个相关样本、配 对观测值或重复测量之间的差异。
稳健性
由于不对数据分布做严格假设, Wilcoxon符号秩检验对于异常值和偏 离正态分布的数据具有较好的稳健性 。
无需正态分布假设
与参数检验(如t检验)不同,
Wilcoxon符号秩检验不需要数据服从

SAS讲义_第二十七课符号检验和Wilcoxon符号秩检验

SAS讲义_第二十七课符号检验和Wilcoxon符号秩检验

SAS讲义_第二十七课符号检验和Wilcoxon符号秩检验第二十七课符号检验和Wilcoxon 符号秩检验在统计推断和假设检验中,传统的检验统计量都叫做参数检验,因为它们都依赖于确定的概率分布,这个分布带有一组自由的参数。

参数检验被认为是依赖于分布假定的。

通常情况下,我们对数据进行分析时,总是假定误差项服从正态分布,这是人们易于接受的事实,因为正态分布的原始出发点就是来自于误差分布,至于当样本相当大时,数据的正态近似,这是由于大样本理论所保证的。

但有些资料不一定满足上述要求,或不能测量具体数值,其观察结果往往只有程度上的区别,如颜色的深浅、反应的强弱等,此时就不适用参数检验的方法,而只能用非参数统计方法(non-parametric statistical analysis )来处理。

这种方法对数据来自的总体不作任何假设或仅作极少的假设,因此在实用中颇有价值,适用面很广。

一、单样本的符号检验符号检验(sign test )是一种最简单的非参数检验方法。

它是根据正、负号的个数来假设检验。

首先需要将原始观察值按设定的规则,转换成正、负号,然后计数正、负号的个数作出检验。

该检验可用于样本中位数和总体中位数的比较,数据的升降趋势的检验,特别适用于总体分布不服从正态分布或分布不明的配对资料,有时当配对比较的结果只能定性的表示,如试验前后比较结果为颜色从深变浅、程度从强变弱,成绩从一般变优秀,即不能获得具体数字,也可用符号检验,例如用正号表示颜色从深变浅,用负号表示颜色从浅变深。

用于配对资料时,符号检验的计算步骤为:首先定义成对数据指定正号或负号的规则,然后计数正号的个数+S 及负号的个数-S ,由于在具体比较配对资料时,可能存在配对资料的前后没有变化,或等于假设中的中位数,此时仅需要将这些观察值从资料中剔除,当然样本大小n 也随之减少,故修正样本大小-++=S S n 。

当样本n 较小时,应使用二项分布确切概率计算法,当样本n 较大时,常利用二项分布的正态近似。

符号检验和符秩检验

符号检验和符秩检验
H0: =0.5。
计 • 建立原假设为真前提下的下列检验统计
量:

Z P 0.5 ~ N (0,1)
0.25
n
天津财经大学 统计学系
三、威尔科克森配对符号秩检验
统 • 以上所介绍的两总体情形下符号检验方法,
用配对观测之间差别的符号进行检验,而不

注重差别的大小,因此对资料的利用不够充 分。

天津财经大学 统计学系
• 如果改革没有引起居民经济情况的变 化,那么居民经济情况的前后差异就

完全是由于各种随机因素的影响形成
的(假定其它重要的影响因素都已控

制不变),于是正差值的个数与负差
值的个数会大体相等。
学 • 把0差值舍去后,对总体(正差值与负 差值组成的总体)作独立重复贝努里
试验,每次试验出现正号的概率是
统 • 取正差值的个数,近似服从正态分布, 所以通常可以将其标准化为标准正态变
计 量,作为检验统计量。 • 建立原假设为真前提下的下列检验统计
学 量:
Z 0.5n ~ N (0,1)
0.25n
天津财经大学 统计学系
二、两总体问题的符号检验
• 两总体符号检验适用于检验配对样本情

形下,两总体分布在位置特征上是否有

符号检验与符秩检验


主讲人:杨贵军 教授
天津财经大学 统计学系
主要内容
统 • 单总体问题的符号检验 计 • 两总体问题的符号检验
• 威尔科克森配对符号秩检验 学
天津财经大学 统计学系
一、单总体问题的符号检验
• 单总体符号检验适用于检验总体中位数 统 是否在某一指定位置。
• 检验时,可根据样本中正号的数目来决 计 定是否拒绝原假设:

SAS的非参数检验

SAS的非参数检验

SAS的非参数检验非参数检验是一种统计方法,用于处理数据不满足正态分布或方差齐性的情况。

它们不依赖于任何概率分布的假设,因此也被称为非参数检验。

SAS(统计分析系统)是一种常用的统计软件,提供了多种非参数检验方法。

本文将介绍一些常见的非参数检验方法及其在SAS中的应用。

1. Wilcoxon符号秩检验(Wilcoxon Signed Rank Test):Wilcoxon符号秩检验是一种用于比较两个相关样本或配对样本的非参数检验方法。

它对于数据不满足正态分布的情况非常有用。

它的原假设是两个样本的中位数不同。

在SAS中,可以使用PROC UNIVARIATE来执行Wilcoxon符号秩检验。

下面是一个示例代码:```proc univariate data=mydata;var x1 x2;wilcoxon signedrank;run;```其中,mydata是数据集名称,x1和x2是要比较的两个变量。

wilcoxon signedrank选项告诉SAS执行Wilcoxon符号秩检验。

2. Mann-Whitney U检验(Mann-Whitney U Test):Mann-Whitney U检验是一种用于比较两个独立样本的非参数检验方法。

它的原假设是两个样本的总体分布相同。

在SAS中,可以使用PROC NPAR1WAY来执行Mann-Whitney U检验。

下面是一个示例代码:```proc npar1way data=mydata;var x;class group;mannwhitney u(x) / wilcoxon;run;```其中,mydata是数据集名称,x是要比较的变量,group是分组变量。

mannwhitney u选项告诉SAS执行Mann-Whitney U检验。

3. Kruskal-Wallis检验(Kruskal-Wallis Test):Kruskal-Wallis检验是一种用于比较三个或更多独立样本的非参数检验方法。

Wilcoxon符号秩检验

Wilcoxon符号秩检验

第二节Wilcoxon符号秩检验Wilcoxon符号秩检验符号检验只用了差的符号,但没有利用差值的大小。

12 3Wilcoxon符号秩检验(Wilcoxon signed-rank test) 把差的绝对值的秩分别按照不同的符号相加作为其检验统计量。

显然,相比较于符号检验,Wilcoxon符号秩检验利用了更多的信息。

Wilcoxon符号秩检验:条件u Wilcoxon符号秩检验需要一点总体分布的性质;它要求假定样本点来自连续对称总体分布;而符号检验不需要知道任何总体分布的性质。

u在对称分布中,总体中位数和总体均值是相等的;因此,对于总体中位数的检验,等价于对于总体均值的检验。

u Wilcoxon符号秩检验实际是对对称分布的总体中位数(或均值)的检验。

Wilcoxon符号秩检验:基本原理u计算差值绝对值的秩。

u分别计算出差值序列里正数的秩和(W+)以及负数的秩和(W-)。

u如果原假设成立,W+与W-应该比较接近。

如果W+和W-过大或过小,则说明原假设不成立。

u将正数的秩和或者负数的秩作为检验统计量,根据其统计分布计算p值,从而可以得出检验的结论。

具体步骤设定原假设和备择假设。

分别计算出差值序列中正数的秩和W+以及负数的秩和W-。

根据W+和W-建立检验统计量,计算p值并得出检验的结论。

在双侧检验中检验统计量可以取为W=min(W+,W-)。

显然,如果原假设成立,W+与W-应该比较接近。

如果二者过大或过小,则说明原假设不成立。

秩的计算注意问题计算差值绝对值的秩时,注意差值等于0值不参与排序。

下面一行R i就是上面一行数据Z i的秩。

Z i159183178513719 R i75918426310数据中相同的数值称为“结”。

结中数字的秩为它们所占位置的平均值Z i159173178513719 R i758.518.5426310关于P值u有了检验统计量W,我们就可根据其统计分布计算p值了,双侧检验的p值等于,式中w为检验统计量的样本观测值。

Wilcoxon符号秩检验

Wilcoxon符号秩检验

第二节Wilcoxon符号秩检验Wilcoxon符号秩检验符号检验只用了差的符号,但没有利用差值的大小。

12 3Wilcoxon符号秩检验(Wilcoxon signed-rank test) 把差的绝对值的秩分别按照不同的符号相加作为其检验统计量。

显然,相比较于符号检验,Wilcoxon符号秩检验利用了更多的信息。

Wilcoxon符号秩检验:条件u Wilcoxon符号秩检验需要一点总体分布的性质;它要求假定样本点来自连续对称总体分布;而符号检验不需要知道任何总体分布的性质。

u在对称分布中,总体中位数和总体均值是相等的;因此,对于总体中位数的检验,等价于对于总体均值的检验。

u Wilcoxon符号秩检验实际是对对称分布的总体中位数(或均值)的检验。

Wilcoxon符号秩检验:基本原理u计算差值绝对值的秩。

u分别计算出差值序列里正数的秩和(W+)以及负数的秩和(W-)。

u如果原假设成立,W+与W-应该比较接近。

如果W+和W-过大或过小,则说明原假设不成立。

u将正数的秩和或者负数的秩作为检验统计量,根据其统计分布计算p值,从而可以得出检验的结论。

具体步骤设定原假设和备择假设。

分别计算出差值序列中正数的秩和W+以及负数的秩和W-。

根据W+和W-建立检验统计量,计算p值并得出检验的结论。

在双侧检验中检验统计量可以取为W=min(W+,W-)。

显然,如果原假设成立,W+与W-应该比较接近。

如果二者过大或过小,则说明原假设不成立。

秩的计算注意问题计算差值绝对值的秩时,注意差值等于0值不参与排序。

下面一行R i就是上面一行数据Z i的秩。

Z i159183178513719 R i75918426310数据中相同的数值称为“结”。

结中数字的秩为它们所占位置的平均值Z i159173178513719 R i758.518.5426310关于P值u有了检验统计量W,我们就可根据其统计分布计算p值了,双侧检验的p值等于,式中w为检验统计量的样本观测值。

威尔可森符号秩检验

威尔可森符号秩检验

威尔可森符号秩检验威尔科克森符号秩检验(Wilcoxon signed-rank test)是一种非参数统计方法,用于比较成对样本的差异。

它基于样本数据的符号秩来进行推断。

以下是威尔科克森符号秩检验的基本步骤:1、假设检验:●零假设(H0):成对样本之间没有差异(即两个样本的中位数相等)。

●对立假设(H1):成对样本之间存在差异(即两个样本的中位数不相等)。

2、计算差异:●对每对成对样本计算差异。

●将这些差异按照绝对值大小进行排序,并为每个差异分配一个符号秩(正负号),如果有相同的差异,则取平均秩。

3、计算符号秩和:分别计算正符号秩和负符号秩的总和。

4、计算检验统计量:使用计算得到的正负符号秩和,计算检验统计量W。

5、根据检验统计量W进行假设检验:●对于小样本(n<30),可以使用查表法或精确法确定临界值,以判断是否拒绝零假设。

●对于大样本,可以使用正态近似法(z检验)进行假设检验。

威尔科克森符号秩检验用于成对样本的非参数分析,并且不要求数据满足正态分布假设。

它适用于样本大小较小或无法满足正态分布假设的情况下使用。

在Matlab中,可以使用signrank函数执行威尔科克森符号秩检验。

以下是一个示例:matlab% 假设有两组成对样本数据group1 = [5, 7, 9, 11, 13];group2 = [4, 6, 10, 12, 14];% 进行威尔科克森符号秩检验[p, h, stats] = signrank(group1, group2);% 显示结果disp(['p值:', num2str(p)]);if hdisp('拒绝零假设');elsedisp('接受零假设');enddisp(['检验统计量W:', num2str(stats.signedrank)]);disp(['样本大小n:', num2str(stats.n)]);在上述示例中,我们假设有两组成对样本数据group1 和group2,并使用signrank 函数进行威尔科克森符号秩检验。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

第二十七课 符号检验和Wilcoxon 符号秩检验在统计推断和假设检验中,传统的检验统计量都叫做参数检验,因为它们都依赖于确定的概率分布,这个分布带有一组自由的参数。

参数检验被认为是依赖于分布假定的。

通常情况下,我们对数据进行分析时,总是假定误差项服从正态分布,这是人们易于接受的事实,因为正态分布的原始出发点就是来自于误差分布,至于当样本相当大时,数据的正态近似,这是由于大样本理论所保证的。

但有些资料不一定满足上述要求,或不能测量具体数值,其观察结果往往只有程度上的区别,如颜色的深浅、反应的强弱等,此时就不适用参数检验的方法,而只能用非参数统计方法(non-parametric statistical analysis )来处理。

这种方法对数据来自的总体不作任何假设或仅作极少的假设,因此在实用中颇有价值,适用面很广。

一、 单样本的符号检验符号检验(sign test )是一种最简单的非参数检验方法。

它是根据正、负号的个数来假设检验。

首先需要将原始观察值按设定的规则,转换成正、负号,然后计数正、负号的个数作出检验。

该检验可用于样本中位数和总体中位数的比较,数据的升降趋势的检验,特别适用于总体分布不服从正态分布或分布不明的配对资料,有时当配对比较的结果只能定性的表示,如试验前后比较结果为颜色从深变浅、程度从强变弱,成绩从一般变优秀,即不能获得具体数字,也可用符号检验,例如用正号表示颜色从深变浅,用负号表示颜色从浅变深。

用于配对资料时,符号检验的计算步骤为:首先定义成对数据指定正号或负号的规则,然后计数正号的个数+S 及负号的个数-S ,由于在具体比较配对资料时,可能存在配对资料的前后没有变化,或等于假设中的中位数,此时仅需要将这些观察值从资料中剔除,当然样本大小n 也随之减少,故修正样本大小-++=S S n 。

当样本n 较小时,应使用二项分布确切概率计算法,当样本n 较大时,常利用二项分布的正态近似。

1. 小样本时的二项分布概率计算当20≤n 时,+S 或-S 的检验p 值由精确计算尺度二项分布的卷积获得。

在比较配对资料试验前后有否变化,或增加或减小的假设检验时,如果我们定义试验后比试验前增加为正号,反之为负号,那么对于原假设:试验前后无变化来说,正号的个数+S 和负号的个数-S 可能性应当相等,即正号出现的概率p =0.5,于是+S 与-S 均服从二项分布)5.0,(n B ,对于太大的+S 相应太小的-S ,或者太大的-S 相应太小的+S ,都将拒绝接受原假设;对于原假设:试验后比试验前有增加来说,正号的个数+S 大于负号的个数-S 的可能性应该大,即正号出现的概率5.0>p ,对于太小的+S 相应太大的-S ,将拒绝接受原假设;对于原假设:试验后比试验前减小来说,正号的个数+S 小于等于负号的个数-S 的可能性应该大,即正号出现的概率5.0≤p ,对于太大的+S 相应太小的-S ,将拒绝接受原假设。

例27.1有一种提高学生某种素质的训练,有人说它是无效的,有人说它是有效的,那么真实情况究竟应该是怎样的呢?随机地选取15名学生作为试验样本,在训练开始前做了一次测验,每个学生的素质按优、良、中、及、差打分,经过三个月训练后,再做一次测试对每个学生打分。

数据见表27.1所示。

我们将素质提高用正号表示,反之用负号表示,没有变化用0表示。

显著性水平取0.1。

表27.1 训练前后的素质比较学生编号训练之前 训练之后 差异符号 1 中 优 + 2 及 良 + 3 良 中 - 4 差 中 + 5 良 良 0 6 中 优 + 7 差 及 + 8 良 优 + 9 中 差 - 10 差 中 + 11 中 优 + 12 及 良 + 13 中 及 - 14 中 优 + 15差中+从表27.1中15名学生训练前后的差异分析可得出:有14名学生有差异,其中+S =11,-S =3。

1名学生无差异(学生编号为5),应该从分析中去掉,所以n =15-1=14。

假设检验为:5.0:0≤p H 即训练之后学生素质没有提高。

5.0:1>p H 即训练之后学生素质有提高。

由于试验的结果只有两种可能,正号或负号,对每一个学生试验出现正号的假定概率为p =0.5,负号为1—p =0.5,这样整个试验的概率是相同的,并且每一个试验是相互独立的。

因此在n =14次独立的试验中,正号出现的次数服从二项分布)5.0,14(B ,见表27.2所示。

表27.2 二项分布的概率和累计概率n =14,p =0.5正号出现的次数正号出现的概率累计概率 0 0.0001 0.0001 1 0.0009 0.0009 2 0.0056 0.0065 30.02220.02874 0.0611 0.08985 0.1222 0.21206 0.1833 0.39537 0.2095 0.60478 0.1833 0.78809 0.1222 0.9102 10 0.0611 0.9713 11 0.0222 0.9935 12 0.0056 0.9991 13 0.0009 0.9999 140.00011.0000从表27.2的累计概率列中我们看到,正号出现的次数大于10的概率为1-0.9713=0.0287,或者换一种方法计算为=0.0001+0.0009+0.0056+0.0222=0.0287,二者的微小差异是因为小数点后舍入问题造成的。

而试验的结果:正号出现的次数为11,大于10,出现的概率不会超过0.0287,我们开始设定的显著性水平为0.1,由于0.0287<0.1,所以我们拒绝原假设,接受备选假设。

如果我们的原假设为p =0.5,既训练前后学生素质相等,那么就是双侧检验,应该加上正号出现的次数小于4的概率0.0287,即2×0.0287=0.0574<0.1,同样是拒绝原假设,接受区间为4次到10次,而拒绝区间为小于等于3次(小于4次)或大于等于11次(大于10 次)。

2. 大样本时的正态近似概率计算当20>n 时,样本可以认为是大样本。

我们可以利用二项分布的正态近似,即对于),(~p n B S ,二项分布的期望均值为np ,方差为)1(p np -,当n 比较大时,且np 和)1(p n -大于5,可以近似地认为)1,0(~)1(N p np np S z --=(27.1)公式中的S 表示正号或者负号的个数,符号检验时,p =0.5代入(27.1)式中,得到大样本时的正态近似统计量)1,0(~5.05.0N nn S z -=(27.2)当S >2/n 时,应该修正S 为S -0.5;当S <2/n 时,应该修正S 为S +0.5。

S 值加或减的0.5是连续性修正因子,目的是为了能将连续分布应用到近似的离散型分布。

二、 配对资料的Wilcoxon 符号秩检验当两组配对资料近似服从正态分布,它们差值的检验可以使用配对t 检验法。

如果配对资料的正态分布的假设不能成立,就可以使用Frank Wilcoxon (1945)符号秩检验,它是一种非参数检验方法,对配对资料的差值采用符号秩方法来检验。

它的基本要求是差值数据设置为最小的序列等级和两组配对资料是相关的(配成对)。

在两组配对资料的差异有具体数值的情况下,符号检验只利用大于0和小于0的信息,即正号和负号的信息,而对差异大小所包含的信息却未加利用,但Wilcoxon 符号秩检验方法既考虑了正、负号,又利用了差值大小,故效率较符号检验法高。

例27.2某制造商想要比较两种不同的生产方法所花费的生产时间是否有差异。

随机地选取了11个工人,每一个工人都分别使用两种不同的生产方法来完成一项相同的任务,每一个工人开始选用的生产方法是随机的,即可以先使用生产方法1再使用生产方法2,也可以先用生产方法2再使用生产方法1。

这样,在样本中的每一个工人都提供了一个配对观察。

数据见表27.3所示。

任务完成时间的正差值表示生产方法1需要更多的时间,负差值表示生产方法2需要更多的时间。

表27.3 两种不同生产方法完成任务的时间(分钟) 工人编号n生产方法M 差值D 绝对差值 秩次 R 符号秩次R M 1 M 2 D =M 1-M 2|D | - + 1 10.2 9.5 0.7 0.7 8 8 2 9.6 9.8 -0.2 0.2 2 2 3 9.2 8.8 0.4 0.4 3.5 3.5 4 10.6 10.1 0.5 0.5 5.5 5.5 5 9.9 10.3 -0.4 0.4 3.5 3.5 6 10.2 9.3 0.9 0.9 10 10 7 10.6 10.5 0.1 0.1 1 1 8 10.0 10.0 0 0 — — — 9 11.2 10.6 0.6 0.6 7 7 10 10.7 10.2 0.5 0.5 5.5 5.5 1110.69.80.80.899 符号秩次总和-T =5.5,+T =49.55.549.5为了比较两种方法的任务完成时间是否有显著差异,假设检验为::0H 任务完成时间的两个总体是相同的。

:1H 任务完成时间的两个总体是不相同的。

使用Wilcoxon 符号秩检验方法的主要步骤见表27.3中每列的计算方法和过程,先求出每对数据的差值D ,按差值绝对值|D |由小到大排列并给秩R ,从秩1开始到秩10,注意工人编号为8的配对数据,由于差值为0,在排秩中丢弃,样本数目修正为n =11-1=10。

在给秩值时,遇到相等|D |,也称为结值(tied ),使用平均秩,如工人编号3和5具有相同的绝对差值0.4,所以平分秩3和秩4,各为秩3.5。

一旦绝对差值的秩值R 给出后,然后将R 分成正和负差值的两个部分秩值+R 和-R ,最后求符号秩和∑++=RT ,∑--=RT ,如-T =2+3.5=5.5。

对于样本数目有n 个,+T 与-T 的最小可能值为0,而最大可能值为(1+2+…+n )=n (n +1)/2。

显然,应当有+T +-T = n (n +1)/2,如本例5.5+49.5=55=10(10+1)/2。

那么符号秩的平均值为n (n +1)/4。

构造Wilcoxon 符号秩统计量为4)1(+-=+n n T S (27.3)显然如果原假设为真,+T 与-T 应该有相同的值,等于n (n +1)/4,因此太大的S 值或太小的S 值都是我们拒绝的依据。

在实际工作中便于计算常取W=min (+T ,-T ),W 服从所谓的Wilcoxon 符号秩分布,对于本例n =10,=S 49.5-10(10+1)/4=22,W = min (49.5,5.5)=5.5,查表可得在显著水平=α0.05,n =10的双侧检验的临界值为8,即W 值的拒绝区域为0到8,接受区域为8到27.5。

相关文档
最新文档