多样本均值比较

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

我有五个样本，用group识别，每个样本都有一个变量x，我要比较这五个样本的x总体是否有显著差异，最好还能比较每两个样本之间是否有差异，请问各位大牛sas有什么proc 过程可以实现，最好能帮我把具体程序贴出来。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

关键词：均值比较 Group 具体程序 两个样本 ROC 样本 程序 最好

1. proc glm;
   
2.   class group;
   
3.   model x=group;
   
4.   means group /snk bon sidak  hovtest=levene(type=abs)  hovtest=bartlett;
   
5. run;

复制代码

楼上正解，另外可用proc anova ，具体与glm的区别在于处理的组是否平衡

对单因素2组或者多组定量资料进行假设检验时,

首先想到的也许是ANOVA,

但对独立样本进行ANOVA时,

1.要检验各组数据是否满足正态性;

2.要检验各组数据方差是否齐;



当这2个条件不满足或不完全满足时,

可以采用秩和检验;



组数为3组或者3组以上时,

若Kruskal-Wallis H检验结果有统计学意义,

*若没有相同秩次时,无需对服从均匀分布的H统计量进行校正,否则需对其进行校正;

并想知道哪两组间有统计学意义,

可进行两两比较.



以3组独立样本的定量数据为例,

来实现上述假设检验.



data score;
  do group = 1 to 3;  *group=分组;
   do id = 1 to 18;     *各组样本均为18;
input score @@;     *score=得分;
output;
end;
  end;
  cards;
0 17 17 17 18 20 22 23 23 23 24 24 24 25 28 29 30 30
22 22 23 23 24 25 25 26 26 28 28 29 29 29 29 29 30 30
23 25 25 26 27 28 28 29 29 29 29 30 30 31 31 33 35 37
;
run;



proc univariate data=score normal; *正态性检验;* 正态性不是检验得到的,或许要结合专业背景;
var score;
class group;

run;



proc anova data=score;
class group;
model score = group;
means group/hovtest = levene;  *方差齐性检验;
   run;



proc npar1way data=score;  *Kruskal-Wallis检验 结果显示 3组的总体中位数不全相等, 接着可进行两两比较;
var score;
class group;
run;



proc rank data=score out=rankscore; *生成秩次;
var score;
ranks rs;
run;



proc anova data=rankscore;
class group;
model rs=group;
means group/snk;   *对各组秩次进行两两比较;
run;

priss111 发表于 2012-5-7 11:58
对单因素2组或者多组定量资料进行假设检验时,

首先想到的也许是ANOVA,

谢谢您的详细指导。

thanks for this. good to learn.

webgu 发表于 2012-5-7 09:22

请问这个适用于三组样本量不均吗