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