get_var_y function

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

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()

Returns

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