if (requireNamespace("MASS", quietly = TRUE)) {
data(Boston, package = "MASS")
# Simple linear regression
model_simple <- reg_lm(formula = medv ~ lstat)
model_simple <- fit(model_simple, Boston)
pred_simple <- predict(model_simple, Boston)
head(pred_simple)
# Polynomial regression (degree 2)
model_poly <- reg_lm(formula = medv ~ poly(lstat, 2, raw = TRUE))
model_poly <- fit(model_poly, Boston)
pred_poly <- predict(model_poly, Boston)
head(pred_poly)
# Multiple regression
model_multi <- reg_lm(formula = medv ~ lstat + rm + ptratio)
model_multi <- fit(model_multi, Boston)
pred_multi <- predict(model_multi, Boston)
head(pred_multi)
}
Run the code above in your browser using DataLab