R：从 GLMNet 获取 AIC/BIC/可能性

Question

我想知道我是否可以从 GLMNet 获得 AIC 和 BIC。我发现 glmnet.cr 似乎能够做到，但我的反应是时间，而不是顺序。我可以根据可能性自己计算它，但 glmnet 也不会返回。

切线：我真的可以返回 l1norm 吗？我觉得它应该只是

fit$norm

但似乎不是。（我知道它说不要把数字拉出来，但我实际上没有使用 R）

在此先感谢您的帮助。

Answer 1

我在如何计算 glmnet 模型的 AIC 和 BIC 的方法上遇到了很多困难。然而，经过大量搜索后，我在谷歌结果的第三页找到了答案。在这里能找到它。我将其发布在这里供未来的读者使用，因为我相信我不是唯一的读者。

最终我按照如下方式实现了AIC和BIC：

fit <- glmnet(x, y, family = "multinomial") 

tLL <- fit$nulldev - deviance(fit)
k <- fit$df
n <- fit$nobs
AICc <- -tLL+2*k+2*k*(k+1)/(n-k-1)
AICc

BIC<-log(n)*k - tLL
BIC