作者norchen (乌金彻令)
看板Statistics
标题Re: [问题] 关於重复量数的SAS语法
时间Thu May 15 18:12:52 2008
我把data 附上并更改所需的语法
资料档的部份版大若是习惯用CARDS指令
就省略INFILE的部份
※ 引述《norchen (乌金彻令)》之铭言：
: 我的实验设计是三因子重复量数
: 三个因子(SP SD PERS)都是within
: 每个因子有两个水准
: 只有一个依变数
三因子重复量数data
sub A1B1C1 A1B1C2 A1B2C1 A1B2C2 A2B1C1 A2B1C2 A2B2C1 A2B2C2
1 888.45 971.45 903.41 909.83 1087.34 1177 1031.95 1054.46
2 830.83 1027.18 1030.05 806.78 1009.63 1051.92 842 1259.05
3 2172.28 1095.28 1297.64 1723.57 1035.58 730.33 754.93 814.86
4 966.92 995.46 667.04 708.5 558.1 545.75 609.82 699
5 876.43 777.54 975.34 856.77 1099.45 1010.54 1020.43 1030.45
6 576.44 574.32 666.43 571.34 982.2 765.32 874.34 666.43
语法如下：
data A ;
INFILE 'C:\PRE-TEST.TXT' FIRSTOBS=2 EXPANDTABS;
INPUT SUB A1B1C1 A1B1C2 A1B2C1 A1B2C2 A2B1C1 A2B1C2 A2B2C1 A2B2C2;
PROC GLM data=A ;
CLASS sub;
MODEL A1B1C1 A1B2C1 A1B1C2 A1B2C2 A2B1C1 A2B2C1 A2B1C2 A2B2C2 = ;
REPEATED SP 2 CONTRAST(2) , SD 2 , PERS 2 /SUMMARY PRINTM PRINTE PRINTH ;
RUN;
QUIT;
=====
结果如下：
The GLM Procedure
Repeated Measures Analysis of Variance
Univariate Tests of Hypotheses for Within Subject Effects
Source DF Type III SS Mean Square F Value Pr > F
SP 1 27956.053 27956.053 0.10 0.7657
Error(SP) 5 1411990.737 282398.147
SD 1 18170.1919 18170.1919 0.57 0.4831
Error(SD) 5 158455.6816 31691.1363
PERS 1 22158.76963 22158.76963 3.64 0.1148
Error(PERS) 5 30452.04247 6090.40849
SP*SD 1 11180.69701 11180.69701 1.02 0.3590
Error(SP*SD) 5 54838.84834 10967.76967
SP*PERS 1 1204.40403 1204.40403 0.06 0.8140
Error(SP*PERS) 5 97918.40902 19583.68180
SD*PERS 1 66693.9210 66693.9210 0.89 0.3894
Error(SD*PERS) 5 375673.7127 75134.7425
SP*SD*PERS 1 12.7927 12.7927 0.00 0.9877
Error(SP*SD*PERS) 5 244726.6323 48945.3265
========
我另外使用单变量的方式跑
DATA如下：
SUB RT SP SD PERS
1 888.45 1 1 1
1 971.45 1 2 1
1 903.41 1 1 2
1 909.83 1 2 2
1 1087.34 2 1 1
1 1177 2 2 1
1 1031.95 2 1 2
1 1054.46 2 2 2
2 830.83 1 1 1
2 1027.18 1 2 1
2 1030.05 1 1 2
2 806.78 1 2 2
2 1009.63 2 1 1
2 1051.92 2 2 1
2 842 2 1 2
2 1259.05 2 2 2
3 2172.28 1 1 1
3 1095.28 1 2 1
3 1297.64 1 1 2
3 1723.57 1 2 2
3 1035.58 2 1 1
3 730.33 2 2 1
3 754.93 2 1 2
3 814.86 2 2 2
4 966.92 1 1 1
4 995.46 1 2 1
4 667.04 1 1 2
4 708.5 1 2 2
4 558.1 2 1 1
4 545.75 2 2 1
4 609.82 2 1 2
4 699 2 2 2
5 876.43 1 1 1
5 777.54 1 2 1
5 975.34 1 1 2
5 856.77 1 2 2
5 1099.45 2 1 1
5 1010.54 2 2 1
5 1020.43 2 1 2
5 1030.45 2 2 2
6 576.44 1 1 1
6 574.32 1 2 1
6 666.43 1 1 2
6 571.34 1 2 2
6 982.2 2 1 1
6 765.32 2 2 1
6 874.34 2 1 2
6 666.43 2 2 2
语法如下：
data A ;
INFILE 'C:\PRE-TEST2.TXT' FIRSTOBS=2 EXPANDTABS;
INPUT SUB RT SP SD PERS;
PROC GLM data=A ;
CLASS SUB SP SD PERS;
MODEL RT = SP SD PERS SUB*SP SUB*SD SUB*PERS SUB*SP*SD SUB*SP*PERS
SUB*SD*PERS SUB*SP*SD*PERS SP*SD SP*PERS SD*PERS SP*SD*PERS ;
LSMEANS SP*SD*PERS/E= SUB*SP*SD*PERS TDIFF PDIFF STDERR;
TEST H= SP*SD SP*PERS SD*PERS SP*SD*PERS E= SUB*SP*SD*PERS;
RUN;
QUIT;
结果如下：
The GLM Procedure
Dependent Variable: RT
Tests of Hypotheses Using the Type III MS for SUB*SP*SD*PERS as an Error Term
Source DF Type III SS Mean Square F Value Pr > F
SP*SD 1 11180.69701 11180.69701 0.23 0.6529
SP*PERS 1 1204.40403 1204.40403 0.02 0.8815
SD*PERS 1 66693.92101 66693.92101 1.36 0.2957
SP*SD*PERS 1 12.79268 12.79268 0.00 0.9877
========
两种方式的结果不一样
或许是因为用了不同的ERROR TERM
我想知道的是
我应该使用那一种分析方式
以及事後比较的部份为何跑不出来
(我自己认为是第一种，但是就是无法跑出事後比较~.~")
========
若版上的先进认为这需要有偿教学
请来信连系
谢谢
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 218.174.37.143