如何在SAS中比较均值（µ1 + µ2 + µ3）/ 3 =（µ3 + µ4）/ 2：使用'ESTIMATE'或'CONTRAST'？

ihb 发表于 Dev

我正在使用5种组装零件技术的数据集，并从每种技术中随机抽取了工人样本，以评估工人完成任务需要多长时间。我想将均值与进行比较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问题吗？

乌鸦

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

让我们简化您的初始方程式。

完整的PROC 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;

并输出：