Learn R Programming

geiger (version 2.0.11)

aicm: Akaike's Information Criterion for MCMC samples (AICM)

Description

Computes Akaike's Information Criterion for MCMC samples (AICM: see Raftery et al. 2007). Can be used to perform model selection using output from fitContinuousMCMC.

Usage

aicm(x)

Value

AICM - Akaike's Information Criterion for the posterior sample

Arguments

x

a vector containing sampled likelihoods from the MCMC chain. Assumes that burn-in has been removed prior to computation of AICM score.

Author

Graham Slater

Details

AICM is one way of comparing model fit using posterior likelihood samples. It has advantages over approaches such as thermodynamic integration in that it uses the chain output directly and thus has little added time cost. Furthermore, it has been shown to perform better than the harmonic mean estimator of the marginal likelihood. However, it is also less robust than stepping-stone or thermodynamic integration approaches and should be used with care.

References

Raftery et al. 2007. Estimating the integrated likelihood via posterior simulation using the harmonic mean identity. In. Bernardo et al. (eds) Bayesian Statistics. Oxford University Press.

Examples

Run this code

## generate a random set of values from a normal distribution, 
## similar to a set of likelihood samples obtained via MCMC.

x <- rnorm(1000, -275, 2);
aicm(x);

Run the code above in your browser using DataLab