library(openxlsx)
library(stargazer)
library(car)

sample1<-read.xlsx("http://kanggc.iptime.org/data/pitcher.xlsx")

money<-sample1$MONEY
win<-sample1$WIN
four<-sample1$FOUR
save<-sample1$SAVE
lose<-sample1$LOSE
year<-sample1$YEAR

ols1<-lm(money~win+four+save+lose+year)
summary(ols1)

jointHo<-c("four","lose")
linearHypothesis(ols1, jointHo)

ols2<-lm(money~win+save+year)
summary(ols2)

stargazer(ols1, ols2, type="text")

new<-data.frame(win=14, save=5, year=10)
predict(lm(money~win+save+year), new, se.fit=T)

pred.w.clim<-predict(lm(money~win+save+year), new, interval="confidence")
pred.w.clim

pred.w.plim<-predict(lm(money~win+save+year), new, interval="prediction")
pred.w.plim