作者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