如何通过配方获得多项式项的边际效应？

Question

我想使用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完全创建模型公式，然后计算边际效应