mclust (version 3.1-0)

bic: BIC for Parameterized Gaussian Mixture Models

Description

Computes the BIC (Bayesian Information Criterion) for parameterized mixture models given the loglikelihood, the dimension of the data, and number of mixture components in the model.

Usage

bic(modelName, loglik, n, d, G, noise=FALSE, equalPro=FALSE, ...)

Arguments

modelName

A character string indicating the model. The help file for mclustModelNames describes the available models.

loglik

The loglikelihood for a data set with respect to the Gaussian mixture model specified in the modelName argument.

n

The number of observations in the data used to compute loglik.

d

The dimension of the data used to compute loglik.

G

The number of components in the Gaussian mixture model used to compute loglik.

noise

A logical variable indicating whether or not the model includes an optional Poisson noise component. The default is to assume no noise component.

equalPro

A logical variable indicating whether or not the components in the model are assumed to be present in equal proportion. The default is to assume unequal mixing proportions.

...

Catches unused arguments in an indirect or list call via do.call.

Value

• The BIC or Bayesian Information Criterion for the given input arguments.

References

C. Fraley and A. E. Raftery (2002). Model-based clustering, discriminant analysis, and density estimation. Journal of the American Statistical Association 97:611:631.

C. Fraley and A. E. Raftery (2006). MCLUST Version 3 for R: Normal Mixture Modeling and Model-Based Clustering, Technical Report no. 504, Department of Statistics, University of Washington.

See Also

nVarParams, mclustBIC, do.call.

Examples

n <- nrow(iris)
d <- ncol(iris)-1
G <- 3

emEst <- me(modelName="VVI", data=iris[,-5], unmap(iris[,5]))
names(emEst)

args(bic)
bic(modelName="VVI", loglik=emEst$loglik, n=n, d=d, G=G)
do.call("bic", emEst)    ## alternative call