关于proc mixed的应用

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

各位路过的前辈们，小的用别人的数据做了proc mixed
但是结论跟他的结论居然是相反的。不知道小的程序是否有问题？
附上数据与程序，求帮忙！
data treat;
input pnt treat sbp0 sbp2 sbp4 sbp6 sbp8@@;
cards;
1 1 20.2 18.0 16.6 16.1 16.0
2 1 19.2 18.0 17.3 17.4 16.6
3 1 20.1 20.0 18.6 19.0 17.6
4 1 18.2 18.0 17.3 17.0 15.8
5 1 19.7 16.4 17.4 16.8 16.9
6 1 21.5 17.3 16.1 . .
7 1 19.7 18.6 17.3 17.3 17.3
8 1 20.8 19.9 18.5 18.0 17.0
9 1 18.9 17.4 17.2 16.6 15.6
10 1 20.6 19.5 19.0 19.0 18.8
11 0 20.8 19.8 18.7 18.1 18.0
12 0 19.7 18.7 17.3 17.7 17.6
13 0 22.1 20.0 18.5 19.1 17.5
14 0 21.1 19.9 18.0 17.3 17.1
15 0 20.7 18.4 18.8 17.9 17.1
16 0 23.5 22.5 23.0 21.5 20.0
17 0 22.5 21.0 21.0 20.0 20.0
18 0 24.0 23.5 23.0 22.5 22.5
19 0 21.5 21.0 20.0 20.0 18.5
20 0 24.6 22.0 21.0 20.5 20.6
;
data b;
set treat;
array t{4} sbp2 sbp4 sbp6 sbp8;
do i=1 to 4;
if i=1 then time=2;
if i=2 then time=4;
if i=3 then time=6;
if i=4 then time=8;
y=t{i};
timepnt=time;
output;
end;
drop i sbp2 sbp4 sbp6 sbp8;
run;
proc mixed data=b covtest method=ml;
class treat pnt timepnt;
model y=treat time /htype=3s;
repeated timepnt/type=un subject=pnt r;
run;

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：proc mixed Mixed mix ROC Repeated treat 程序

I do not see mistakes from your code. But you can change your options for METHOD and TYPE. Also, what happened to you sbp0 ? I cannot type Chinese.

你好，谢谢你的回复，我试看看换个参数，
试验是这个的，研究某一降压药，随机分两组，观察0，2，4，6，8周收缩压的变化。

feifeiw09
我找到跟例子的区别了，他多了以下几步
y1=sbp1-sbp0;
y2=sbp2-sbp0;
y3=sbp3-sbp0;
y4=sbp4-sbp0;
做混合线性模型的时候治疗后测量值需要扣掉治疗前的吗？
不可以直接比较后4次测量结果的差异？