如何在R中迭代移除n个glmnet模型选中变量的截距项

你可以通过对Variables.Selected列表的每个元素进行过滤，批量移除其中的(Intercept)项，以下是两种贴合你现有代码风格的实现方式：

方法1：使用purrr管道风格（与你现有代码逻辑一致） #

# 移除截距项，生成仅含自变量的选中变量列表
Variables.Selected.NoIntercept <- Variables.Selected |>
  map(\(var_list) var_list[var_list != "(Intercept)"])

方法2：使用base R的lapply #

如果你更习惯base R语法，也可以这样写：

Variables.Selected.NoIntercept <- lapply(Variables.Selected, function(var_list) {
  var_list[var_list != "(Intercept)"]
})

验证结果 #

运行后查看前4个元素，就能看到截距项已经被移除，格式与Structural_Variables完全匹配：

> head(Variables.Selected.NoIntercept, n = 4)
[[1]]
 [1] "X1"  "X6"  "X7"  "X8"  "X10" "X11" "X13" "X15" "X17" "X20" "X22" "X24" "X26" "X27" "X28" "X29" "X30"

[[2]]
 [1] "X3"  "X5"  "X8"  "X9"  "X13" "X14" "X16" "X19" "X20" "X24" "X25" "X26" "X29" "X30"

[[3]]
 [1] "X1"  "X4"  "X5"  "X10" "X12" "X13" "X14" "X19" "X20" "X21" "X24" "X25" "X27" "X29"

[[4]]
 [1] "X3"  "X4"  "X5"  "X9"  "X10" "X11" "X14" "X17" "X18" "X22" "X24" "X26" "X27" "X28"

之后你就可以直接和Structural_Variables进行对比，比如用setequal()判断变量选择是否完全一致，或者计算交集、并集来评估模型性能。

内容的提问来源于stack exchange，提问作者Marlen