r - R：来自 glmnet 岭回归的系数不匹配

Question

在岭回归中，系数具有解析形式：

其中 lambda 是一个正数。

使用，我使用,glmnet执行岭回归。但是，为什么不匹配那些解析公式呢？alpha =0lambda = 1fit1$beta

library(glmnet)
set.seed(3)
x = matrix(rnorm(100 * 20), 100, 20)
y = rnorm(100)
lambda <- 1
fit1 = glmnet(x, y, alpha = 0, lambda = lambda)
> head(fit1$beta)
6 x 1 sparse Matrix of class "dgCMatrix"
             s0
V1 -0.016685052
V2  0.024373749
V3 -0.008228917
V4  0.035215872
V5 -0.008518129
V6 -0.001925812
> head(solve(t(x) %*% x + lambda * diag(p)) %*% t(x) %*% y)
            [,1]
[1,] -0.05338655
[2,]  0.04285142
[3,] -0.03454318
[4,]  0.08704770
[5,] -0.02599216
[6,] -0.02031802

标签： rregressionglmnet