假设我想进行多元多元回归分析 并测试(使用单个测试)假设,即两个解释变量的回归参数为0.
在R
中,我会执行以下操作:
(假设数据是一个data.frame,列名为x1 x2(对于DV)和x3 x4(对于IV)
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。
答案 0 :(得分:1)
试试这个:
PROC GLM DATA=thedata;
MODEL x1 x2 = x3 x4 /SS3 NOUNI;
MANOVA H=_ALL_/PRINTE PRINTH;
RUN;
答案 1 :(得分: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;