losglmnet el paquete usa un rango deLASSO parámetros de ajustelambda escalado desde el máximolambda_max bajo el cual no se seleccionan predictores. Quiero saber comoglmnet calcula estolambda_max valor. Por ejemplo, en un conjunto de datos trivial:

set.seed(1)
library("glmnet")
x <- matrix(rnorm(100*20),100,20)
y <- rnorm(100)
fitGLM <- glmnet(x,y)
max(fitGLM$lambda)
# 0.1975946

La viñeta del paquete (http://www.jstatsoft.org/v33/i01/paper) describe en la sección 2.5 que calcula este valor de la siguiente manera:

sx <- as.matrix(scale(x))
sy <- as.vector(scale(y))
max(abs(colSums(sx*sy)))/100
# 0.1865232

Que claramente está cerca pero no tiene el mismo valor. Entonces, ¿qué causa esta diferencia? Y en una pregunta relacionada, ¿cómo podría calcularlambda_max para una regresión logística?