LinearRegression

# Read Data
vehicle <- read.csv("https://raw.githubusercontent.com/bkrai/Statistical-Modeling-and-Graphs-with-R/main/vehicle.csv", header = TRUE)

# Data Review and Data Preparation
vehicle$lh[vehicle$lh==0] <- mean(vehicle$lh)
vehicle$lc[vehicle$lc==0] <- mean(vehicle$lc)

# Data Partition
set.seed(1234)
ind <- sample(2, nrow(vehicle), 
              replace = TRUE, 
              prob = c(0.7, 0.3))
training <- vehicle[ind==1,]
testing <- vehicle[ind==2,]

# Multiple Linear Regression
model <- lm(lc~lh, data=training)
model
summary(model)
plot(lc~lh, training)
abline(model, col = "blue")

# Model Diagnostics
par(mfrow=c(2,2))
plot(model)
vehicle[1620,]

# Prediction
pred <- predict(model, testing)
predict(model, data.frame(lh=10))