SAS中文论坛

 找回密码
 立即注册

扫一扫,访问微社区

查看: 756|回复: 8
打印 上一主题 下一主题

平均数t检定的问题?

[复制链接]

49

主题

76

帖子

1462

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
1462
楼主
 楼主| 发表于 2011-2-19 11:23:41 | 只看该作者

平均数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
回复 支持 反对

使用道具 举报

49

主题

76

帖子

1462

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
1462
沙发
 楼主| 发表于 2011-2-19 23:12:45 | 只看该作者

Re: 平均数t检定的问题?

难道不是pairwise comparison里不做MULTIPLICITY ADJUSTMENT的Student’s for Individual? 比如:
[code:1f27gqjq]PROC GLM;
...
LSMEANS/PDIFF =control ('1') ADJUST=T ;
QUIT;[/code:1f27gqjq]
回复 支持 反对

使用道具 举报

49

主题

76

帖子

1462

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
1462
板凳
 楼主| 发表于 2011-2-20 11:27:55 | 只看该作者

Re: 平均数t检定的问题?

但我只要想要比较各组之间的平均值t检定而已,看看是有具差异性,那应怎么做呢??
回复 支持 反对

使用道具 举报

49

主题

76

帖子

1462

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
1462
地板
 楼主| 发表于 2011-2-20 12:04:58 | 只看该作者

Re: 平均数t检定的问题?

SAS 有个语句proc ttest
proc anova,或glm也可以的。
回复 支持 反对

使用道具 举报

49

主题

76

帖子

1462

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
1462
5#
 楼主| 发表于 2011-2-20 16:07:56 | 只看该作者

Re: 平均数t检定的问题?

需要将数据集格式转换一下,再用,proc ttest的paired语句
回复 支持 反对

使用道具 举报

49

主题

76

帖子

1462

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
1462
6#
 楼主| 发表于 2011-2-21 01:39:20 | 只看该作者

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。即你对下一个问题的解答。
个人理解,探讨为主。京剧
回复 支持 反对

使用道具 举报

49

主题

76

帖子

1462

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
1462
7#
 楼主| 发表于 2011-2-21 13:12:22 | 只看该作者

Re: 平均数t检定的问题?

@jingju11
呵呵,对这个题按LZ每组对比的意思,我的想法是根据group将原数据集transpose成return1...return5,5列,再用如下一题的paired的方法。
回复 支持 反对

使用道具 举报

49

主题

76

帖子

1462

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
1462
8#
 楼主| 发表于 2011-2-21 14:07:15 | 只看该作者

Re: 平均数t检定的问题?

请你讲一讲用paired t test的理由。对于我们来说,即使统计不是精准的科学,你的做法也必应有你的理由。你何以猜得Benny是要dependent t test而不是independent one?
回复 支持 反对

使用道具 举报

49

主题

76

帖子

1462

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
1462
9#
 楼主| 发表于 2011-2-21 17:40:28 | 只看该作者

Re: 平均数t检定的问题?

恩,前辈教训的是。之前小弟草率主观判断。做任何统计处理、分析都要有充分的理由,不能为了处理数据而处理数据。
感谢指正!
[quote="jingju11":1cwjyr5o]请你讲一讲用paired t test的理由。对于我们来说,即使统计不是精准的科学,你的做法也必应有你的理由。你何以猜得Benny是要dependent t test而不是independent one?[/quote:1cwjyr5o]
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|小黑屋|手机版|Archiver|SAS中文论坛  

GMT+8, 2025-6-12 18:43 , Processed in 0.070099 second(s), 21 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.

快速回复 返回顶部 返回列表