2014-02-13 53 views
1

背景:我有一个分类变量,X,与我适合作为单独的虚拟变量四个级别。因此,有三个虚拟变量代表x = 1,x = 2,x = 3(x = 0是基线)。SAS对比:“权值计算”在PROC GLM参数估计的线性组合

问题:我想测试模型参数的线性组合的显着性,例如:2 * B1 + 2 * B2 + B3 = 0。

在Stata,第一个问题可以很容易的模型后做的是配合使用以下:

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

现在,如果我想使用CONTRAST语句PROC GLM为此在SAS,我例如,如果在一个不相关的例子中,我想测试以下四个连续变量: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声明的功能,专门为虚拟变量编码的分类变量,如问题中最初所述。

回答

2

使用PROC REG。

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

谢谢。我很难跟踪SAS各种回归程序支持的内容。我没有想到PROC REG,因为我之前曾尝试过一个简单的*计算*(而不是测试)线性组合参数(即将SAS用作“计算器”和存储的beta估计值),就像Stata的' lincom'命令和PROC REG不支持。这实际上也是我需要做的事情,我会将它添加到我的问题中。 – Meg

+0

没关系 - 我认为它更适合作为自己的问题。现在会发布新的问题。 – Meg

+1

附加的,有点相关的问题:http://stackoverflow.com/questions/21758486/sas-automatically-calculate-linear-combination-of-parameter-estimates-with-pr – Meg