R embed step_pca_sparse 稀疏PCA信号提取

step_pca_sparse() 创建配方步骤的规范，该步骤将数值数据转换为一个或多个可以具有一些零系数的主成分。

用法

step_pca_sparse(
  recipe,
  ...,
  role = "predictor",
  trained = FALSE,
  num_comp = 5,
  predictor_prop = 1,
  options = list(),
  res = NULL,
  prefix = "PC",
  keep_original_cols = FALSE,
  skip = FALSE,
  id = rand_id("pca_sparse")
)

参数

recipe

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

...

一个或多个选择器函数用于选择将使用哪些变量来计算组件。有关更多详细信息，请参阅selections()。对于tidy 方法，当前未使用这些。

role

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

trained

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

num_comp

保留作为新预测变量的组件数量。如果num_comp大于列数或可能组件的数量，则将使用较小的值。如果设置了 num_comp = 0 ，则不会进行任何转换，并且所选变量将保持不变，无论 keep_original_cols 的值如何。

predictor_prop

每个 PCA 分量可以具有非零系数的原始预测变量的最大数量(通过正则化)。

options

irlba::ssvd() 默认方法的选项列表。

res

此预处理步骤经过 prep() 训练后的旋转矩阵。

prefix

将作为结果新变量的前缀的字符串。请参阅下面的注释。

keep_original_cols

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

skip

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

id

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

值

recipe 的更新版本，其中新步骤添加到现有步骤(如果有)的序列中。对于 tidy 方法，包含列 terms (选择的选择器或变量)、 value (加载)和 component 的小标题。

细节

此步骤需要 irlba 包。如果未安装，则在定义步骤时将提示用户安装。 irlba::ssvd() 函数用于鼓励稀疏性；该文档包含有关此方法的详细信息。

参数 num_comp 控制将保留的组件数量(用于派生组件的原始变量将从数据中删除)。新组件的名称以 prefix 和一系列数字开头。变量名称用零填充。例如，如果 num_comp < 10 ，它们的名称将为 PC1 - PC9 。如果是 num_comp = 101 ，则名称将为 PC1 - PC101 。

整理

当您tidy()此步骤时，将返回包含terms(选定的选择器或变量)、value 和component 列的小标题。

调整参数

此步骤有 2 个调整参数：

• num_comp : # 组件(类型：整数，默认值：5)

• predictor_prop：预测变量的比例(类型：double，默认值：1)

箱重

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

也可以看看

例子

library(recipes)
library(ggplot2)

data(ad_data, package = "modeldata")

ad_rec <-
  recipe(Class ~ ., data = ad_data) %>%
  step_zv(all_predictors()) %>%
  step_YeoJohnson(all_numeric_predictors()) %>%
  step_normalize(all_numeric_predictors()) %>%
  step_pca_sparse(
    all_numeric_predictors(),
    predictor_prop = 0.75,
    num_comp = 3,
    id = "sparse pca"
  ) %>%
  prep()

tidy(ad_rec, id = "sparse pca") %>%
  mutate(value = ifelse(value == 0, NA, value)) %>%
  ggplot(aes(x = component, y = terms, fill = value)) +
  geom_tile() +
  scale_fill_gradient2() +
  theme(axis.text.y = element_blank())