作者switchfish (生活=投稿+做计画+写论文)
看板Statistics
标题[问题] SAS跑 Mixed model 置入变项的问题
时间Sun Nov 7 12:27:26 2010
Dear 统计界的先进大家好
小弟最近利用SAS软体处理ㄧ个介入性的重复量测数据(实验组与对照组各有六次的检测值)
我想要看Treatment effect, Time effect, Treatment*Time(interaction)
在Mixed model中经常看到同一个模式放入 Group Time Group*Time
然後根据 Group Time Group*Time 三个P值来决定他们的Effect,请问这样判断是对的吗?
还是要类似跑线性回归时, 先在模式中个别单一放入 Group 或 Time
然後探讨interaction时, 另一个模式同时放入 Group Time Group*Time
以下大概是我用SAS建档的 Data step
DATA test;
INPUT NO. Group Time X(检测值);
CARDS;
1 I 1 22
1 I 2 32
1 I 3 23
1 I 4 45
1 I 5 31
1 I 6 45
1 C 1 32
1 C 2 22
1 C 3 21
1 C 4 24
1 C 5 33
1 C 6 36
2 I 1 42
2 I 2 42
2 I 3 33
2 I 4 55
2 I 5 41
2 I 6 42
2 C 1 30
2 C 2 32
2 C 3 41
2 C 4 42
2 C 5 36
2 C 6 38
RUN;
PROC MIXED DATA=test;
CLASS NO. Group Time;
MODEL X=Group Time Group*Time;
REPEATED Time / SUBJECT = NO. TYPE = cs RCORR;
RUN;
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 114.46.196.214