calculate_bic {pharmr}

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_individals)

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

Usage

calculate_bic(model, likelihood, type = NULL)

Arguments

model

(Model) Pharmpy model object

likelihood

(numeric) -2LL to use

type

(str (optional)) 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)