Я хотел бы использовать recipes для построения модельной формулы и терминов, а затем вычислять предельные эффекты. Эти термины могут включать взаимодействия и полиномы.
Например,
library(recipes)
rec <- recipe(mtcars, mpg ~ cyl + hp) %>%
step_num2factor(cyl) %>%
step_dummy(cyl) %>%
step_poly(hp) %>%
step_interact(terms = ~ starts_with("cyl"):starts_with("hp")) %>%
step_normalize(all_predictors(), all_outcomes()) %>%
prep()
fit <- lm(formula = formula(rec), data = juice(rec))
Можно ли получить предельные эффекты для hp из этой подгонки модели? Я знаю, как рассчитать предельные эффекты, если взаимодействия и полиномиальные термины созданы внутри lm, но мне интересно, есть ли способ создать формулу модели полностью с recipes, а затем рассчитать предельные эффекты.