attach(airquality) res.3 <- lm(Ozone ~ Solar.R+Wind+Temp) res.2 <- lm(Ozone ~ Solar.R+Wind) lambda <- -2*(logLik(res.2)-logLik(res.3)) 1-pchisq(lambda,1) detach(airquality)