PLNmodels (version 1.3.1)

BIC.PLNfit: Bayesian Information Criterion for a fitted PLN model

Description

Computes the variational BIC as loglik - 0.5 * log(n) * nb_param (larger is better). This follows the maximization convention used throughout PLNmodels.

Usage

# S3 method for PLNfit
BIC(object, ...)

Value

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

Arguments

object

an R6 object with class PLNfit

...

additional parameters for S3 compatibility. Not used

Examples

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