I thought that the scores of the variables returned by step_select_linear with a binary outcome were the coefficients of the logistic regression model, but that is not the case.
library(recipes)
library(parsnip)
library(colino)
# load the example iris dataset
data(cells, package = "modeldata")
# create a preprocessing recipe
rec <-
recipe(class ~ ., data = cells[, -1]) %>%
step_select_linear(
all_predictors(),
outcome = "class",
threshold = 0.9
)
prepped <- prep(rec, new_data = cells[, -1])
prepped
step_1 <- prepped$steps[[1]]
step_1$scores %>% arrange(desc(score))
glm(class ~ inten_cooc_max_ch_3, data=cells, family=binomial )