vpcSim: VPC simulation

In nlmixr2est: Nonlinear Mixed Effects Models in Population PK/PD, Estimation Routines

VPC simulation

Description

VPC simulation

Usage

vpcSim(
  object,
  ...,
  keep = NULL,
  n = 300,
  pred = FALSE,
  seed = 1009,
  nretry = 50,
  normRelated = TRUE
)

Arguments

object	This is the nlmixr2 fit object
...	Other arguments sent to 'rxSolve()'
keep	Keep character vector
n	Number of simulations
pred	Should predictions be added to the simulation
seed	Seed to set for the VPC simulation
nretry	Number of times to retry the simulation if there is NA values in the simulation
normRelated	should the VPC style simulation be for normal related variables only

Value

data frame of the VPC simulation

Author(s)

Matthew L. Fidler

Examples

one.cmt <- function() {
 ini({
   ## You may label each parameter with a comment
   tka <- 0.45 # Log Ka
   tcl <- log(c(0, 2.7, 100)) # Log Cl
   ## This works with interactive models
   ## You may also label the preceding line with label("label text")
   tv <- 3.45; label("log V")
   ## the label("Label name") works with all models
   eta.ka ~ 0.6
   eta.cl ~ 0.3
   eta.v ~ 0.1
   add.sd <- 0.7
 })
 model({
   ka <- exp(tka + eta.ka)
   cl <- exp(tcl + eta.cl)
   v <- exp(tv + eta.v)
   linCmt() ~ add(add.sd)
 })
}

fit <- nlmixr(one.cmt, theo_sd, est="focei")

head(vpcSim(fit, pred=TRUE))

nlmixr2est documentation built on Oct. 8, 2023, 9:06 a.m.