SAS中文论坛

标题: 平均数t检定的问题? [打印本页]

作者: shiyiming    时间: 2011-2-19 11:23
标题: 平均数t检定的问题?
请教各位前辈,我现在有一个关于平均数t检定的问题想请教大家,我的数据中可分四组(1、2、3、4、5),每组数量不一样。我现在想将其报酬两两做比较(也就是说1和2,1和3,1和4,1和5),看是否具有差异性,想请教各位前辈,若是要完成此项工作,其代码应该要如何设计呢,请各位前辈帮帮忙吧!
group         return
1        12
1        15
1        10
1        22
1        21
2        30
2        22
2        15
2        23
2        15
3        10
3        9
3        13
3        15
4        20
4        22
4        23
4        25
4        20
4        26
5              30
5              23
5              25
5              21
作者: shiyiming    时间: 2011-2-19 23:12
标题: Re: 平均数t检定的问题?
难道不是pairwise comparison里不做MULTIPLICITY ADJUSTMENT的Student’s for Individual? 比如:
[code:1f27gqjq]PROC GLM;
...
LSMEANS/PDIFF =control ('1') ADJUST=T ;
QUIT;[/code:1f27gqjq]
作者: shiyiming    时间: 2011-2-20 11:27
标题: Re: 平均数t检定的问题?
但我只要想要比较各组之间的平均值t检定而已,看看是有具差异性,那应怎么做呢??
作者: shiyiming    时间: 2011-2-20 12:04
标题: Re: 平均数t检定的问题?
SAS 有个语句proc ttest
proc anova,或glm也可以的。
作者: shiyiming    时间: 2011-2-20 16:07
标题: Re: 平均数t检定的问题?
需要将数据集格式转换一下,再用,proc ttest的paired语句
作者: shiyiming    时间: 2011-2-21 01:39
标题: Re: 平均数t检定的问题?
To benny10:
看样子你的意思类似如下:
[code:1pcrjs13]proc ttest data =had; where group in (1 2);
class group;
var return;
run;[/code:1pcrjs13]
然后做(1 3)(1 4)(1 5)。也就是说每次比较你只需要两个组别的信息而舍弃其他的组。如果如此,做4次t test即可。也可以改变data利用by。

[code:1pcrjs13]data have;
        set had end =Eof; by group notsorted;
        if group ^=1 then do;        ByClass =group; output;
                if last. group then do p =1 to nobs;
                        set had point =p nobs =nobs; if group =1 then output;
                end;
        end;
        if Eof then stop;
run;
proc ttest data =have; by ByClass notsorted;
        class group; var return;
run;[/code:1pcrjs13]
但是据我的理解两两t test经常意指如下:
[code:1pcrjs13]proc glm data =had;        class group;
        model return =group;
        lsmeans group/pdiff =control('1') tdiff adjust =t;
run; quit;[/code:1pcrjs13]
to HSSNOW:
Paired T test通常针对同组的test。也就是说差异为0的单组test。即你对下一个问题的解答。
个人理解,探讨为主。京剧
作者: shiyiming    时间: 2011-2-21 13:12
标题: Re: 平均数t检定的问题?
@jingju11
呵呵,对这个题按LZ每组对比的意思,我的想法是根据group将原数据集transpose成return1...return5,5列,再用如下一题的paired的方法。
作者: shiyiming    时间: 2011-2-21 14:07
标题: Re: 平均数t检定的问题?
请你讲一讲用paired t test的理由。对于我们来说,即使统计不是精准的科学,你的做法也必应有你的理由。你何以猜得Benny是要dependent t test而不是independent one?
作者: shiyiming    时间: 2011-2-21 17:40
标题: Re: 平均数t检定的问题?
恩,前辈教训的是。之前小弟草率主观判断。做任何统计处理、分析都要有充分的理由,不能为了处理数据而处理数据。
感谢指正!
[quote="jingju11":1cwjyr5o]请你讲一讲用paired t test的理由。对于我们来说,即使统计不是精准的科学,你的做法也必应有你的理由。你何以猜得Benny是要dependent t test而不是independent one?[/quote:1cwjyr5o]




欢迎光临 SAS中文论坛 (http://mysas.net/forum/) Powered by Discuz! X3.2