全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 SAS专版
3759 3
2011-04-18
各位路过的前辈们,小的用别人的数据做了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;
二维码

扫码加我 拉你入群

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

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

全部回复
2011-4-19 01:20:37
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.
二维码

扫码加我 拉你入群

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

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

2011-4-19 12:30:14
你好,谢谢你的回复,我试看看换个参数,
试验是这个的,研究某一降压药,随机分两组,观察0,2,4,6,8周收缩压的变化。
二维码

扫码加我 拉你入群

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

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

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

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群