x2<-c(1,2,3,2) x3<-c(2,1,1,2) y<-c(1,1,2,3) ols<-lm(y~x2+x3) summary(ols) confint(ols) resid(ols) predict(lm(y~x2+x3)) new<-data.frame(x2=3, x3=2) predict(lm(y~x2+x3), new, se.fit = TRUE) pred.w.clim<-predict(lm(y~x2+x3), new, interval = "confidence") pred.w.clim pred.w.plim<-predict(lm(y~x2+x3), new, interval = "prediction") pred.w.plim library(car) library(stargazer) x_2<-c(1,2,3,2,4) x_3<-c(2,1,1,2,3) x_4<-c(2,3,2,4,2) y_n<-c(1,1,2,3,3) m<-length(y_n) k=4 ur.lm<-lm(y_n~x_2+x_3+x_4) summary(ur.lm) ur.r2<-summary(ur.lm)$r.squared ur.r2 r.lm<-lm(y_n~x_2) r.r2<-summary(r.lm)$r.squared r.r2 stargazer(ur.lm, r.lm, type="text") Fstat1<-((ur.r2-r.r2)/2)/((1-ur.r2)/(m-k)) Fstat1 pval<-1-pf(Fstat1, 2, m-k) pval jointHo<-c("x_3","x_4") linearHypothesis(ur.lm, jointHo)