# Load caret librarylibrary(caret)# Load Boston Housing Datasetboston = read.csv("boston.csv")# Use 10-fold X-validationcrtl = trainControl(method="cv",number=10)# Model 1: Plain Vanilla Regressionmodel1 = train(MEDV~., data=boston, method="lm", trControl=ctrl)# Show resultsmodel1