作者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/cn.aspx?n=bbs/Statistics/M.1415049307.A.C94.html