pkbuild: Automatic PK model building

In MarcLavielle/Rsmlx: R Speaks 'Monolix'

Automatic PK model building

Description

Fit several structural PK models and select the best one based on a Bayesian Information Criterion. Models to compare can be defined by rate constants and/or clearances and can include or not nonlinear elimination models. See https://monolix.lixoft.com/rsmlx/pkbuild/ for more details.

Usage

pkbuild(
  data = NULL,
  project = NULL,
  stat = FALSE,
  param = "clearance",
  new.dir = ".",
  MM = FALSE,
  linearization = F,
  criterion = "BICc",
  level = NULL,
  settings.stat = NULL
)

Arguments

data	a list with fields dataFile: path of a formatted data file headerTypes: a vector of strings administration ("iv", "bolus", "infusion", "oral", "ev"): route of administration
project	a Monolix project
stat	(FALSE, TRUE): the statistical model is also built (using buildmlx) (default=FALSE)
param	("clearance", "rate", "both): parametrization (default="clearance")
new.dir	name of the directory where the created files are stored (default is the current working directory) )
MM	(FALSE, TRUE): tested models include or not Michaelis Menten elimination models (default=FALSE)
linearization	TRUE/FALSE whether the computation of the likelihood is based on a linearization of the model (default=FALSE)
criterion	penalization criterion to optimize c("AIC", "BIC", "BICc", gamma) (default="BICc")
level	an integer between 1 and 9 (used by setSettings)
settings.stat	list of settings used by buildmlx (only if stat=TRUE)

Value

A list of results

Examples

## Not run: 
# Build a PK model for the warfarin PK data. 
# By default, only models using clearance (and inter compartmental clearances) are used
warf.pk1 <- pkbuild(data=warfarin)

# Models using elimination and transfer rate constants are used, 
# as well as nonlinear elimination models
warf.pk2 <- pkbuild(data=warfarin,  new.dir="warfarin", param="rate", MM=TRUE)

# Both models using clearances and rates are used. 
# Level is set to 7 in order to get accurate results.
warf.pk3 <- pkbuild(data=warfarin,  new.dir="warfarin", param="both", level=7)

## End(Not run)

MarcLavielle/Rsmlx documentation built on June 11, 2024, 8:22 a.m.