R: Obtenir AIC/BIC/Probabilité de GLMNet

Question

Je me demande si je peux obtenir AIC et BIC de GLMNet. J'ai trouvé glmnet.cr qui semble être capable de le faire mais ma réponse est le temps, pas ordinale. Je pourrais le calculer moi-même à partir de la probabilité, mais glmnet ne le renvoie pas non plus. Tangentiel: puis-je retourner le l1norm? Je me sens comme il devrait juste être

fit$norm 

mais il ne semble pas être. (Je sais qu'il est dit de ne pas tirer les chiffres, mais je n'utilise pas le R)

Merci d'avance pour l'aide.

Answer 1

Je luttais beaucoup avec un moyen de calculer AIC et BIC pour les modèles glmnet. Cependant, après beaucoup de recherche, j'ai trouvé sur la troisième page de résultats google la réponse. Il peut être trouvé here. Je l'affiche ici pour les futurs lecteurs car je crois que je ne peux pas être le seul.

En fin de compte, je mis en œuvre l'AIC et BIC de la manière suivante:

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