bic_score: Calculates the BIC score.

In sqyu/ZiDAG: Directed Graphical Models and Causal Discovery for Zero-Inflated Data

Description

Calculates the BIC score.

Usage

bic_score(avgnll, n, numpar)

Arguments

avgnll	Average negative log likelihood over the sample.
n	The sample size.
numpar	The degree of freedom.

Details

Returns 2*avgnll+log(n)/n*numpar.

Value

A number, the BIC score.

Examples

Y <- matrix(rnorm(500) * rbinom(500, 1, 0.7), nrow=50)
res <- zi_fit(Y != 0, Y, 1, 2:5, value_only=FALSE)
bic_score(res$nll, res$n, res$effective_df)
bic_from_fit(res)

sqyu/ZiDAG documentation built on Jan. 19, 2021, 4:11 p.m.