配对设计秩和检验
配对设计的情况:(1)同对的受试者分别给予两种处理,如同窝别、体重相近的两只动物配对;(2)同一受试对象处理前后的比较。
配对设计资料一般采用配对t检验,但若配对资料的差值不满足正态分布,则采用非参数的方法进行检验。
在H0成立时,如果样本比较多,正差值和负差值的秩和理论上应相等,即使有差别也是因随机因素造成的。换言之,如果H0成立,随机样本出现正差值和负差值的秩和悬殊较大的可能性比较小,如果悬殊比较大,则我们有理由拒绝H0,接受H1,则认为处理效应不同。
一、检验步骤说明
1、 建立检验假设,确定检验水准
H0:差值的总体中位数等于0
H1:差值的总体中位数不等于0
α=0.05
2、 计算检验统计量T值
(1) 求差值
(2) 编秩
(3) 分别求正、负秩和,记为T+、T-
(4) T值=min(T+,T-)(双侧检验。取绝对值较小者)
3、 统计推断
(1) n<=0,查T界值表
(2) T分布以均值为中心对称的非连续分布。随着n的增大,T的分布逐渐逼近均数为n(n+1)/4,方差为n(n+1)(2n+n)/24的正态分布。n>25,T分布近似正态分布。若n>50,则可用Z检验。当相同差值较多时,采用校正的Z检验。
P<0.05,拒绝H0,接受H1;反之则,不拒绝H0。
二、结果
三、代码
data a; /创建数据集,配对设计资料/
input x1 x2 ;
d = x1 - x2;
cards;
39 55
42 54
51 55
43 47
55 53
45 63
22 52
48 44
40 48
45 55
40 32
49 57
;
run;
proc univariate data = a; /检验/
var d;
run;