数学建模社区-数学中国

标题: 多个样本的非参数检验的两两比较 [打印本页]
作者: 厚积薄发    时间: 2010-8-14 19:57
标题: 多个样本的非参数检验的两两比较
由于各个样本的总体呈偏态分布或方差不齐,无法使用方差分析中的两两比较，检验其总体分布是否相同，常用的非参数秩和检验方法是Kruskal-Wallis 法,在SAS 软件中实现的过程步有以下三种方法: NPAR1 WAY 过程、FREQ过程以及RAN K 和ANOVA 两过程的结合，而样本间两两比较则可以通过RAN K 和ANOVA 两过程的结合,采用MEANS 语句来实现，即先对原始数据进行排序，求相应的秩次，然后对秩进行参数的多重比较。

# T* n5 @1 E$ g# i8 r/ J2 l3 L
   调用FREQ 过程产生的第二个CMH统计量” row mean scores differ” (Kruskal-Wallis结果)、NPAR1 WAY 过程产生的卡方统计量以及ANOVA 过程产生的R-Square 与υ总(总自由度)之积,即为Kruskal-Wallis 检验结果。本文种ANOVA 过程步中的MEANS 规定采用LSD 法进行两两比较，也可使用其他方法。
data sample ;
do group = 1 to 3 ;
input x @@;
output ;
end ;
cards ;
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 ;
tables group*x/ scores = rank cmh2 noprint;
run;
proc npar1way wilcoxon;
class group;
var x;
run;
proc rank data =sample out = a ;
var x;
ranks r;
proc anova;
class group ;
model r = group ;
means group/ lsd snk ;
run;
quit;

  ^+ H$ T; a+ R( g
  成组设计的等级资料或频数表资料多个样本及其两两比较
  这种类型的资料为成组设计的等级资料或频数表资料,但频数表资料时各个样本的总体呈偏态分布或方差不齐,检验其总体分布是否相同，常用的非参数秩和检验方法也是Kruskal-Wallis 法,在SAS 软件中实现的过程步同上,程序略有差别。
data sample2;
do effect = 1 to 4 ;
do group = 1 to 3 ;
input f @@;
output ;
end ;
end ;
cards ;
10 24 48 17 41 65 19 33 36 4 7 8
;
proc freq;
tables group*effect/ scores = rank cmh2 noprint;
weight f;
0 }* J# V4 O  Hrun;
data a;
set sample2;
do i = 1 to f ;
output;
end ;
run;
proc npar1way wilcoxon data=a;
var effect;
class group;
run;
proc rank data=a out=b;
var effect ;
ranks r;
run;
proc glm;
class group;
model r = group ;
means group/lsd snk tukey dunnett ;
run;
quit;

- v8 c/ o2 G& G7 J: M

0 X: A) N  t% M2 H$ E9 t5 G随机区组设计的定量资料多个样本及其两两比较
各个样本的总体呈偏态分布或方差不齐,检验其总体分布是否相同,常用的非参数秩和检验方法是Friedman 法,在SAS 软件中实现的过程步有以下两种方法: FREQ 过程以及RAN K和GLM 两过程的结合。各个样本两两比较可通过RAN K 和GLM 两过程的结合, 采用LSMEANS 语句来实现

; B8 a  @2 H8 U8 k: q3 O( r调用FREQ 过程产生的第二个CMH统计量和GLM 过程产生的R2 与(υ总-υ区组) 之积,即
为Friedman 检验结果。LSMEANS 语句规定对调整后的样本均数进行两两比较。
data sample3;
do block = 1 to 7 ;
do treat = 1 to 4 ;
input x @@;
output;
end;
end ;
cards ;
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 ;
tables block*treat*x/ scores = rank cmh2 noprint;
run;
proc sort data=sample3 out=a;
by block ;
run;
0 e0 y+ G6 _# l; Z: T( k
proc rank data=a out=b;
var x;
by block;
ranks r ;
run;
proc glm data=b ;
class treat block ;
model r = treat block/ ss1 ;
lsmeans treat/ stderr pdiff ;
run;
quit;
! H- h9 W- b0 z( U5 t
来源：http://www.shumo.cn/item/show.asp?m=1&d=699
% `# B# M% W. z6 C/ S% q
作者: 紫辰    时间: 2010-8-14 23:50
版主很好的啊支持你！
作者: foreveringxq    时间: 2010-8-15 07:38
支持~~~~~~~~~~~~~~~
作者: shif    时间: 2010-8-15 08:59
ding~~~ding
作者: shif    时间: 2010-8-15 08:59
zaiding!dfgsdfg
作者: shif    时间: 2010-8-15 09:00
这个帖子真的很好~激动~
作者: yitaosc    时间: 2010-8-15 10:04
ding ^^^^^^^^^^^^^^^^^66
作者: lifelight    时间: 2010-8-26 11:33
偶的天啊!这是什么样的贴子呢？
作者: 帮忙    时间: 2010-8-26 11:34
楼主，你写得实在是太好了。我惟一能做的，就只有把这个帖子顶上去这件事了
作者: lsrj    时间: 2010-8-26 11:35
试试运气啦~~~~~~~~~~~
作者: ApolloYang    时间: 2010-8-26 11:36
我要把这个帖子一直往上顶，往上顶！
作者: ksyy024    时间: 2010-8-26 11:38
我要把这个帖子一直往上顶，往上顶！
作者: rancid    时间: 2010-8-26 11:45
不错不错,我喜欢看  
作者: xiaoyang1115    时间: 2010-8-26 13:52
我来了~~~~~~~~~ 闪人~~~~~~~~~~~~~~~~  
作者: fangruiscience    时间: 2010-8-26 18:58
我来了~~~~~~~~~ 闪人~~~~~~~~~~~~~~~~  
作者: li    时间: 2010-8-26 19:16
我要把这个帖子一直往上顶，往上顶！
作者: banana001    时间: 2010-8-27 00:00
哦～～
作者: as    时间: 2010-8-27 08:00
楼主，你写得实在是太好了。我惟一能做的，就只有把这个帖子顶上去这件事了
作者: tiger123    时间: 2010-8-27 12:00
我要把这个帖子一直往上顶，往上顶！
作者: whrichardliu    时间: 2010-8-27 15:00
楼主的帖子实在是写得太好了。可是我立刻想到，这么好的帖子，倘若别人看不到，那么不是浪费楼主的心血吗？经过痛苦的思想斗争，我终于下定决心，牺牲小我，奉献大我。我要拿出这帖子奉献给世人赏阅，我要把这个帖子一直往上顶，往上顶！顶到所有人都看到为止！  
作者: TomYu    时间: 2010-8-27 20:00
哦～～
作者: icq0607    时间: 2010-8-28 08:00
呵呵 大家好奇嘛 来观看下~~~~  
作者: 米菲兔Laura    时间: 2010-8-28 10:48
？？没看懂~
作者: inforworm    时间: 2010-8-28 12:00
留个脚印```````
作者: yuwengchen    时间: 2010-8-28 15:00
哦～～
作者: 皮皮雅    时间: 2010-8-28 20:00
我要把这个帖子一直往上顶，往上顶！
作者: keaianzi    时间: 2010-8-28 23:59
顶顶更健康，越顶吃的越香。
作者: douzi    时间: 2010-8-29 08:00
楼主，你写得实在是太好了。我惟一能做的，就只有把这个帖子顶上去这件事了
作者: phenix    时间: 2010-8-29 12:00
偶的天啊!这是什么样的贴子呢？
作者: hqmailes    时间: 2010-8-29 15:00
强人，佩服死了。呵呵，不错啊
作者: 西海岸的帆    时间: 2010-8-29 20:00
哦～～
作者: huanglonghaoliu    时间: 2010-8-30 08:00
试试运气啦~~~~~~~~~~~
作者: frilin    时间: 2010-8-30 12:00
呵呵 大家好奇嘛 来观看下~~~~  
作者: dlight    时间: 2010-8-30 15:00
声明一下：本人看贴和回贴的规则，好贴必看，精华贴必回。
作者: 太玄    时间: 2010-8-30 20:00
哦～～
作者: kunnba    时间: 2010-8-31 08:00
楼主那种裂纸欲出的大手笔，竟使我忍不住一次次的翻开楼主的帖子……   
作者: zjjweiwei88    时间: 2010-8-31 12:00
留个脚印```````
作者: benzheren    时间: 2010-8-31 15:00
楼主，你写得实在是太好了。我惟一能做的，就只有把这个帖子顶上去这件事了
作者: brighter    时间: 2010-8-31 20:00
我要把这个帖子一直往上顶，往上顶！
作者: huangzhu54    时间: 2010-9-1 08:00
(*^__^*) 指点系词……激扬文字……  
作者: xhw2007    时间: 2010-9-1 12:00
我回不抢呢 考虑再三 还是不抢了吧 ^_^
作者: aaaaaaaa    时间: 2010-9-1 12:00
来报道!!!!!!!!!!!
作者: hongda1227    时间: 2010-9-1 15:00
我要把这个帖子一直往上顶，往上顶！
作者: lxy519    时间: 2010-9-1 20:00
楼主的帖子实在是写得太好了。可是我立刻想到，这么好的帖子，倘若别人看不到，那么不是浪费楼主的心血吗？经过痛苦的思想斗争，我终于下定决心，牺牲小我，奉献大我。我要拿出这帖子奉献给世人赏阅，我要把这个帖子一直往上顶，往上顶！顶到所有人都看到为止！  
作者: loveaixing    时间: 2010-9-2 08:00
看起来好~~像啊~~~~~
作者: linqi2004    时间: 2010-9-2 12:00
我来了~~~~~~~~~ 闪人~~~~~~~~~~~~~~~~  
作者: maoloverme    时间: 2010-9-2 15:00
哦～～
作者: bing    时间: 2010-9-2 20:00
楼主，你写得实在是太好了。我惟一能做的，就只有把这个帖子顶上去这件事了
作者: 小痒    时间: 2010-9-3 08:00
强人，佩服死了。呵呵，不错啊
作者: tkttsai    时间: 2010-9-3 12:00
试试运气啦~~~~~~~~~~~
作者: eflying    时间: 2010-9-3 15:00
鉴定完毕！  
作者: bourbaki    时间: 2010-9-3 20:00
顶顶更健康，越顶吃的越香。
作者: lgauss    时间: 2010-9-4 08:00
强人，佩服死了。呵呵，不错啊
作者: jmwen    时间: 2010-9-4 12:00
(*^__^*) 指点系词……激扬文字……  
作者: darkuard    时间: 2010-9-4 15:00
声明一下：本人看贴和回贴的规则，好贴必看，精华贴必回。
作者: zhouhping    时间: 2010-9-4 20:00
留个脚印```````
作者: uoykcah    时间: 2010-9-5 12:00
留个脚印```````
作者: 4402103    时间: 2010-9-5 15:00
我要把这个帖子一直往上顶，往上顶！
作者: 54hero007    时间: 2010-9-5 20:00
强烈支持。楼主万岁
作者: kekonglin    时间: 2010-9-6 08:00
留个脚印```````
作者: fire007sha    时间: 2010-9-6 12:00
强人，佩服死了。呵呵，不错啊
作者: wyl    时间: 2010-9-6 15:00
声明一下：本人看贴和回贴的规则，好贴必看，精华贴必回。
作者: 东风21    时间: 2010-9-6 20:00
我要把这个帖子一直往上顶，往上顶！
作者: xxxx    时间: 2010-9-7 08:00
声明一下：本人看贴和回贴的规则，好贴必看，精华贴必回。
作者: zbh228    时间: 2010-9-7 12:00
看起来好~~像啊~~~~~
作者: kevin_nit    时间: 2010-9-7 15:00
试试运气啦~~~~~~~~~~~
作者: matdrchen    时间: 2010-9-7 20:00
我来了~~~~~~~~~ 闪人~~~~~~~~~~~~~~~~  
作者: meseeme    时间: 2010-9-8 08:00
楼主那种裂纸欲出的大手笔，竟使我忍不住一次次的翻开楼主的帖子……   
作者: pcm    时间: 2010-9-8 12:00
留个脚印```````
作者: class77787    时间: 2010-9-8 15:00
顶顶更健康，越顶吃的越香。
作者: wenfeng018    时间: 2010-9-8 20:00
我要把这个帖子一直往上顶，往上顶！
作者: kamma    时间: 2010-9-9 08:00
(*^__^*) 指点系词……激扬文字……  
作者: dddddddd    时间: 2010-9-9 08:00
我基本上是采用看英语文章的办法，先泛读，再精读，再一句一句看，最后再提纲挈领，总算是明白一点了，当然，也可能还是领悟错了。最后要说的一句话是：楼主，你很牛叉，希望你不是真的有病。   
作者: julius2004    时间: 2010-9-9 12:00
试试运气啦~~~~~~~~~~~
作者: yuro    时间: 2010-9-9 15:00
试试运气啦~~~~~~~~~~~
作者: zys402    时间: 2010-9-9 20:00
试试运气啦~~~~~~~~~~~
作者: 两个人    时间: 2010-9-9 20:00
我基本上是采用看英语文章的办法，先泛读，再精读，再一句一句看，最后再提纲挈领，总算是明白一点了，当然，也可能还是领悟错了。最后要说的一句话是：楼主，你很牛叉，希望你不是真的有病。   
作者: pgfzhy    时间: 2010-9-10 11:40
来报道!!!!!!!!!!!
作者: amocilin    时间: 2010-9-10 11:52
我要把这个帖子一直往上顶，往上顶！
作者: xuyisen    时间: 2010-9-10 12:00
留个脚印```````
作者: BMSHAO    时间: 2010-9-10 15:00
我基本上是采用看英语文章的办法，先泛读，再精读，再一句一句看，最后再提纲挈领，总算是明白一点了，当然，也可能还是领悟错了。最后要说的一句话是：楼主，你很牛叉，希望你不是真的有病。   
作者: 双子座    时间: 2010-9-10 20:00
声明一下：本人看贴和回贴的规则，好贴必看，精华贴必回。
作者: liuwei    时间: 2010-9-11 15:01
强烈支持。楼主万岁
作者: jixue    时间: 2010-9-11 20:00
哦～～
作者: chenzhiyongdlut    时间: 2010-9-12 15:02
声明一下：本人看贴和回贴的规则，好贴必看，精华贴必回。
作者: likaiwudi    时间: 2010-9-12 20:00
声明一下：本人看贴和回贴的规则，好贴必看，精华贴必回。
作者: hiv    时间: 2010-9-13 15:00
试试运气啦~~~~~~~~~~~
作者: 半生蕉    时间: 2012-8-14 12:48
不错。。。。。
作者: 杨小丁92    时间: 2013-7-18 20:43




欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5