是否可以计算通过glmnet拟合的岭/套索逻辑回归模型的边际效应？

Question

似乎不支持应用

margins

包从 glmnet 模型中提取

平均边际效应

。我无法从

glm

复制相同的行为，因为我正在拟合岭回归，并且（我认为）后者不支持此计算。是否有任何选项可以计算逻辑岭回归模型的边际效应？

Answer 1

marginaleffects

的开发版本（即将在CRAN上发布0.16.0版本）支持

tidymodels

。首先，安装软件包：

remotes::install_github("vincentarelbundock/marginaleffects")

然后，完全重新启动

R

。

最后，请注意，仅当通过

glmnet

的公式接口（不是矩阵形式）指定拟合时才支持

tidymodels

：

library(tidymodels)
library(marginaleffects)

penguins <- modeldata::penguins |> na.omit()

mod <- linear_reg(mode = "regression", penalty = 1) |>
    set_engine("glmnet") |>
    fit(bill_length_mm ~ ., data = penguins)

avg_slopes(mod, newdata = penguins)
# 
#               Term           Contrast Estimate
#  bill_depth_mm     dY/dX               0.00000
#  bill_length_mm    dY/dX               0.00000
#  body_mass_g       dY/dX               0.00142
#  flipper_length_mm dY/dX               0.14098
#  island            Dream - Biscoe      0.00000
#  island            Torgersen - Biscoe  0.00000
#  sex               male - female       0.00000
#  species           Chinstrap - Adelie  5.35902
#  species           Gentoo - Adelie     0.00000
# 
# Columns: term, contrast, estimate 
# Type:  numeric