gmgm (version 1.1.2)

BIC: Compute the Bayesian Information Criterion (BIC) of a Gaussian mixture model or graphical model

Description

This function computes the Bayesian Information Criterion (BIC) of a Gaussian mixture model or graphical model: $$BIC=logLik-\frac{\log (n_{obs})}{2}{n}_{par}$$ where $logLik$ is the log-likelihood, $n_{obs}$ the number of observations in the data and $n_{par}$ the number of free parameters.

Usage

# S3 method for gmm
BIC(object, data, y = NULL, regul = 0.01, ...)
# S3 method for gmbn
BIC(object, data, col_seq = NULL, ...)
# S3 method for gmdbn
BIC(object, data, col_seq = NULL, ...)

Value

If object is a gmm object, a numeric value corresponding to the BIC.

If object is a gmbn or gmdbn object, a list with elements:

global

A numeric value corresponding to the global BIC.

local

For a gmbn object, a numeric vector containing the local conditional BICs. For a gmdbn object, a list of numeric vectors containing these values for each gmbn element.

Arguments

object

An object of class gmm, gmbn or gmdbn.

data

A data frame containing the data used to compute the BIC. Its columns must explicitly be named after the variables (or nodes) of object. If object is a gmm object, a numeric matrix can be passed.

y

A character vector containing the dependent variables if a conditional BIC is computed. If NULL (the default), the joint BIC is computed.

regul

A positive numeric value corresponding to the regularization constant if a penalty term is added for Bayesian regularization. If NULL, no penalty term is added. If a conditional BIC is computed, this argument is ignored.

...

Unused arguments from the generic function.

col_seq

A character vector containing the column names of data that describe the observation sequence. If NULL (the default), all the observations belong to a single sequence. If object is a temporal gmbn or gmdbn object, the observations of a same sequence must be ordered such that the $t$th one is related to time slice $t$ (note that the sequences can have different lengths). If object is a non-temporal gmbn object, this argument is ignored.

See Also

AIC, logLik

Examples

data(gmm_body, data_body)
bic_1 <- BIC(gmm_body, data_body)
bic_2 <- BIC(gmm_body, data_body, y = "WAIST")

data(gmbn_body, data_body)
bic_3 <- BIC(gmbn_body, data_body)

data(gmdbn_air, data_air)
bic_4 <- BIC(gmdbn_air, data_air, col_seq = "DATE")