Learn R Programming

qpcR (version 1.3-7.1)

akaike.weights: Calculation of Akaike weights/relative likelihoods/delta-AICs

Description

Calculates Akaike weights from a vector of AIC values.

Usage

akaike.weights(x)

Arguments

x
a vector containing the AIC values.

Value

  • A list containing the following items:
  • deltaAICthe $\Delta(AIC)$ values.
  • rel.LLthe relative likelihoods.
  • weightsthe Akaike weights.

encoding

latin1

Details

Although Akaike's Information Criterion is recognized as a major measure for selecting models, it has one major drawback: The AIC values lack intuitivity despite higher values meaning less goodness-of-fit. For this purpose, Akaike weights come to hand for calculating the weights in a regime of several models. Additional measures can be derived, such as $\Delta(AIC)$ and relative likelihoods that demonstrate the probability of one model being in favor over the other. This is done by using the following formulas: delta AICs: $$\Delta_i(AIC) = AIC_i - min(AIC)$$ relative likelihood: $$L \propto exp\left(-\frac{1}{2}\Delta_i(AIC)\right)$$ Akaike weights: $$w_i(AIC) = \frac{exp\left(-\frac{1}{2}\Delta_i(AIC)\right)}{\sum_{k=1}^K exp\left(-\frac{1}{2}\Delta_k(AIC)\right)}$$

References

Classical literature: Akaike Information Criterion Statistics. Sakamoto Y, Ishiguro M and Kitagawa G. D. Reidel Publishing Company (1986). Model selection and inference: a practical information-theoretic approach. Burnham KP & Anderson DR. Springer Verlag, New York, USA (2002). A good summary: AIC model selection using Akaike weights. Wagenmakers EJ & Farrell S. Psychonomic Bull Review (2004), 11: 192-196.

See Also

AIC, logLik.

Examples

Run this code
## apply a list of different sigmoidal models to data
## and analyze GOF statistics with Akaike weights
## on 6 different sigmoidal models 
modList <- list(l5, l4, l3, b5, b4, b3)
aics <- sapply(modList, function(x) AIC(pcrfit(reps, 1, 2, x))) 
akaike.weights(aics)$weights

Run the code above in your browser using DataLab