tests/dev/test-timeprofile-groups.R

In Open-Systems-Pharmacology/TLF-Library: TLF Library

# The next line simply remove the variables from the environment
# Similar to clear from Matlab
rm(list = ls())

library(ggplot2)
library(tlf)

# Set your folder to test folder:
# setwd("C:/Design2Code/TLF-Playground/TLF-Playground/R/")
useTheme(tlfTheme)

# -------------------------------------------------
# Get the data and metadata for PK Ratio plot

nPopulation <- 20

# -------------------------------------------------

testData <- data.frame(
  IndivdualID = c(rep(1, 24), rep(2, 24)),
  Time = c(seq(1, 24), seq(1, 24)),
  Mean = c(10 * exp(-0.06 * seq(1, 24)), 10 * exp(-0.1 * seq(1, 24))),
  Min = c(5 * exp(-0.06 * seq(1, 24)), 8 * exp(-0.1 * seq(1, 24))),
  Max = c(15 * exp(-0.06 * seq(1, 24)), 12 * exp(-0.1 * seq(1, 24)))
)



testMetaData <- list(
  IndivdualID = list("unit" = "", "dimension" = ""),
  Time = list("unit" = "min", "dimension" = "Time"),
  Mean = list("unit" = "mg/L", "dimension" = "Concentration"),
  Min = list("unit" = "mg/L", "dimension" = "Concentration"),
  Max = list("unit" = "mg/L", "dimension" = "Concentration")
)

# -------------------------------------------------
# Define Default plot Configuration & Mapping from R6 class for PK Ratio
group <- GroupMapping$new(color = "IndivdualID")

meanDataMapping <- TimeProfileDataMapping$new(
  x = "Time",
  y = "Mean",
  groupMapping = group
)

errorDataMapping <- TimeProfileDataMapping$new(
  x = "Time",
  y = "Mean",
  ymin = "Min",
  ymax = "Max",
  groupMapping = group
)

group <- GroupMapping$new(fill = "IndivdualID")
rangeDataMapping <- TimeProfileDataMapping$new(
  x = "Time",
  ymin = "Min",
  ymax = "Max",
  groupMapping = group
)
# -------------------------------------------------
# Plot PK Ratio using the previously defined variables
meanPlot <- plotTimeProfile(data = testData, metaData = testMetaData, dataMapping = meanDataMapping)

errorPlot <- plotTimeProfile(data = testData, metaData = testMetaData, dataMapping = errorDataMapping)

rangePlot <- plotTimeProfile(data = testData, metaData = testMetaData, dataMapping = rangeDataMapping)