R:从 GLMNet 获取 AIC/BIC/Likelihood

Question

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

Tangential:我真的可以返回 l1norm 吗？我觉得应该是这样

fit$norm

但似乎并非如此。 (我知道它说不要提取数字，但实际上我没有使用 R)

预先感谢您的帮助。

Answer 1

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

最终我通过以下方式实现了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