作者open2356 (Shing)
看板Statistics
標題[問題] 請教有關SAS 重複測量中的多重比較
時間Tue Nov 4 05:15:05 2014
各位前輩好，小弟剛入研究所接觸SAS沒多久，目前有一些疑問想請教一下..
我的試驗是不同濃度(B 0 5 10 15 20)的試劑放進產品，然後了解其pH值的變化，
那試驗中要測量不同時間點(0 2 4 6 8週)的data
目前小弟的程度判斷應該是用重複測量的方式下去跑數據，
系上習慣使用duncan來做多重比較，用混合模式沒辦法用means，
而lsmeans也沒辦法使用duncan校正，
比較麻煩的是使用tukey法，試驗結果幾乎沒有差異，用LSD比較又太丟臉，
這些狀況導致論文還沒辦法寫QQ
現在弟寫出這些語法，但有遇到一個瓶頸。
───────────────────────────────────
data abc;
input trt$ week1 week2 week3 week4 week5;
datalines;
B 6.275 6.370 6.390 6.415 6.565
B 5.945 5.775 6.215 6.135 6.130
B 5.970 5.775 6.195 6.120 6.130
0 6.305 6.410 6.410 6.435 6.585
0 6.405 6.425 6.665 6.580 6.570
0 6.535 5.430 6.750 6.690 6.570
5 6.350 6.445 6.445 6.470 6.625
5 6.455 6.525 6.730 6.625 6.635
5 6.440 6.525 6.725 6.650 6.635
10 6.375 6.490 6.510 6.565 6.660
10 6.395 6.475 6.670 6.590 6.590
10 6.390 6.475 6.640 6.580 6.590
15 6.360 6.435 6.485 6.520 6.660
15 6.380 6.465 6.665 6.570 6.580
15 6.405 6.465 6.680 6.600 6.580
20 6.345 6.425 6.445 6.475 6.630
20 6.420 6.470 6.670 6.595 6.590
20 6.365 6.470 6.680 6.620 6.590
;
proc print data=abc;
var trt week1 week2 week3 week4 week5;
proc glm data=abc;
class trt;
model week1-week5 =trt / nouni;
repeated time 5 / summary printe;
means trt / duncan;
────────────────────────────────────
目前跑出來的結果可以知道在不同時間點(如week1，也就是第0週)中各組間的差異，
不過我也想要知道同一個處理組(如B組)在不同時間點的差異顯著性，
但現在太淺，也爬了好陣子文仍然不知道怎麼處理，
想請教各位前輩們，不知道是該增加特別的語法還是分組的方式來解決，
還是有更好的比較法能解決這個問題，小弟真的不知道該怎麼辦才好，
在這裡先向大家跪謝了QQ
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.127.12.175
※ 文章網址: http://webptt.com/m.aspx?n=bbs/Statistics/M.1415049307.A.C94.html