Linear model with 4 variables in R

Linear_model.R

#-------------------------------------------------------------------------------
#  Another linear model using more information y = a*x1 + b*x2 + c*x3 + d*x4
#-------------------------------------------------------------------------------

# Plot data
plot(disp,wt,type='p',xlab='Disp',ylab='Wt',main='Linear regression')

# Add a legend
legend("topleft",c("Observ.","Predicted"), col=c("black","red"), lwd=3)

model4 <- lm(wt ~ disp+mpg+hp+drat)
summary(model4)
coef(model4)

predicted <- fitted(model4)
original <- wt
mat <- cbind(predicted,original)

points(disp,predicted, type="p", col="red")