report-NO7.R

#' @title report-NO7
#' @description
#' Report from Test NO7 initial RE testing

# Clean environent and (re-)load the ospsuite.reportingengine package
rm(list=ls())
library(ospsuite.reportingengine)
ospsuite::removeAllUserDefinedPKParameters()
ospsuite::clearMemory(clearSimulationsCache = TRUE)

# Define inputs for workflow
simulationFile <- "Models/Theophylline (Boeckmann 1994) PO.pkml"
reportFolder <- "tests/Reports/Test-NO7"
reportTitle <- "Report NO7"

# Define data source
dataSource <- DataSource$new(
  dataFile = "Data/test-NO7/data.nmdat",
  metaDataFile = "Data/test-NO7/dict.csv",
  caption = "Boeckmann 1994"
)

# Define the simulation outputs possibly with their PK parameters
pkCmax <- PkParameterInfo$new(pkParameter = "C_max", displayUnit = "mg/l")
pkTmax <- PkParameterInfo$new(pkParameter = "t_max",  displayUnit = "h")
pkCtend <- PkParameterInfo$new(pkParameter = "C_tEnd",  displayUnit = "mg/l")
pkAuCtend <- PkParameterInfo$new(pkParameter = "AUC_tEnd",  displayUnit = "mg*h/l")

output <- Output$new(
  path = "Organism|PeripheralVenousBlood|Theophylline|Plasma (Peripheral Venous Blood)",
  displayName = "Theophylline",
  displayUnit = "mg/l",
  dataDisplayName = "Boeckmann 1994",
  dataSelection = "ALL",
  pkParameters = c(pkCmax, pkTmax, pkCtend, pkAuCtend)
)

# Define the simulation sets possibly mapping simulations to observed data
simulationSet <- SimulationSet$new(
  simulationFile = simulationFile,
  simulationSetName = "PO320mg",
  outputs = output,
  dataSource = dataSource
)

# Define the workflow and its result folder
workflow <- MeanModelWorkflow$new(
  simulationSets = simulationSet,
  workflowFolder = reportFolder,
  reportTitle = reportTitle,
  # To speed up the report generation, we skip the Word report here
  createWordReport = FALSE
)

# Define the workflow tasks
workflow$inactivateTasks()
workflow$activateTasks(tasks = c(
  "simulate",
  "calculatePKParameters",
  "plotTimeProfilesAndResiduals",
  "plotPKParameters"
  ))

# Run the workflow
workflow$runWorkflow()