一个简单的问题 R中什么命令同时检验a1=a2=0呢？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

比如y=a1*x1+a2*x2+a3*x3+a4*x4+u
在R中用什么命令实现检验a1=a2=0呢？高手请不吝赐教！
先表示感谢了！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：检验

做拟合吗？拟合的话在括号里面加fixed（NA，NA，0,0,0）

lenovoq 发表于 2012-6-6 23:31
做拟合吗？拟合的话在括号里面加fixed（NA，NA，0,0,0）

做检验！

你可以去做Q(m)检验。不过你拟合出来的结果里面已经有p值了，若p>显著水平，则不拒绝原假设，不显著，即可以视为0。然后需要重新拟合模型，就是用2楼的命令。

lm=lm(y ~ ., data = freeny)
Rmat=rbind(c(0,1,0,0,0),c(0,0,1,0,0))
rvec = c(0,0)
bhat = coef(lm)
avarRbhat = Rmat%*%vcov(lm)%*%t(Rmat)
wald.stat = t(Rmat%*%bhat-rvec)%*%solve(avarRbhat)%*%(Rmat%*%bhat-rvec)
as.numeric(wald.stat)             #  41.18802
p.wald = 1 - pchisq(wald.stat,2)  #  1.137979e-09
p.wald

##########wald test in e-views
Wald Test:                        
Equation: EQ01                        
                        
Test Statistic        Value                      df            Probability
                        
F-statistic        20.59401             (2, 34)          0.0000
Chi-square        41.18802                   2          0.0000