### LESSON 5 EXERCISE OF SIMPLE REGRESSION ANALYSIS  (27/4/2021)

library(haven)
setwd("S:/_STAFF/Gulser/R NOTES/LESSON_5")
dataset<-read_dta('all752806.dta')
attach(dataset)
dataset
summary(dataset)

shapiro.test(statura)
hist(statura)

shapiro.test(peso)
hist(peso)

###Correlation

cor(peso,statura)
cor.test(peso,statura)

plot(peso,statura)

lm(peso~statura)

##Summary
summary(lm(peso~statura))

###Scatter Plot

plot(peso~statura)
abline(lm(peso~statura))

lm.vel <- lm(peso~statura)
lm.vel

predict(lm.vel)
predict(lm.vel,int="c")

resid(lm.vel)
resi<-resid(lm.vel)
resi

shapiro.test(resi)
plot(resi)