PLNmodels (version 1.3.1)

ICL: Integrated Classification Likelihood

Description

Generic function to compute the Integrated Classification Likelihood (ICL) of a fitted model. ICL = BIC - entropy of the variational distribution (larger is better).

ICL.PLNfit: ICL for a fitted PLNfit.

ICL.ZIPLNfit: ICL for a fitted ZIPLNfit.

Usage

ICL(object, ...)
# S3 method for PLNfit
ICL(object, ...)
# S3 method for ZIPLNfit
ICL(object, ...)

Value

A scalar: the variational ICL (larger is better).

Arguments

object

an R6 object with class ZIPLNfit

...

additional parameters passed to methods

Examples

data(trichoptera)
trichoptera <- prepare_data(trichoptera$Abundance, trichoptera$Covariate)
model <- PLN(Abundance ~ 1, data = trichoptera)
ICL(model)
data(trichoptera)
trichoptera <- prepare_data(trichoptera$Abundance, trichoptera$Covariate)
model <- ZIPLN(Abundance ~ 1, data = trichoptera)
ICL(model)