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

Question

我修改了这篇文章中的函数（在 R 中的 glmnet 图中添加曲线上的标签）以将图例添加到图中，如下所示：

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
}

Answer 2

您可以在plot调用中指定颜色以匹配您为图例选择的颜色：

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", col=1:dim(coef(fit))[1])
lbs_fun(fit)