In Open-Systems-Pharmacology/OSPSuite-R: R package to manipulate OSPSuite Models

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>",
  fig.showtext = TRUE
)

Unit conversion

Every entity - a molecule, a parameter, or an observer - has a certain dimension, like Amount, Concentration, or Volume. The dimension is a property of an entity:

library(ospsuite)

# Load a simulation
simFilePath <- system.file("extdata", "Aciclovir.pkml", package = "ospsuite")
sim <- loadSimulation(simFilePath)

# Get the parameter volume of the liver.
livParam <- getParameter("Organism|Liver|Volume", sim)
print(livParam)

# Dimension of the parameter
livParam$dimension

The values of a certain dimension can be presented in different units - for example, l or ml for the dimension Volume, or mol and µmol for the dimension Amount. The list of all available units for an entity can be obtained using allUnits() method:

# Dimension of the parameter
livParam$dimension

# Units of the parameter
livParam$allUnits

Internally, OSPS works with the base units, and all the values that are shown or passed to functions are in base units by default. These base units are often different from the units that are displayed by default in PK-Sim (and MoBi). The list of base and default display units can be found in the documentation.

As an example, the parameter BMI is given in the default unit kg/dm², while the default display unit is the more convenient kg/m².

The {ospuite} package provides a set of methods for conversion between different units. The methods toUnit(), toBaseUnit(), and toDisplayUnit() require the quantity to get the correct dimension and units; however, it does not change the value of the quantity!

# Get the BMI parameter
heightParam <- getParameter("Organism|Height", sim)
print(heightParam)

# Print the base and the default display units
heightParam$unit
heightParam$displayUnit

# Convert the value from the base into the default display unit
toDisplayUnit(quantity = heightParam, values = heightParam$value)

# Convert the value to the base unit, that can be used. e.g. for setting new parameter value
toBaseUnit(quantity = heightParam, values = 180, unit = "cm")

liverVolume <- getParameter("Organism|Liver|Volume", sim)
print(liverVolume)

liverVolume$allUnits

# Convert from base volume unit to µl
toUnit(quantity = liverVolume, values = c(1, 2, 3, 4), targetUnit = "ml")

Open-Systems-Pharmacology/OSPSuite-R documentation built on May 8, 2024, 11:36 a.m.