r - R 中 glmnet 图的图例标签错误

Question

我修改了这篇文章 (Adding labels on curves in glmnet plot in R) 中的函数，将图例添加到图中，如下所示:

library(glmnet)
fit = glmnet(as.matrix(mtcars[-1]), mtcars[,1])

lbs_fun <- function(fit, ...) {
        L <- length(fit$lambda)
        x <- log(fit$lambda[L])
        y <- fit$beta[, L]
        labs <- names(y)
        text(x, y, labels=labs, ...)
        legend('topright', legend=labs, col=1:length(labs), lty=1) # <<< ADDED BY ME
}
plot(fit, xvar="lambda")
lbs_fun(fit)

但是，我发现绘图和图例中的文本标签不匹配。变量“am”的颜色显然不正确。错误在哪里？感谢您的帮助。

Answer 1

plot(fit, xvar="lambda") 利用函数 matplot。默认情况下，matplot 使用 6 种颜色并回收它们。所以你必须相应地创建图例:

lbs_fun <- function(fit, ...) {
        L <- length(fit$lambda)
        x <- log(fit$lambda[L])
        y <- fit$beta[, L]
        labs <- names(y)
        text(x, y, labels=labs, ...)
        legend('topright', legend=labs, col=1:6, lty=1) # only 6 colors
}