run_vpc — run_vpc • nonmem.utils

Run a VPC analysis on the provided data.

run_vpc(data, x = "TIME", group = "All", bins = 7, ci = 0.8, lloq = 10)

Arguments

data: A data.frame of the pk data
x: Name of x variable, usually "TIME"
group: Name of group variable
bins: Number of bins for the VPC
ci: Confidence interval for the VPC, usually 0.8 (10th-90th percentiles)
lloq: Lower limit of quantification

Value

A data.frame

Examples


run_vpc(data_501, x = "TIME")
#> # A tibble: 4 × 82
#> # Groups:   Bins [4]
#>   Bins  time_min_med time_min_min time_min_max obs_min_med obs_min_min
#>   <fct>        <dbl>        <dbl>        <dbl>       <dbl>       <dbl>
#> 1 1                1            1            1        4.89        4.89
#> 2 5                5            5            5       19.1        19.1 
#> 3 12              12           12           12       12.5        12.5 
#> 4 24              24           24           24        3.28        3.28
#> # ℹ 76 more variables: obs_min_max <dbl>, dv_min_med <dbl>, dv_min_min <dbl>,
#> #   dv_min_max <dbl>, pc_obs_min_med <dbl>, pc_obs_min_min <dbl>,
#> #   pc_obs_min_max <dbl>, pc_dv_min_med <dbl>, pc_dv_min_min <dbl>,
#> #   pc_dv_min_max <dbl>, npde_obs_min_med <dbl>, npde_obs_min_min <dbl>,
#> #   npde_obs_min_max <dbl>, npde_dv_min_med <dbl>, npde_dv_min_min <dbl>,
#> #   npde_dv_min_max <dbl>, blq_obs_min_med <dbl>, blq_obs_min_min <dbl>,
#> #   blq_obs_min_max <dbl>, blq_dv_min_med <dbl>, blq_dv_min_min <dbl>, …