tests/testthat/test-theta-lhs-as.expression.R

In lotri: A Simple Way to Specify Symmetric, Block Diagonal Matrices

test_that("as.expression handling; lhs of theta parameters", {
  .lotri <- loadNamespace("lotri")
  
  x <- lotri({
    tka <- 0.45; label("Log Ka")
    tcl <- 1; label("Log Cl")
    tv <- 3.45; label("Log V")
    eta.ka ~ 0.6
    eta.cl ~ 0.3
    eta.v ~ 0.1
    add.err <- 0.7
  })
  
  xdf <- as.data.frame(x)
  x1 <- xdf[1, ]
  
  expect_equal(.lotri$.lotriLhsExpressionFromDf1(x1), quote(0.45))
  
  x1$fix <- TRUE
  expect_equal(.lotri$.lotriLhsExpressionFromDf1(x1),quote(fix(0.45)))
  
  x1$upper <- 1
  expect_equal(.lotri$.lotriLhsExpressionFromDf1(x1), quote(fix(-Inf, 0.45, 1)))
  
  x1$fix <- FALSE
  expect_equal(.lotri$.lotriLhsExpressionFromDf1(x1), quote(c(-Inf, 0.45, 1)))
  
  x1$upper <- Inf
  x1$lower <- 0
  expect_equal(.lotri$.lotriLhsExpressionFromDf1(x1), quote(c(0, 0.45)))
  
  x1$fix <- TRUE
  expect_equal(.lotri$.lotriLhsExpressionFromDf1(x1), quote(fix(0, 0.45)))
  
  x1$upper <- 1
  expect_equal(.lotri$.lotriLhsExpressionFromDf1(x1), quote(fix(0, 0.45, 1)))
  expect_equal(.lotri$.lotriAssignmentExpressionFromDf1(x1), quote(tka <- fix(0, 0.45, 1)))
  
  x1$fix <- FALSE
  expect_equal(.lotri$.lotriLhsExpressionFromDf1(x1), quote(c(0, 0.45, 1)))
  expect_equal(.lotri$.lotriAssignmentExpressionFromDf1(x1), quote(tka <- c(0, 0.45, 1)))
  
  expect_equal(.lotri$.lotriBackTransformFromDf1(x1), NULL)
  
  expect_equal(.lotriExpressionLinesFromDf1(x1),
               list(quote(tka <- c(0, 0.45, 1)),
                    quote(label("Log Ka"))))
  
  x1$backTransform <- "exp"
  expect_equal(.lotri$.lotriBackTransformFromDf1(x1), list(quote(backTransform("exp"))))
  
  expect_equal(.lotriExpressionLinesFromDf1(x1),
               list(quote(tka <- c(0, 0.45, 1)),
                    quote(backTransform("exp")),
                    quote(label("Log Ka"))))
  
  expect_equal(.lotri$.lotriLabelFromDf1(x1), list(quote(label("Log Ka"))))
  
  
  x1$label <- "This is a fun label\"'"
  expect_equal(.lotri$.lotriLabelFromDf1(x1), list(quote(label("This is a fun label\"'"))))
  
  x1$label <- NA_character_
  expect_equal(.lotri$.lotriLabelFromDf1(x1), NULL)
  
  expect_equal(.lotri$.lotriExpressionLinesFromDf1(x1),
               list(quote(tka <- c(0, 0.45, 1)),
                    quote(backTransform("exp"))))
  
  x1$backTransform <- NA_character_
  expect_equal(.lotri$.lotriExpressionLinesFromDf1(x1),
               list(quote(tka <- c(0, 0.45, 1))))
  
  expect_equal(.lotri$.lotriGetPopLinesFromDf(xdf),
               list(
                 quote(tka <- 0.45),
                 quote(label("Log Ka")),
                 quote(tcl <- 1),
                 quote(label("Log Cl")),
                 quote(tv <- 3.45),
                 quote(label("Log V")),
                 quote(add.err <- 0.7)))
  
  x1 <- lotri({
    et5 ~ 1
    et2 + et3 ~ c(
      1,
      2, 3
    )
    et1 ~ fix(3)
  })
  
  expect_equal(.lotri$.lotriGetEtaMatrixElements(x1),
               list(quote(et5 ~ 1),
                    quote(et2 + et3 ~ c(1, 2, 3)),
                    quote(et1 ~ fix(3))))
  
  fix2 <- lotri({
    f+g ~ fix(1,
              0.5, 1) | occ
    m+n ~ c(2,
            0.5, 1)
  })
  
  expect_equal(.lotri$.lotriGetEtaMatrixElements(fix2),
               list(quote(m + n ~ c(2, 0.5, 1)),
                    quote(f + g ~ fix(1, 0.5, 1) | occ)))
  
  x <- lotri({
    tka <- 0.45; label("Log Ka")
    tcl <- 1; label("Log Cl")
    tv <- 3.45; label("Log V")
    eta.ka ~ 0.6
    eta.cl ~ 0.3
    eta.v ~ 0.1
    add.err <- 0.7
  })
  
  df <- as.data.frame(x)
  
  expect_equal(lotriDataFrameToLotriExpression(df),
               quote(lotri({
                 tka <- 0.45
                 label("Log Ka")
                 tcl <- 1
                 label("Log Cl")
                 tv <- 3.45
                 label("Log V")
                 add.err <- 0.7
                 eta.ka ~ 0.6
                 eta.cl ~ 0.3
                 eta.v ~ 0.1
               })))
  
  
  expect_equal(lotriDataFrameToLotriExpression(df, useIni=TRUE),
               quote(ini({
                 tka <- 0.45
                 label("Log Ka")
                 tcl <- 1
                 label("Log Cl")
                 tv <- 3.45
                 label("Log V")
                 add.err <- 0.7
                 eta.ka ~ 0.6
                 eta.cl ~ 0.3
                 eta.v ~ 0.1
               })))
  
  # You may also call as.expression directly from the lotri object
  
  expect_equal(as.expression(x),
               quote(lotri({
                 tka <- 0.45
                 label("Log Ka")
                 tcl <- 1
                 label("Log Cl")
                 tv <- 3.45
                 label("Log V")
                 add.err <- 0.7
                 eta.ka ~ 0.6
                 eta.cl ~ 0.3
                 eta.v ~ 0.1
               })))
  
  expect_equal(as.expression(x, useIni=TRUE),
               quote(ini({
                 tka <- 0.45
                 label("Log Ka")
                 tcl <- 1
                 label("Log Cl")
                 tv <- 3.45
                 label("Log V")
                 add.err <- 0.7
                 eta.ka ~ 0.6
                 eta.cl ~ 0.3
                 eta.v ~ 0.1
               })))
  
  expect_error(lotriDataFrameToLotriExpression(rnorm))
  
})