全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 SAS专版
4314 0
2015-04-26
求助。如题。我的数据及程序如下。time1 分三次测,1、2、3;测量值为count,是个2分类变量。
data a;
/*time 1-pre 2-MMT 3-during count 1-no 0-yes */
input id time1 count @@;
cards;
1001        1        0
1001        3        0
1001        2        0
1002        1        0
1002        3        0
1002        2        0
1003        1        0
1003        3        0
1003        2        0
1004        1        0
1004        3        0
1004        2        0
1005        1        0
1005        3        0
1005        2        0
1006        1        0
1006        3        0
1006        2        0
1007        1        0
1007        3        0
1007        2        0
1008        1        0
1008        3        0
1008        2        0
1009        1        1
1009        3        0
1009        2        0
1010        1        0
1010        3        0
1010        2        0
1011        1        0
1011        3        0
1011        2        0
1012        1        1
1012        3        0
1012        2        0
1013        1        0
1013        3        0
1013        2        0
1014        1        0
1014        3        0
1014        2        0
1015        1        1
1015        3        0
1015        2        1
1016        1        0
1016        3        0
1016        2        0
1017        1        1
1017        3        1
1017        2        1
1018        1        0
1018        3        0
1018        2        0
1019        1        0
1019        3        0
1019        2        0
1020        1        0
1020        3        0
1020        2        0
1021        1        0
1021        3        0
1021        2        0
1022        1        1
1022        3        0
1022        2        0
1023        1        0
1023        3        0
1023        2        1
1024        1        0
1024        3        0
1024        2        0
1025        1        1
1025        3        1
1025        2        1
1026        1        1
1026        3        0
1026        2        0
1027        1        0
1027        3        1
1027        2        0
1028        1        1
1028        3        1
1028        2        1
1029        1        0
1029        3        0
1029        2        0
1030        1        0
1030        3        0
1030        2        0
1031        1        1
1031        3        0
1031        2        0
1032        1        0
1032        3        0
1032        2        0
1033        1        0
1033        3        0
1033        2        0
1034        1        1
1034        3        0
1034        2        0
1035        1        0
1035        3        0
1035        2        0
1036        1        1
1036        3        0
1036        2        0
1037        1        0
1037        3        0
1037        2        0
1038        1        0
1038        3        0
1038        2        0
1039        1        0
1039        3        0
1039        2        0
1040        1        0
1040        3        0
1040        2        0
1041        1        0
1041        3        0
1041        2        0
1042        1        0
1042        3        0
1042        2        0
1043        1        0
1043        3        0
1043        2        0
1044        1        0
1044        3        0
1044        2        0
1045        1        0
1045        3        0
1045        2        0
1046        1        0
1046        3        0
1046        2        0
1047        1        0
1047        3        0
1047        2        0
1048        1        0
1048        3        0
1048        2        0
1049        1        0
1049        3        0
1049        2        0
1050        1        0
1050        3        0
1050        2        0
1051        1        0
1051        3        0
1051        2        0
1052        1        0
1052        3        0
1052        2        0
1053        1        0
1053        3        0
1053        2        0
1054        1        1
1054        3        0
1054        2        0
1055        1        0
1055        3        0
1055        2        0
1056        1        0
1056        3        0
1056        2        0
1057        1        0
1057        3        0
1057        2        0
1058        1        0
1058        3        0
1058        2        0
1059        1        0
1059        3        0
1059        2        0
1060        1        0
1060        3        0
1060        2        0
1061        1        0
1061        3        0
1061        2        0
1062        1        0
1062        3        0
1062        2        0
1063        1        0
1063        3        0
1063        2        0
1064        1        0
1064        3        0
1064        2        0
1065        1        0
1065        3        0
1065        2        0
1066        1        0
1066        3        0
1066        2        0
1067        1        0
1067        3        0
1067        2        0
1068        1        0
1068        3        0
1068        2        0
1069        1        0
1069        3        0
1069        2        0
1070        1        1
1070        3        1
1070        2        1
1071        1        1
1071        3        0
1071        2        0
1072        1        0
1072        3        0
1072        2        0
1073        1        0
1073        3        0
1073        2        0
1074        1        0
1074        3        0
1074        2        0
1075        1        0
1075        3        0
1075        2        0
1076        1        0
1076        3        0
1076        2        0
1077        1        0
1077        3        0
1077        2        0
1078        1        0
1078        3        0
1078        2        0
1079        1        0
1079        3        0
1079        2        0
1080        1        0
1080        3        0
1080        2        0
1081        1        0
1081        3        0
1081        2        0
1082        1        0
1082        3        0
1082        2        0
1083        1        0
1083        3        0
1083        2        0
1084        1        0
1084        3        0
1084        2        0
1085        1        0
1085        3        0
1085        2        0
1086        1        0
1086        3        0
1086        2        0
1087        1        0
1087        3        0
1087        2        0
1088        1        0
1088        3        0
1088        2        0
1089        1        0
1089        3        0
1089        2        0
1090        1        0
1090        3        0
1090        2        0
1091        1        0
1091        3        0
1091        2        0
1092        1        0
1092        3        0
1092        2        0
1093        1        0
1093        3        0
1093        2        0
1094        1        0
1094        3        0
1094        2        0
1095        1        0
1095        3        0
1095        2        0
1096        1        0
1096        3        0
1096        2        0
1097        1        0
1097        3        0
1097        2        0
1098        1        0
1098        3        0
1098        2        0
1099        1        0
1099        3        0
1099        2        0
1100        1        0
1100        3        0
1100        2        0
;
run;

proc mixed method=reml;
class  time1;
model count= id time1 ;
repeated time1/type=un subject=id r;
contrast 'time1' time1 1 -1 0;
contrast 'time1' time1 1 0 -1;
contrast 'time1' time1 0 1 -1;
lsmeans time1;
run;

问 contrast与lsmeans要比较的是否一致。都是比较的什么?为何contrast与lsmeans比较的结果不一致?
先谢谢大侠了!

二维码

扫码加我 拉你入群

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

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

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

说点什么

分享

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