2014-03-25 132 views
1

假设我想执行一个多变量回归分析 和测试(使用单个测试)的假说,即回归参数对于两个解释变量是0多元在SAS多元回归

R我会做如下: (假设海图是列名X1 X2一个data.frame(为DVS)和X3 X4(对于IVS)

fit = manova(cbind(x1,x2) ~ cbind(x3,x4),data=thedata) 
summary(fit) 

这给了我Pillais'跟踪和大约关联的F 我问题是:什么程序会产生同样的分析在SAS? 我已经试过GLM的以下调用的变化:

PROC GLM DATA=thedata; 
MODEL x1 x2 = x3 x4 /SS3 NOUNI; 
MANOVA H=_ALL_; 
RUN; 

但我似乎无法获得超越测试X3和X4分别。

回答

1

我找到了答案,以我自己的问题用PROC REG

PROC REG DATA=thedata; 
MODEL x1 x2 = x3 x4; 
MTEST x3,x4; 
RUN; 

此外,在SAS论坛用户指出该解决方案对我说:

PROC GLM DATA=thedata; 
MODEL x1 x2 = x3 x4; 
CONTRAST 'All Zero' x3 1,x4 1; 
MANOVA; 
RUN; 
1

试试这个:

PROC GLM DATA=thedata; 
MODEL x1 x2 = x3 x4 /SS3 NOUNI; 
MANOVA H=_ALL_/PRINTE PRINTH; 
RUN; 
+0

SAS中的E矩阵会给我正确的错误矩阵,我原则上可以使用E中的E矩阵来计算Pillais轨迹通过运行GLM与MODEL x1 x2 =/SS3 NOUNI;因为这是我寻找的H矩阵。但我正在寻找一种方法让SAS执行此计算。 – surlol

+0

我不知道该怎么做;它应该是可能的,甚至可能已经完成,但我没有看到它。 –