Nothing
library(testthat)
context("Test the conversion to RxODE")
source(paste0("", "testUtils.R"))
test_that("Export method works (2-cpt model, comments)", {
regFilename = "2_cpt_model_abs_comments"
model <- model_suite$testing$nonmem$advan4_trans4
# Add some comments/line breaks programmatically
main <- model@model %>% getByName("MAIN")
main <- main %>% replace(Equation("KA", "THETA_KA*exp(ETA_KA)", comment="EQUATION COMMENT"))
main <- main %>% add(LineBreak())
main <- main %>% add(Comment("THIS IS A COMMENT"))
model <- model %>% replace(main)
model <- model %>% add(Bioavailability(compartment=1, rhs="0.75"))
model <- model %>% add(LagTime(compartment=2, rhs="2"))
model <- model %>% add(InfusionDuration(compartment=2, rhs="2"))
model <- model %>% add(InfusionRate(compartment=1, rhs="2"))
model <- model %>% add(InitialCondition(compartment=2, rhs="50"))
# Code comparison
rxmod <- model %>% export(dest="RxODE")
rxodeNonRegTest(rxmod, regFilename)
# Theta vector comparison
expect_equal(c(THETA_KA=1, THETA_CL=5, THETA_V2=80, THETA_V3=20, THETA_Q=4), rxmod@theta)
# Omega matrix comparison
omega <- diag(rep(0.025, 5))
rownames(omega) <- c("ETA_KA", "ETA_CL", "ETA_V2", "ETA_V3", "ETA_Q")
colnames(omega) <- c("ETA_KA", "ETA_CL", "ETA_V2", "ETA_V3", "ETA_Q")
expect_equal(omega, rxmod@omega)
# Sigma matrix comparison
sigma <- matrix(0.025)
rownames(sigma) <- c("EPS_PROP")
colnames(sigma) <- c("EPS_PROP")
expect_equal(sigma, rxmod@sigma)
})
test_that("Export method works (2-cpt model, if-statements)", {
regFilename = "2_cpt_model_abs_if_statements"
model <- model_suite$testing$nonmem$advan4_trans4
model <- model %>% delete(Equation("KA"))
model <- model %>% add(Equation("KA", "0"))
model <- model %>% add(IfStatement("OCC==1", Equation("KA", "THETA_KA*1.5*exp(ETA_KA)")))
model <- model %>% add(IfStatement("OCC==2", Equation("KA", "THETA_KA*0.5*exp(ETA_KA)")))
model <- model %>% add(IfStatement("OCC==3", Equation("KA", "THETA_KA*0.1*exp(ETA_KA)")))
# Code comparison
rxmod <- model %>% export(dest="RxODE")
rxodeNonRegTest(rxmod, regFilename)
})
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.