data(eduperform)
model1 <- lm(pi ~ hs, data = eduperform)
model2 <- lm(pi ~ hs + ps, data = eduperform)
model3 <- lm(pi ~ hs + ps + sh, data = eduperform)
models <- list(model1, model2, model3)
select_optimal_model(models, criterion= "AIC")
Run the code above in your browser using DataLab