AIC.seqModel: Information criteria for a sequence of regression models In robustHD: Robust Methods for High-Dimensional Data

Description

Compute the Akaike or Bayes information criterion for for a sequence of regression models, such as submodels along a robust least angle regression sequence, or sparse least trimmed squares regression models for a grid of values for the penalty parameter.

Usage

 1 2 3 4 5 6 7 8 9 10 11 ## S3 method for class 'seqModel' AIC(object, ..., k = 2) ## S3 method for class 'sparseLTS' AIC(object, ..., fit = c("reweighted", "raw", "both"), k = 2) ## S3 method for class 'seqModel' BIC(object, ...) ## S3 method for class 'sparseLTS' BIC(object, ...)

Arguments

 object the model fit for which to compute the information criterion. ... for the BIC method, additional arguments to be passed down to the AIC method. For the AIC method, additional arguments are currently ignored. k a numeric value giving the penalty per parameter to be used. The default is to use 2 as in the classical definition of the AIC. fit a character string specifying for which fit to compute the information criterion. Possible values are "reweighted" (the default) for the information criterion of the reweighted fit, "raw" for the information criterion of the raw fit, or "both" for the information criteria of both fits.

Details

The information criteria are computed as n (log(2 pi) + 1 + log(sigma^2)) + df k, where n denotes the number of observations, sigma is the robust residual scale estimate, df is the number of nonzero coefficient estimates, and k is penalty per parameter. The usual definition of the AIC uses k = 2, whereas the BIC uses k = log(n). Consequently, the former is used as the default penalty of the AIC method, whereas the BIC method calls the AIC method with the latter penalty.

Value

A numeric vector or matrix giving the information criteria for the requested model fits.

Note

Computing information criteria for several objects supplied via the ... argument (as for the default methods of AIC and BIC) is currently not implemented.

Andreas Alfons

References

Akaike, H. (1970) Statistical predictor identification. Annals of the Institute of Statistical Mathematics, 22(2), 203–217.

Schwarz, G. (1978) Estimating the dimension of a model. The Annals of Statistics, 6(2), 461–464.