R/mean-model-workflow.R

In Open-Systems-Pharmacology/OSPSuite.ReportingEngine: R package to create reports for the Open Systems Pharmacology models

#' @title MeanModelWorkflow
#' @description  R6 class for Reporting Engine Mean Model Workflow
#' @field simulate `SimulationTask` object for time profile simulations
#' @field calculatePKParameters `CalculatePKParametersTask` object for PK parameters calculation
#' @field calculateSensitivity `SensitivityAnalysisTask` object for sensitivity analysis
#' @field plotTimeProfilesAndResiduals `PlotTask` object for goodness of fit plots
#' @field plotMassBalance `PlotTask` object for mass balance plot
#' @field plotAbsorption `PlotTask` object for absorption plot
#' @field plotPKParameters `PlotTask` object for PK parameters plot
#' @field plotSensitivity `PlotTask` object for sensitivity plot
#' @export
#' @import tlf
#' @family workflows
MeanModelWorkflow <- R6::R6Class(
  "MeanModelWorkflow",
  inherit = Workflow,
  public = list(
    simulate = NULL,
    calculatePKParameters = NULL,
    calculateSensitivity = NULL,
    plotTimeProfilesAndResiduals = NULL,
    plotMassBalance = NULL,
    plotAbsorption = NULL,
    plotPKParameters = NULL,
    plotSensitivity = NULL,

    #' @description
    #' Create a new `MeanModelWorkflow` object.
    #' @param ... input parameters inherited from R6 class object `Workflow`.
    #' @return A new `MeanModelWorkflow` object
    initialize = function(...) {
      super$initialize(...)
      logCatch({
        self$simulate <- loadSimulateTask(self)
        self$calculatePKParameters <- loadCalculatePKParametersTask(self)
        self$calculateSensitivity <- loadCalculateSensitivityTask(self)

        self$plotTimeProfilesAndResiduals <- loadPlotTimeProfilesAndResidualsTask(self)
        self$plotMassBalance <- loadPlotMassBalanceTask(self)
        self$plotAbsorption <- loadPlotAbsorptionTask(self)
        self$plotPKParameters <- loadPlotPKParametersTask(self)
        self$plotSensitivity <- loadPlotSensitivityTask(self)

        self$taskNames <- enum(self$getAllTasks())
      })
    },

    #' @description
    #' Run mean model workflow tasks for all simulation sets if tasks are activated
    #' The order of tasks is as follows:
    #' # 1) Run simulations
    #' # 2) Perform PK and sensitivity analyses
    #' # 3) Perform plot tasks
    #' ## 3.a) time profiles and residual plots
    #' ## 3.b) absorption plots
    #' ## 3.c) mass balance plots
    #' ## 3.d) PK and sensitivity analyses tables and plots
    #' # 4) Render report
    #' @return All results and plots as a structured output in the workflow folder
    runWorkflow = function() {
      # Prevent crashes if folder was deleted before (re) running a worflow
      dir.create(self$workflowFolder, showWarnings = FALSE, recursive = TRUE)
      # In case other logs were saved before running workflow
      setLogFolder(self$workflowFolder)
      actionToken1 <- re.tStartMetadataCapture(metaDataCapture = TRUE)
      actionToken2 <- re.tStartAction(actionType = "Run")

      logInfo(messages$runStarting("Mean Model Workflow"))
      t0 <- tic()

      logCatch({
        if (self$simulate$active) {
          self$simulate$runTask(self$simulationStructures)
        }

        if (self$calculatePKParameters$active) {
          self$calculatePKParameters$runTask(self$simulationStructures)
        }

        if (self$calculateSensitivity$active) {
          self$calculateSensitivity$runTask(self$simulationStructures)
        }

        for (plotTask in self$getAllPlotTasks()) {
          if (self[[plotTask]]$active) {
            self[[plotTask]]$runTask(self$simulationStructures)
          }
        }

        for (userDefinedTask in self$userDefinedTasks) {
          if (userDefinedTask$active) {
            userDefinedTask$runTask(self$simulationStructures)
          }
        }
        # Merge appendices into final report
        appendices <- c(
          as.character(sapply(self$getAllPlotTasks(), function(taskName) {
            self[[taskName]]$fileName
          })),
          as.character(sapply(self$userDefinedTasks, function(userDefinedTask) {
            userDefinedTask$fileName
          }))
        )
        appendices <- appendices[file.exists(appendices)]
        if (length(appendices) > 0) {
          initialReportPath <- file.path(self$workflowFolder, self$reportFileName)
          mergeMarkdownFiles(appendices, initialReportPath)
          renderReport(
            file.path(self$workflowFolder, self$reportFileName),
            createWordReport = self$createWordReport,
            numberSections = self$numberSections,
            intro = getIntroFromReportTitle(self$reportTitle),
            wordConversionTemplate = self$wordConversionTemplate
          )
          # Move report if a non-default path is provided
          copyReport(from = initialReportPath, to = self$reportFilePath, copyWordReport = self$createWordReport, keep = TRUE)
        }

        re.tStoreFileMetadata(access = "write", filePath = file.path(self$workflowFolder, defaultFileNames$logInfoFile()))
        re.tStoreFileMetadata(access = "write", filePath = file.path(self$workflowFolder, defaultFileNames$logDebugFile()))
        if (file.exists(file.path(self$workflowFolder, defaultFileNames$logErrorFile()))) {
          re.tStoreFileMetadata(access = "write", filePath = file.path(self$workflowFolder, defaultFileNames$logErrorFile()))
        }

        re.tEndAction(actionToken = actionToken2)
        re.tEndMetadataCapture(outputFolder = "./", actionToken = actionToken1)
      })
      logInfo(messages$runCompleted(getElapsedTime(t0), "Mean Model Workflow"))
      # Stop logging messages in workflowFolder after run is completed
      # Prevents potential logging of new messages in previous workflowFolder
      setLogFolder()
    }
  )
)