Skip to contents

Carry out and visualize sensitivity analysis (with OSPSuite)

Source: R/sensitivity-calculation.R
sensitivityCalculation.Rd

Carry out and visualize sensitivity analysis (with OSPSuite)

Usage

sensitivityCalculation(
  simulation,
  outputPaths,
  parameterPaths,
  variationRange = c(seq(0.1, 1, by = 0.1), seq(2, 10, by = 1)),
  pkParameters = c("C_max", "t_max", "AUC_inf"),
  pkDataFilePath = NULL,
  simulationRunOptions = NULL
)

Arguments

simulation

An object of type Simulation.

outputPaths

Path (or a vector of paths) to the output(s) for which the sensitivity will be analyzed.

parameterPaths

A single or a vector of the parameter path(s) to be varied.

variationRange

Optional numeric vector defining the scaling of the parameters. The same variation range is applied to all specified parameters. If not specified, the following vector will be used: c(0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10).

pkParameters

A vector of names of PK parameters for which the sensitivities will be calculated. For a full set of available standard PK parameters, run names(ospsuite::StandardPKParameter). By default, only the following parameters will be considered: "C_max", "t_max", "AUC_inf". If NULL, all available PK-parameters will be calculated. You can also specify custom PK parameters.

pkDataFilePath

Path to excel file in which PK-parameter data should be saved. If a file already exists, it will be overwritten. Default is NULL, meaning the data will not be saved to a spreadsheet.

simulationRunOptions

Optional instance of a SimulationRunOptions used during the simulation run

Value

A list containing following objects:

  • SimulationResults

  • specified output paths

  • specified parameter paths

  • A data frame of PK parameters

See also

Other sensitivity-calculation: sensitivitySpiderPlot(), sensitivityTimeProfiles()

Examples

if (FALSE) {
simPath <- system.file("extdata", "Aciclovir.pkml", package = "ospsuite")
simulation <- loadSimulation(simPath)
outputPaths <- "Organism|PeripheralVenousBlood|Aciclovir|Plasma (Peripheral Venous Blood)"
parameterPaths <- c(
  "Aciclovir|Lipophilicity",
  "Applications|IV 250mg 10min|Application_1|ProtocolSchemaItem|Dose",
  "Neighborhoods|Kidney_pls_Kidney_ur|Aciclovir|Glomerular Filtration-GFR|GFR fraction"
)

# extract the results into a list of dataframes
sensitivityCalculation(
  simulation = simulation,
  outputPaths = outputPaths,
  parameterPaths = parameterPaths
)
}