多个样本的非参数检验的两两比较[教育]
多个样本的非参数检验的两两比较 [SS]
2013-02-23 17:26:31
由于各个样本的总体呈偏态分布或方差不齐,无法使用方差分析中的两两比较,检验其总体分布是否相同,常用的非参数秩和检验方法是Kruskl-Wllis
法,在SS 软件中实现的过程步有以下三种方法: NPR1 WY 过程、FREQ过程以及RN K 和NOV 两过程的结合,而样本间两两比较则可以通过
RN K 和NOV 两过程的结合,采用MENS 语句来实现,即先对原始数据进行排序,求相应的秩次,然后对秩进行参数的多重比较。
调用FREQ 过程产生的第二个CMH统计量” row men scores differ” (Kruskl-Wllis结果)、NPR1 WY 过程产生的卡方统计量以及NOV 过程产生的
R-Squre 与υ总(总自由度)之积,即为Kruskl-Wllis 检验结果。本文种NOV 过程步中的MENS 规定采用LSD 法进行两两比较,也可使用其他方
法。
dt smple ;
do group = 1 to 3 ;
input x @@;
output ;
end ;
crds ;
9.8 0.6 0.4 10.2 1.2 1.9 10.6 2.0 2.2 13.0 2.4 2.5 14.0 3.1 2.8 14.8 4.1 3.1 15.6 5.0 3.7 15.6 5.9 3.9 21.6 7.4 4.6 24.0 13.6 7.0 ;
proc freq ;
tbles group*x/ scores = rnk cmh2 noprint; run;
proc npr1wy wilcoxon;
clss group;
vr x;
run;
proc rnk dt =smple out = ;
vr x;
rnks r;
proc nov;
clss group ;
model r = group ;
mens group/ lsd snk ;
run;
quit;
成组
设计
领导形象设计圆作业设计ao工艺污水处理厂设计附属工程施工组织设计清扫机器人结构设计
的等级资料或频数表资料多个样本及其两两比较
这种类型的资料为成组设计的等级资料或频数表资料,但频数表资料时各个样本的总体呈偏态分布或方差不齐,检验其总体分布是否相同,常用的非
参数秩和检验方法也是Kruskl-Wllis 法,在SS 软件中实现的过程步同上,程序略有差别。 dt smple2;
do effect = 1 to 4 ;
do group = 1 to 3 ;
input f @@;
output ;
end ;
end ;
crds ;
10 24 48 17 41 65 19 33 36 4 7 8 ;
proc freq;
tbles group*effect/ scores = rnk cmh2 noprint;
weight f;
run;
dt ;
set smple2;
do i = 1 to f ;
output;
end ;
run;
proc npr1wy wilcoxon dt=; vr effect;
clss group;
run;
proc rnk dt= out=b;
vr effect ;
rnks r;
run;
proc glm;
clss group;
model r = group ;
mens group/lsd snk tukey dunnett ; run;
quit;
随机区组设计的定量资料多个样本及其两两比较
各个样本的总体呈偏态分布或方差不齐,检验其总体分布是否相同,常用的非参数秩和检验方法是Friedmn 法,在SS 软件中实现的过程步有以下两
种方法: FREQ 过程以及RN K和GLM 两过程的结合。各个样本两两比较可通过RN K 和GLM 两过程的结合, 采用LSMENS 语句来实现
调用FREQ 过程产生的第二个CMH统计量和GLM 过程产生的R2 与(υ总-υ区组) 之积,即
为Friedmn 检验结果。LSMENS 语句规定对调整后的样本均数进行两两比较。 dt smple3;
do block = 1 to 7 ;
do tret = 1 to 4 ;
input x @@;
output;
end;
end ;
crds ;
63 188 138 54 90 238 220 144 54 300 83 92 45 140 213 100 54 175
150 36 72 300 163 90 64 207 185 87
;
proc freq ;
tbles block*tret*x/ scores = rnk cmh2 noprint;
run;
proc sort dt=smple3 out=; by block ;
run;
proc rnk dt= out=b;
vr x;
by block;
rnks r ;
run;
proc glm dt=b ;
clss tret block ;
model r = tret block/ ss1 ; lsmens tret/ stderr pdiff ; run;
quit;