VPC simulation

VPC simulation

Usage

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

Arguments

object

This is the nlmixr2 fit object

...

Other arguments sent to `rxSolve()`

keep

Column names to keep in the output simulated dataset

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

minN

With retries, the minimum number of studies to restimulate (by default 10)

normRelated

should the VPC style simulation be for normal related variables only

Value

data frame of the VPC simulation

Author

Matthew L. Fidler

Examples

# \donttest{

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")
#>  
#>  
#>  
#>  
#> ℹ parameter labels from comments are typically ignored in non-interactive mode
#> ℹ Need to run with the source intact to parse comments
#> → Calculating residuals/tables
#> ✔ done
#> → compress origData in nlmixr2 object, save 5952
#> → compress parHistData in nlmixr2 object, save 5944

head(vpcSim(fit, pred=TRUE))
#>  
#>  
#>   sim.id id time    ipred        sim  tad        depot  central nlmixrRowNums
#> 1      1  1 0.00 0.000000 -0.2599637 0.00 319.99200000   0.0000             2
#> 2      1  1 0.25 4.249815  3.8667868 0.25 170.36831700 148.3259             3
#> 3      1  1 0.57 6.836970  7.2090448 0.57  76.03060828 238.6220             4
#> 4      1  1 1.12 8.203210  7.2515453 1.12  18.99931015 286.3061             5
#> 5      1  1 2.02 8.219526  8.2583056 2.02   1.96443108 286.8756             6
#> 6      1  1 3.82 7.386852  7.6542381 3.82   0.02100075 257.8138             7
#>   rxLambda rxYj rxLow rxHi     pred
#> 1        1    2     0    1 0.000000
#> 2        1    2     0    1 3.263294
#> 3        1    2     0    1 5.832159
#> 4        1    2     0    1 7.870179
#> 5        1    2     0    1 8.513144
#> 6        1    2     0    1 7.628606

# }