Gracias por la respuesta de @ TommyL, pero su respuesta no es directa sobre la construcción de e . De alguna manera "resuelvo" esto yo mismo. Primero, cuando aumenta, no aumentará cuando cada disminuya monotónicamente. Esto sucede cuando es ortonormal, en el que tenemosXyλ∥β∗∥2β∗iX
β∗i=sign(βLSi)(βLSi−λ)+
Geométricamente, en esta situación, mueve perpendicularmente al contorno de la norma , por lo que no puede aumentar.β∗ℓ1∥β∗∥2
En realidad, Hastie et al. mencionado en el artículo La regresión progresiva por etapas y el lazo monótono , una condición necesaria y suficiente de la monotonicidad de las rutas de perfil:
En la Sección 6 del artículo, construyeron un conjunto de datos artificiales basado en funciones de base lineal por partes que viola la condición anterior, mostrando la no monotonicidad. Pero si tenemos suerte, también podemos crear un conjunto de datos aleatorios que demuestren un comportamiento similar pero de una manera más simple. Aquí está mi código R:
library(glmnet)
set.seed(0)
N <- 10
p <- 15
x1 <- rnorm(N)
X <- mat.or.vec(N, p)
X[, 1] <- x1
for (i in 2:p) {X[, i] <- x1 + rnorm(N, sd=0.2)}
beta <- rnorm(p, sd=10)
y <- X %*% beta + rnorm(N, sd=0.01)
model <- glmnet(X, y, family="gaussian", alpha=1, intercept=FALSE)
Deliberadamente dejé que las columnas de estuvieran altamente correlacionadas (lejos del caso ortonormal), y el verdadero tiene grandes entradas positivas y negativas. Aquí está el perfil de (no es sorprendente que solo se activen 5 variables):Xββ∗
y la relación entre y :λ∥β∗∥2
Entonces podemos ver que durante algún intervalo de , aumenta a medida que aumenta.λ∥β∗∥2λ