# The Living Thing / Notebooks : Information criteria, degrees of freedom, complexity penalties

Understanding the “Degrees of freedom” of a model. Estimating that trace penalty matrix.

TODO: Explain AIC, $$C_p$$, SURE, and BIC-type degrees of freedom, and whatever variants there are out there.

As seen in robust estimation Akaike-type information criteria, these quantities crop up in model selection. (i.e. choosing the complexity of model appropriate to your data.)

Efron (Efro04) is an excellent introduction, compressing 30 years of theory into 2 pages. Massart (Mass00) seems more modern in flavour to me:

The reader who is not familiar with model selection via [complexity] penalization can legitimately ask the question: where does the idea of penalization come from? It is possible to answer this question at two different levels:

• at some intuitive level by presenting the heuristics of one of the first criterion of this kind which has been introduced by Akaike (1973);
• at some technical level by explaining why such a strategy of model selection has some chances to succeed.

Yuan and Lin (YuLi06) are an example of the kind of argumentation I need to use to use linear model approximation for general application of DOF in sparse model selection.

Zou et al (ZoHT07):

Degrees of freedom is a familiar phrase for many statisticians. In linear regression the degrees of freedom is the number of estimated predictors. Degrees of freedom is often used to quantify the model complexity of a statistical modeling procedure (Hastie and Tibshirani HaTi90). However, generally speaking, there is no exact correspondence between the degrees of freedom and the number of parameters in the model (Ye, Ye98). […] Stein’s unbiased risk estimation (SURE) theory (Stei81) gives a rigorous definition of the degrees of freedom for any fitting procedure. […] Efron (Efro04) showed that $$C_p$$ is an unbiased estimator of the true prediction error, and in some settings it offers substantially better accuracy than cross-validation and related nonparametric methods. Thus degrees of freedom plays an important role in model assessment and selection. Donoho and Johnstone (DoJo95) used the SURE theory to derive the degrees of freedom of soft thresholding and showed that it leads to an adaptive wavelet shrinkage procedure called SureShrink. Ye (Ye98) and Shen and Ye (ShYe02) showed that the degrees of freedom can capture the inherent uncertainty in modeling and frequentist model selection. Shen and Ye (ShYe02) and Shen, Huang and Ye (ShHY04) further proved that the degrees of freedom provides an adaptive model selection criterion that performs better than the fixed-penalty model selection criteria.

## Information Criteria

Akaike and friends. With M-estimation, (e.g. maximum likelihood estimation and robust estimation) these are marvelous and general shortcuts to do model selection. (i.e. choosing the complexity of model appropriate to your data) without resorting to computionally expensive cross validation.

For all of these, a thing called the number of effective degrees of freedom is important. There are several different definitions for that, and they only sometimes coincide, so I leave that for a different notebook.

Information criteria can effectively then e.g. do the same thing cross-validation at a small fraction of the computational cost. In fact, they are asymptotically the same - see below.

Estimated cross entropy (KL-divergence between model and… what?) seems to be the flavour of the minute in machine learning. (Because cross validation is unbearably slow?) Here is Chrosopher Olah’s excellent visual explanation of it These do relate, no?

To learn:

• How this interacts with robust estimators
• How to use AIC with nonparametric or high dimensional methods (GIC)
• How it relates to minimum description length (e.g. BHLL08)

Influential current English-language texts in this area are BuAn02, ClHj08 and KoKi08. The first of these is highly cited and brought the AIC method into the mainstream in the West from where it had been on the specalised fringes. The latter two focus on extensions such as TIC and GIC.

TBD: general description.

TBD: clarify relationship to Minimum Description Length, Rissanen-style.

In the literature, selection criteria are usually classified into two categories: consistent (e.g., the Bayesian information criterion BIC, Schwarz, 1978) and efficient (e.g., the Akaike information criterion AIC, Akaike, 1974; the generalized cross-validation GCV, Craven and Wahba, 1979). A consistent criterion identifies the true model with a probability that approaches 1 in large samples when a set of candidate models contains the true model. An efficient criterion selects the model so that its average squared error is asymptotically equivalent to the minimum offered by the candidate models when the true model is approximated by a family of candidate models. Detailed discussions on efficiency and consistency can be found in Shibata (1981, 1984), Li (1987), Shao (1997) and McQuarrie and Tsai (1998).

## Efficient: AIC etc

The classic.

### Takeuchi Information Criterion (TIC)

Apparently this one was influential in Japan, but untranslated into English, so only belately common in the west. Good explanations are in ClHj08 and KoKi08. Relaxes the assumption that the model is Fisher efficient (i.e. that the true generating process is included in your model, and with enough data you’d discover that.)

### Konishi and Kitegawa’s Generalised Information Criterion (GIC)

Taking information criteria to general (e.g. robust, penalised) M-estimation instead of purely ML estimation; also relaxing the assumption that we even have the “true” model in our class. (KoKi96); C&C BuNo95, probably others. In paricular, you are no longer trying to fit the midel by minimising least-squares errors, for example. ClHj08 mention the “Robustified Information Criterion” in passing, which may relate?

TBD: Explain my laborious reasoning that generalised Akaike information criteria for penalised regression don’t seem work when the penalty term is not differentiable, (cross validation works fine though, abd possibly also BIC) and the issues that therefore arise in model selection for such models in the sparse case.

### Focussed information criterion (FIC)

Claeskens and Hjort define this (ClHj08, chapter 6):

The model selection methods presented earlier (such as AIC and the BIC) have one thing in common: they select one single ‘best model’, which should then be used to explain all aspects of the mechanisms underlying the data and predict all future data points. The tolerance discussion in chapter 5 showed that sometimes one model is best for estimating one type of estimand, whereas another model is best for another estimand. The point of view expressed via the [FIC] is that a ‘best model’ should depend on the parameter under focus, such as the mean, or the variance, or the particular covariate values etc Thus the FIC allows and encourages different models to be selected for different parameters of interest.

This sounds very logical; of course, then one must do more work to make it go.

### Network information criterion

MuYA94: “an estimator of the expected loss of a loss function $$\ell(\theta)+\lambda H(\theta)$$ where $$\H(\theta)$$ is a regularisation term”.

### Regularization information criterion

Shib89 - is this distinct from GIC?

### Bootstrap information criterion

A compromise between the computational cheapness of information criteria and the analytic simplicity of cross-validation.

KoKi08 ch 8. See ClHj08 6.3 for a bootstrap-FIC.

### Consistency of model order selected - AIC

Akaike Information criteria are not asymptotically consistent (see KoKi08.) in the sense that if there is a true model, you do not asymptotically select in the large sample limit with P=1. However, the distribution of model orders does not get worse as n increases. BuAn02 6.3 and KoKi08 3.5.2 discuss this; In a sense it would be surprising if it did do especially well in selecting model order; since our criteria is rather designed to minimise prediction error, not model selection error. Model order is more or less a nuisance parameter in this framework.

TBC.

### Cross-validation equivalence

KoKi08, 10.1.4 discuss the asymptotic equivalence of AIC/TIC/GIC and cross validation under various circumstances, attributing the equivalence results to Ston77 and Shib89. ClHj08 proves a similar result.

### Automatic GIC

TBD; I know that KoKi96 give formulae for loss functions for ANY M-estimation and penalisation procedure, but in general the degrees of freedom matrix trace calculation is nasty, and only in-principle estimable from the data, requiring a matrix product of the hessian at every data point. This is not necessarily computationally tractable - I know of formulae only for GLMs and robust regression with $$\ell_2$$ penalties. Can we get such penalties for more general ML fits?

### GIC& the LASSO

I thought this didn’t work because we needed the second derivative of the penalty; but see ZhLT10.

### Information criteria at scale

Big-data information criteria. AIC is already computationally cheaper than cross validation. What about when my data is so large that I would like to select my mode before looking at all of it with such-and-such a guarantee of goodness? Can I do AIC at scale? if I am fitting a model using SGD, can I estimate my model order using partial data? How? I’m interested in doing this in a way that preserves the property of being computationally cheaper than cross-validating.

Here’s an example… Bondell et al (BoKG10):

In order to avoid complete enumeration of all possible $$2^{p+q}$$ models, Wolfinger (1993) and Diggle, Liang and Zeger (1994) recommended the Restricted Information Criterion (denoted by REML.IC), in that, by using the most complex mean structure, selection is first performed on the variance-covariance structure by computing the AIC and/or BIC. Given the best covariance structure, selection is then performed on the fixed effects. Alternatively, Pu and Niu (2006) proposed the EGIC (Extended GIC), where using the BIC, selection is first performed on the fixed effects by including all of the random effects into the model. Once the fixed effect structure is chosen, selection is then performed on the random effects.

In general I’d like to avoid enumerating the model as much as possible and simply select relevent predictors with high probability, compresssive-sensing style.

## Consistent: Bayesian Information Criteria

a.k.a. Schwarz Information Criterion. Also co-invented by the unstoppable Akaike. (Schw78, Akai78)

This is a different family to the original AIC. This has a justification in terms of MDL and of Bayes risk? Different regularity conditions, something something…

How would this work with regularisation? Apparently Mach93 gives some extensions to introduce robust fitting extensions much as the GIC introduces such to the AIC. ClHj08 give an easy summary and more general settings.

## Consistent and/or efficient: Nishii’s Generalised Information Criterion

Nish84, commended by ZhLT10 as a unifying formalism for these efficient/consistent others, which includes efficient and consistent-type information penalties as special cases I don’t know much about this, but like everything

## Refs

Akai73a
Akaike, H. (1973a) Information Theory and an Extension of the Maximum Likelihood Principle. In P. F. Caski (Ed.), Proceeding of the Second International Symposium on Information Theory (pp. 199–213). Budapest: Akademiai Kiado
Akai73b
Akaike, H. (1973b) Maximum likelihood identification of Gaussian autoregressive moving average models. Biometrika, 60(2), 255–265. DOI.
Akai78
Akaike, H. (1978) A new look at the Bayes procedure. Biometrika, 65(1), 53–59. DOI.
Akai81
Akaike, H. (1981) Likelihood of a model and information criteria. Journal of Econometrics, 16(1), 3–14. DOI.
AnKI08
Ando, T., Konishi, S., & Imoto, S. (2008) Nonlinear regression modeling via regularized radial basis function networks. Journal of Statistical Planning and Inference, 138(11), 3616–3633. DOI.
Barr86
Barron, A. R.(1986) Entropy and the Central Limit Theorem. The Annals of Probability, 14(1), 336-342. DOI.
BHLL08
Barron, A. R., Huang, C., Li, J. Q., & Luo, X. (2008) MDL, penalized likelihood, and statistical risk. In Information Theory Workshop, 2008. ITW’08. IEEE (pp. 247–257). IEEE DOI.
BaRY98
Barron, A., Rissanen, J., & Yu, B. (1998) The minimum description length principle in coding and modeling. IEEE Transactions on Information Theory, 44(6), 2743–2760. DOI.
BaHy01
Bashtannyk, D. M., & Hyndman, R. J.(2001) Bandwidth selection for kernel conditional density estimation. Computational Statistics & Data Analysis, 36(3), 279–298. DOI.
BiMa06
Birgé, L., & Massart, P. (2006) Minimal Penalties for Gaussian Model Selection. Probability Theory and Related Fields, 138(1–2), 33–73. DOI.
BoKG10
Bondell, H. D., Krishna, A., & Ghosh, S. K.(2010) Joint Variable Selection for Fixed and Random Effects in Linear Mixed-Effects Models. Biometrics, 66(4), 1069–1077. DOI.
BuBA97
Buckland, S. T., Burnham, K. P., & Augustin, N. H.(1997) Model Selection: An Integral Part of Inference. Biometrics, 53(2), 603–618. DOI.
Bune04
Bunea, F. (2004) Consistent covariate selection and post model selection inference in semiparametric regression. The Annals of Statistics, 32(3), 898–927. DOI.
BuNo95
Burman, P., & Nolan, D. (1995) A general Akaike-type criterion for model selection in robust regression. Biometrika, 82(4), 877–886. DOI.
BuAn02
Burnham, K. P., & Anderson, D. R.(2002) Model selection and multimodel inference: a practical information-theoretic approach. (2nd ed.). New York: Springer
BuAn04
Burnham, K. P., & Anderson, D. R.(2004) Multimodel Inference Understanding AIC and BIC in Model Selection. Sociological Methods & Research, 33(2), 261–304. DOI.
Cava97
Cavanaugh, J. E.(1997) Unifying the derivations for the Akaike and corrected Akaike information criteria. Statistics & Probability Letters, 33(2), 201–208. DOI.
CaSh98
Cavanaugh, J. E., & Shumway, R. H.(1998) An Akaike information criterion for model selection in the presence of incomplete data. Journal of Statistical Planning and Inference, 67(1), 45–65. DOI.
ChCh08
Chen, J., & Chen, Z. (2008) Extended Bayesian information criteria for model selection with large model spaces. Biometrika, 95(3), 759–771. DOI.
ClHj08
Claeskens, G., & Hjort, N. L.(2008) Model selection and model averaging. . Cambridge ; New York: Cambridge University Press
ClKO09
Claeskens, G., Krivobokova, T., & Opsomer, J. D.(2009) Asymptotic properties of penalized spline estimators. Biometrika, 96(3), 529–544. DOI.
DoJo95
Donoho, D. L., & Johnstone, I. M.(1995) Adapting to Unknown Smoothness via Wavelet Shrinkage. Journal of the American Statistical Association, 90(432), 1200–1224. DOI.
Efro86
Efron, B. (1986) How biased is the apparent error rate of a prediction rule?. Journal of the American Statistical Association, 81(394), 461–470. DOI.
Efro04
Efron, B. (2004) The Estimation of Prediction Error. Journal of the American Statistical Association, 99(467), 619–632. DOI.
FaLi01
Fan, J., & Li, R. (2001) Variable Selection via Nonconcave Penalized Likelihood and its Oracle Properties. Journal of the American Statistical Association, 96(456), 1348–1360. DOI.
HaTi90
Hastie, T. J., & Tibshirani, R. J.(1990) Generalized additive models. (Vol. 43). CRC Press
HuCB08
Huang, C., Cheang, G. L. H., & Barron, A. R.(2008) Risk of penalized least squares, greedy selection and l1 penalization for flexible function libraries.
HMXZ09
Huang, J., Ma, S., Xie, H., & Zhang, C.-H. (2009) A group bridge approach for variable selection. Biometrika, 96(2), 339–355. DOI.
HuST98
Hurvich, C. M., Simonoff, J. S., & Tsai, C.-L. (1998) Smoothing Parameter Selection in Nonparametric Regression Using an Improved Akaike Information Criterion. Journal of the Royal Statistical Society. Series B (Statistical Methodology), 60(2), 271–293.
HuTs89
Hurvich, C. M., & Tsai, C.-L. (1989) Regression and time series model selection in small samples. Biometrika, 76(2), 297–307. DOI.
ImKo99
Imoto, S., & Konishi, S. (1999) Estimation of B-spline Nonparametric Regression Models using Information.
JaFH15
Janson, L., Fithian, W., & Hastie, T. J.(2015) Effective degrees of freedom: a flawed metaphor. Biometrika, 102(2), 479–485. DOI.
KaRo14
Kaufman, S., & Rosset, S. (2014) When does more regularization imply fewer degrees of freedom? Sufficient conditions and counterexamples. Biometrika, 101(4), 771–784. DOI.
KoKi96
Konishi, S., & Kitagawa, G. (1996) Generalised information criteria in model selection. Biometrika, 83(4), 875–890. DOI.
KoKi03
Konishi, S., & Kitagawa, G. (2003) Asymptotic theory for information criteria in model selection—functional approach. Journal of Statistical Planning and Inference, 114(1–2), 45–61. DOI.
KoKi08
Konishi, S., & Kitagawa, G. (2008) Information criteria and statistical modeling. . New York: Springer
LeBa06
Leung, G., & Barron, A. R.(2006) Information Theory and Mixing Least-Squares Regressions. IEEE Transactions on Information Theory, 52(8), 3396–3410. DOI.
LiBa00
Li, J. Q., & Barron, A. R.(2000) Mixture Density Estimation. In S. A. Solla, T. K. Leen, & K. Müller (Eds.), Advances in Neural Information Processing Systems 12 (pp. 279–285). MIT Press
Li87
Li, K.-C. (1987) Asymptotic Optimality for $C_p, C_L$, Cross-Validation and Generalized Cross-Validation: Discrete Index Set. The Annals of Statistics, 15(3), 958–975. DOI.
LiLi08
Li, R., & Liang, H. (2008) Variable selection in semiparametric regression modeling. The Annals of Statistics, 36(1), 261–286. DOI.
Mach93
Machado, J. A. F.(1993) Robust Model Selection and M-Estimation. Econometric Theory, 9(3), 478–493. DOI.
Mass00
Massart, P. (2000) Some applications of concentration inequalities to statistics. In Annales de la Faculté des sciences de Toulouse: Mathématiques (Vol. 9, pp. 245–303).
Mass07
Massart, P. (2007) Concentration inequalities and model selection: Ecole d’Eté de Probabilités de Saint-Flour XXXIII - 2003. . Berlin ; New York: Springer-Verlag
MuYA94
Murata, N., Yoshizawa, S., & Amari, S. (1994) Network information criterion-determining the number of hidden units for an artificial neural network model. IEEE Transactions on Neural Networks, 5(6), 865–872. DOI.
Nish84
Nishii, R. (1984) Asymptotic Properties of Criteria for Selection of Variables in Multiple Regression. The Annals of Statistics, 12(2), 758–765. DOI.
QiHa96
Qian, G., & Hans, R. K.(1996) Some notes on Rissanen’s stochastic complexity.
QiKü98
Qian, G., & Künsch, H. R.(1998) On model selection via stochastic complexity in robust linear regression. Journal of Statistical Planning and Inference, 75(1), 91–116. DOI.
RaWu01
Rao, C. R., & Wu, Y. (2001) On model selection. In Institute of Mathematical Statistics Lecture Notes - Monograph Series (Vol. 38, pp. 1–57). Beachwood, OH: Institute of Mathematical Statistics
RaWu89
Rao, R., & Wu, Y. (1989) A strongly consistent procedure for model selection in a regression problem. Biometrika, 76(2), 369–374. DOI.
Riss78
Rissanen, J. (1978) Modeling by shortest data description. Automatica, 14(5), 465–471. DOI.
SKWG14
Saefken, B., Kneib, T., Waveren, C.-S. van, & Greven, S. (2014) A unifying approach to the estimation of the conditional Akaike information in generalized linear mixed models. Electronic Journal of Statistics, 8(1), 201–225. DOI.
Schw78
Schwarz, G. (1978) Estimating the Dimension of a Model. The Annals of Statistics, 6(2), 461–464. DOI.
ShHu06
Shen, X., & Huang, H.-C. (2006) Optimal Model Assessment, Selection, and Combination. Journal of the American Statistical Association, 101(474), 554–568. DOI.
ShHY04
Shen, X., Huang, H.-C., & Ye, J. (2004) Adaptive Model Selection and Assessment for Exponential Family Distributions. Technometrics, 46(3), 306–317. DOI.
ShYe02
Shen, X., & Ye, J. (2002) Adaptive Model Selection. Journal of the American Statistical Association, 97(457), 210–221. DOI.
Shib89
Shibata, R. (1989) Statistical Aspects of Model Selection. In P. J. C. Willems (Ed.), From Data to Model (pp. 215–240). Springer Berlin Heidelberg DOI.
Stei81
Stein, C. M.(1981) Estimation of the Mean of a Multivariate Normal Distribution. The Annals of Statistics, 9(6), 1135–1151. DOI.
Ston77
Stone, M. (1977) An Asymptotic Equivalence of Choice of Model by Cross-Validation and Akaike’s Criterion. Journal of the Royal Statistical Society. Series B (Methodological), 39(1), 44–47.
Sugi78
Sugiura, N. (1978) Further analysts of the data by Akaike’ s Information Criterion and the finite corrections. Communications in Statistics - Theory and Methods, 7(1), 13–26. DOI.
Taddy, M. (2013) One-step estimator paths for concave regularization. arXiv:1308.5623 [Stat].
Tarm11
Tarmaratram, K. (2011) Robust Estimation and Model Selection in Semiparametric Regression Models. . Katholieke Universiteit Leuven
ThCl13
Tharmaratnam, K., & Claeskens, G. (2013) A comparison of robust versions of the AIC based on M-, S- and MM-estimators. Statistics, 47(1), 216–235. DOI.
Tibs96
Tibshirani, R. (1996) Regression Shrinkage and Selection via the Lasso. Journal of the Royal Statistical Society. Series B (Methodological), 58(1), 267–288.
Ye98
Ye, J. (1998) On Measuring and Correcting the Effects of Data Mining and Model Selection. Journal of the American Statistical Association, 93(441), 120–131. DOI.
YuLi06
Yuan, M., & Lin, Y. (2006) Model selection and estimation in regression with grouped variables. Journal of the Royal Statistical Society: Series B (Statistical Methodology), 68(1), 49–67. DOI.
ZhLT10
Zhang, Y., Li, R., & Tsai, C.-L. (2010) Regularization Parameter Selections via Generalized Information Criterion. Journal of the American Statistical Association, 105(489), 312–323. DOI.
ZoHT07
Zou, H., Hastie, T., & Tibshirani, R. (2007) On the “degrees of freedom” of the lasso. The Annals of Statistics, 35(5), 2173–2192. DOI.