calculate_bic: calculate_bic

In pharmpy/pharmr: Interface to the 'Pharmpy' 'Pharmacometrics' Library

calculate_bic

Description

Calculate BIC

Different variations of the BIC can be calculated:

• | mixed (default) | BIC = -2LL + n_random_parameters * log(n_individuals) + | n_fixed_parameters * log(n_observations)

• | fixed | BIC = -2LL + n_estimated_parameters * log(n_observations)

• | random | BIC = -2LL + n_estimated_parameters * log(n_individuals)

• | iiv | BIC = -2LL + n_estimated_iiv_omega_parameters * log(n_individuals)

Usage

calculate_bic(model, likelihood, type = "mixed")

Arguments

model	(Model) Pharmpy model object
likelihood	(numeric) -2LL to use
type	(str) Type of BIC to calculate. Default is the mixed effects.

Value

(numeric) BIC of model fit

Examples

## Not run: 
model <- load_example_model("pheno")
results <- load_example_modelfit_results("pheno")
ofv <- results$ofv
calculate_bic(model, ofv)
calculate_bic(model, ofv, type='fixed')
calculate_bic(model, ofv, type='random')
calculate_bic(model, ofv, type='iiv')

## End(Not run)

pharmpy/pharmr documentation built on Feb. 15, 2025, 3:22 a.m.