SAS CONTRAST:PROC GLM中参数估计的线性组合的“加权”

时间:2014-02-13 02:50:54

标签: testing sas contrast categorical-data

背景:我有一个分类变量 X ,有四个级别,我适合作为单独的虚拟变量。因此,存在三个总虚拟变量,表示x = 1,x = 2,x = 3(x = 0是基线)。

问题/问题:我想测试模型参数的线性组合的重要性,例如:2 * B1 + 2 * B2 + B3 = 0.

在Stata中,使用以下内容完成模型后,可以轻松完成第一个问题:

test 2*B1 + 2*B2 + B3 = 0

现在,如果我想使用PROC GLM语句在SAS中为CONTRAST执行此操作,我知道我的“权重”(缺少更好的术语)必须总和为0.例如,如果,在一个不相关的例子中,我想测试以下四个连续变量:C1 + C2 = C3 + C4,我的对比语句如下:

CONTRAST 'Contrast1' C1 0.5 C2 0.5 C3 -0.5 C4 -0.5

在这种情况下,很明显每个变量应该如何加权。但是,当我想将上面模型中给出的系数(2 * B1 + 2 * B2 + B3 = 0)与这些权重相结合时,我不清楚如何在CONTRAST语句中对函数进行加权,特别是对于虚拟的可变编码分类变量,如问题中最初所述。

1 个答案:

答案 0 :(得分:2)

使用PROC REG。

proc reg data=mydata;
model y = b1 b2 b3;
test 2*b1+2*b2+b3=0;
run;
quit;