R recipes step_spline_nonnegative 非负样条

step_spline_nonnegative() 创建创建非负样条线特征的配方步骤的规范。

用法

step_spline_nonnegative(
  recipe,
  ...,
  role = "predictor",
  trained = FALSE,
  deg_free = 10,
  degree = 3,
  complete_set = FALSE,
  options = NULL,
  keep_original_cols = FALSE,
  results = NULL,
  skip = FALSE,
  id = rand_id("spline_nonnegative")
)

参数

recipe

一个菜谱对象。该步骤将添加到此配方的操作序列中。

...

一个或多个选择器函数用于为此步骤选择变量。有关更多详细信息，请参阅selections()。

role

对于此步骤创建的模型项，应为其分配什么分析角色？默认情况下，此步骤根据原始变量创建的新列将用作模型中的预测变量。

trained

指示预处理数量是否已估计的逻辑。

deg_free

b-spline 的自由度。随着 b-spline 自由度的增加，可以生成更灵活和复杂的曲线。

degree

指定分段多项式次数的非负整数。三次样条的默认值为 3。分段常数基函数允许零度。

complete_set

如果TRUE，将返回完整的基础矩阵。否则，第一个基础将从输出中排除。这映射到intercept相应函数的参数样条2包并具有相同的默认值。

options

splines2::mSpline() 的选项列表，不应包含 x 、 df 、 degree 、 periodic 或 intercept 。

keep_original_cols

将原始变量保留在输出中的逻辑。默认为 FALSE 。

results

训练步骤后创建的对象列表。

skip

一个合乎逻辑的。当bake() 烘焙食谱时是否应该跳过此步骤？虽然所有操作都是在 prep() 运行时烘焙的，但某些操作可能无法对新数据进行(例如处理结果变量)。使用skip = TRUE时应小心，因为它可能会影响后续操作的计算。

id

该步骤特有的字符串，用于标识它。

值

具有类 "step_spline_nonnegative" 和 "step" 的对象。

细节

样条变换采用数字列并创建多个特征，当在模型中使用这些特征时，可以估计列和某些结果之间的非线性趋势。自由度决定了向数据添加多少新特征。

该函数生成M-splines(Curry 和 Schoenberg 1988)，它是非负的并且具有有趣的统计特性(例如积分为 1)。 zero-degree M-spline 生成框/阶跃函数，而一次基函数是三角形的。

在传递给 options 的列表中设置 periodic = TRUE ，使用样条的周期版本。

如果所选列的样条扩展失败，该步骤将删除该列的结果(但将保留原始数据)。使用tidy() 方法确定使用了哪些列。

整理

当您tidy()此步骤时，将返回带有列terms(将受影响的列)的tibble。

调整参数

此步骤有 2 个调整参数：

• deg_free：样条曲线自由度(类型：整数，默认值：10)

• degree ：多项式次数(类型：整数，默认值：3)

箱重

底层操作不允许使用案例权重。

参考

库里 (H.B.)、勋伯格 (I.J.) (1988)。关于 Polya 频率函数 IV：基本样条函数及其极限。见：de Boor, C.(编)I. J. Schoenberg 论文选集。当代数学家。伯克豪瑟，波士顿，马萨诸塞州

Ramsay, J. O.“单调回归样条的实际应用”。统计科学，卷。 3、没有。 4，1988 年，第 425-41 页

也可以看看

例子

library(tidyr)
library(dplyr)

library(ggplot2)
data(ames, package = "modeldata")

spline_rec <- recipe(Sale_Price ~ Longitude, data = ames) %>%
  step_spline_nonnegative(Longitude, deg_free = 6, keep_original_cols = TRUE) %>%
  prep()

tidy(spline_rec, number = 1)
#> # A tibble: 1 × 2
#>   terms     id                      
#>   <chr>     <chr>                   
#> 1 Longitude spline_nonnegative_KVV9Z

# Show where each feature is active
spline_rec %>%
  bake(new_data =  NULL,-Sale_Price) %>%
  pivot_longer(c(starts_with("Longitude_")), names_to = "feature", values_to = "value") %>%
  mutate(feature = gsub("Longitude_", "feature ", feature)) %>%
  filter(value > 0) %>%
  ggplot(aes(x = Longitude, y = value)) +
  geom_line() +
  facet_wrap(~ feature)