求助。如题。我的数据及程序如下。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比较的结果不一致?
先谢谢大侠了!