{base}
{stats}
{tidymodels}
parsnip
{dials}
workflow_set
stacks
bean_rec_trained %>% step_pls(all_numeric_predictors(), outcome = "class", num_comp = 4) %>% plot_validation_results() + ggtitle("Partial Least Squares")