x<-c(2,3,4,5,6)
y<-c(4,4,6,6,10)
data<-data.frame(x,y)
x;y

lm(y~x)
ols<-lm(y~x)
summary(ols)
confint(ols)

dx<-x-mean(x)
dy<-y-mean(y)
lm(dy~dx-1)
ols1<-lm(dy~dx-1)

library(stargazer)
stargazer(ols,ols1,type="text")

plot(y~x, xlim=c(0,7),ylim=c(0,10))
abline(lm(y~x))

library(predict3d)
ggPredict(ols,show.error=T)

library(tigerstats)
par(mfrow=c(2,1))
ptGC(c(-3.656,3.656),region="outside",df=3,graph=T)
ptGC(c(-3.182,3.182),region="outside",df=3,graph=T)

predict(lm(y~x))
new<-data.frame(x = 7)
predict(lm(y~x), new, se.fit = TRUE)
pred.w.plim <- predict(lm(y ~ x), new, interval = "prediction")
pred.w.plim
pred.w.clim<-predict(lm(y ~ x), new, interval = "confidence")
pred.w.clim