calculate_bic function

calculate_bic

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)

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.

Returns

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