tests/testthat/helper_marginals.R

In openanalytics/BIGL: Biochemically Intuitive Generalized Loewe Model

library(BIGL)

data("directAntivirals", package = "BIGL")
data <- subset(directAntivirals, experiment == 1)[, c("effect", "d1", "d2")]

transforms <- list(
  "BiolT" = function(y, args) with(args, N0*exp(y*time.hours)),
  "InvBiolT" = function(T, args) with(args, 1/time.hours*log(T/N0)),
  "PowerT" = function(y, args) with(args, log(y)),
  "InvPowerT" = function(T, args) with(args, exp(T)),
  "compositeArgs" = list("N0" = 1,
                         "time.hours" = 72)
)

## Initial parameters for estimation without transformations
initPars <- initialMarginal(data)
initPars <- initPars + seq(0, 1, length.out = 7)

## Initial parameters for estimations with transformations
initParsT <- initialMarginal(data, transforms = transforms)
initParsT <- initParsT + c(-0.1, 0.2, 0.3, 0.4, 0.6, -0.4, 0.3)

names(initPars) <- names(initParsT) <- c("h1", "h2", "b", "m1", "m2", "e1", "e2")