get_var_y {PKPDsim}

Get expected variance/sd/ci of dependent variable based on PKPDsim model, parameters, and regimen

Description

Get expected variance/sd/ci of dependent variable based on PKPDsim model, parameters, and regimen

Usage

get_var_y(
  model = NULL,
  parameters = list(),
  regimen = list(),
  t_obs = c(1:48),
  obs_comp = NULL,
  obs_variable = NULL,
  omega = c(0.1, 0.05, 0.1),
  omega_full = NULL,
  n_ind = NULL,
  ruv = NULL,
  y = NULL,
  rel_delta = 1e-04,
  method = "delta",
  sequence = NULL,
  auc = FALSE,
  sd = TRUE,
  q = NULL,
  in_parallel = FALSE,
  n_cores = 3,
  return_all = FALSE,
  ...
)

Arguments

model

model, created using 'PKPDsim::new_ode_model()'

parameters

parameters list

regimen

regimen, as created using 'PKPDsim::new_regimen()'

t_obs

vector of observation times

obs_comp

observation compartment. If NULL will be "obs" (default)

obs_variable

observation variable. If NULL, will be ignored, otherwise will override 'obs_comp'.

omega

triangle omega block

omega_full

full omega block

n_ind

number of individuals to simulate with sim method

ruv

residual variability, supplied as a named list, ex: 'list(prop = 0, add = 0, exp = 0)'

y

vector of observations. If NULL, then a new simulation will be performed.

rel_delta

rel_delta

method

method, 'delta' or 'sim'

sequence

for simulations, if not NULL the pseudo-random sequence to use, e.g. "halton" or "sobol". See 'mvrnorm2' for more details.

auc

is AUC?

sd

return as standard deviation ('TRUE') or variance ('FALSE')

q

return vector of quantiles instead of sd/var. Will return parametric quantiles when delta-method is used, non-parametric for simulation-based methods.

in_parallel

run simulations in parallel?

n_cores

if run in parallel, on how many cores?

return_all

return object with all relevant information?

...

passed on to 'sim_ode()'

Value

Vector of standard deviations or variances (or quantiles thereof) for dependent value variable