AICvlm {VGAM}R Documentation

Akaike's An Information Criterion

Description

Calculates the Akaike information criterion for a fitted model object for which a log-likelihood value has been obtained.

Usage

AICvlm(object, ..., k = 2)

Arguments

object Some VGAM object, for example, having class vglmff-class.
... Other possible arguments fed into logLik in order to compute the log-likelihood.
k Numeric, the penalty per parameter to be used; the default is the classical AIC.

Details

The following formula is used for VGLMs: -2*log-likelihood + k*npar, where npar represents the number of parameters in the fitted model, and k = 2 for the usual AIC. One can assign k = log(n) (n the number of observations) for the so-called BIC or SBC (Schwarz's Bayesian criterion).

This code relies on the log-likelihood being defined, and computed, for the object. When comparing fitted objects, the smaller the AIC, the better the fit. The log-likelihood and hence the AIC is only defined up to an additive constant.

Any estimated scale parameter (in GLM parlance) is used as one parameter.

For VGAMs the nonlinear effective degrees of freedom for each smoothed component is used. This is heuristic.

Value

Returns a numeric value with the corresponding AIC (or BIC, or ..., depending on k).

Warning

This code has not been double-checked. The general applicability of AIC for the VGLM/VGAM classes has not been developed fully. In particular, AIC should not be run on some VGAM family functions because of violation of certain regularity conditions, etc.

Note

AIC has not been defined for QRR-VGLMs yet.

Author(s)

T. W. Yee

See Also

VGLMs are described in vglm-class; RR-VGLMs are described in rrvglm-class; AIC.

Examples

# Fit a proportional odds model
data(pneumo)
pneumo = transform(pneumo, let=log(exposure.time))
(fit1 = vglm(cbind(normal, mild, severe) ~ let,
           cumulative(parallel=TRUE, reverse=TRUE), pneumo))
coef(fit1, matrix=TRUE)
AIC(fit1)
(fit2 = vglm(cbind(normal, mild, severe) ~ let,
           cumulative(parallel=FALSE, reverse=TRUE), pneumo))
coef(fit2, matrix=TRUE)
AIC(fit2)

[Package VGAM version 0.7-8 Index]