如何比较SAS中的均值(µ1 +µ2 +µ3)/3 = (µ3 +µ4)/2 :使用'ESTIMATE‘还是'CONTRAST'？

Question

我正在使用5种装配部件的技术的数据集，从每种技术中随机抽取工人样本，以评估工人需要多长时间才能完成任务。我希望将means与t-test进行比较，但我正在努力获得正确的代码，因为我对SAS非常陌生。

可以使用以下代码找到该数据集：

data Ex1;
input technique time @@;
lines;
1 45.6
1 41
1 46.4
1 50.7
1 47.9
1 44.6
2 41
2 49.1
2 49.2
2 54.8
2 45
3 51.7
3 60.1
3 52.6
3 58.6
3 59.8
3 52.6
3 53.8
4 67.5
4 57.7
4 58.2
4 60.6
4 57.3
4 58.3
4 54.8
5 57.1
5 69.6
5 62.7
;
run;

我希望使用PROC GLM来测试(µ1 +µ2 +µ3)/3 = (µ3 +µ4)/2的零假设，而不是这些均值不相等的替代方案。我有以下用于此操作的代码，但在运行它时遇到错误：

proc glm data=Ex1;
class technique;
model time=technique/NOINT SOLUTION E;

CONTRAST 'M1+M2+M3=M3+M4' technique 1 1 1 0 0/DIVISOR=3, technique 0 0 1 1 0/DIVISOR=3;
run;

会产生以下输出错误：

1          OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK;
 72         
 73         proc glm data=Ex1;
 74         class technique;
 75         model time=technique/NOINT SOLUTION E;
 76         
 77         CONTRAST 'M1+M2+M3=M3+M4' technique 1 1 1 0 0/DIVISOR=3, technique 0 0 1 1 0/DIVISOR=3;
                                                          _______
                                                          22
                                                          76
 NOTE: The previous statement has been deleted.
 ERROR 22-322: Syntax error, expecting one of the following: ;, E, EST, ETYPE, SINGULAR.  
 ERROR 76-322: Syntax error, statement will be ignored.
 78         run;
 
 NOTE: Due to the presence of CLASS variables, an intercept is implicitly fitted.  R-Square has been corrected for the mean.
 79         
 80         OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK;
 92 

这是ESTIMATE问题还是CONTRAST问题？

Answer 1

在GLM程序中没有用于对比度的选项除数。

让我们简化你的初始方程。

​

​

完整进程GLM调用：

proc glm data=Ex1;
  class technique;
  model time=technique/NOINT SOLUTION E;

  ESTIMATE '(M1+M2+M3)/3-(M3+M4)/2' technique 2 2 -1 -3 0 / divisor=6;
  CONTRAST '(M1+M2+M3)/3=(M3+M4)/2' technique 0.33333 0.333333 -0.166667 -0.5 0 / e ;
run;

和输出：

​

​