有关GLM里contrast语句的应用

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

请教各位，GLM里contrast语句的写法，
比如 想比较 A与B很简单  contrast  'A VS B' time 1 -1;
但是如果再复杂一些.比如我有4个时间点，每个时间点分为4组，我想在每个时间都进行两两比较，
用GLM里的contrast 如何写？？

参考数据：
data muscles;
do Time=1 to 4;
do Current=1 to 4;
input MuscleWeight @@;
output;
end;
end;
datalines;
72 74 69 61 61 65 62 65 70 85 76 61
67 52 62 60 55 59 64 65 64 67 72 60
57 66 72 72 43 43 63 66 72 56 75 92
57 56 78 60 63 58 61 79 68 73 86 71
46 74 58 60 64 52 71 64 71 53 65 66
44 58 54 57 55 51 62 61 79 60 78 82
53 50 61 56 57 56 56 56 71 56 58 69
46 55 64 56 55 57 64 66 62 59 58 88
;run;

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：CONTRAST contr RAS con CURRENT contrast

没人回答吗？顶起来！！！

data muscles;
do Time=1 to 4;
do Current=1 to 4;
input MuscleWeight @@;
output;
end;
end;
datalines;
72 74 69 61 61 65 62 65 70 85 76 61
67 52 62 60 55 59 64 65 64 67 72 60
57 66 72 72 43 43 63 66 72 56 75 92
57 56 78 60 63 58 61 79 68 73 86 71
46 74 58 60 64 52 71 64 71 53 65 66
44 58 54 57 55 51 62 61 79 60 78 82
53 50 61 56 57 56 56 56 71 56 58 69
46 55 64 56 55 57 64 66 62 59 58 88
;run;
/*我理解的意思是在每个时间点，对不同的current进行两两比较*/
proc glm;
class time current;
model muscleweight=time|current;
contrast 'when time=1  current=1 vs current=2'
current 1 -1 0 0
time*current
1 -1 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
contrast 'when time=1  current=1 vs current=3'
current 1 0 -1 0
time*current
1 0 -1 0 0 0 0 0 0 0 0 0 0 0 0 0;
run;

/*这么复杂的contrast我也没写过，不敢保证一定正确，可以@jingju大侠*/

感谢ziyenano的回复！
有谁能哪位大侠能介绍一下contrast书写的技巧！

http://www.doc88.com/p-893576945977.html可以看看这个，我就是看的这个。

同问，有没有朋友有更详尽从解答呀？谢谢