Running 'testthat.R' [70s]
Running the tests in 'tests/testthat.R' failed.
Complete output:
  > Sys.setenv("R_TESTS" = "")
  > library(testthat)
  > library(campsis)
  Loading required package: campsismod
  
  Attaching package: 'campsismod'
  
  The following object is masked from 'package:utils':
  
      find
  
  The following objects are masked from 'package:base':
  
      Position, as.data.frame, replace, sort, toString, write
  
  
  Attaching package: 'campsis'
  
  The following object is masked from 'package:stats':
  
      simulate
  
  The following object is masked from 'package:base':
  
      sample
  
  > test_check("campsis")
  
  ---:: stderr ::---------------------------------------------
  
  ------------------------------------------------------------
  
  ---:: stderr ::---------------------------------------------
  
  ------------------------------------------------------------
  [ FAIL 39 | WARN 0 | SKIP 29 | PASS 389 ]
  
  ══ Skipped tests (29) ══════════════════════════════════════════════════════════
  • empty test (29): 'testDefaultPlots.R:8:1', 'testDefaultPlots.R:53:1',
    'testDefaultPlots.R:80:1', 'testSimulateEvents.R:8:1',
    'testSimulateEvents.R:32:1', 'testSimulateEvents.R:55:1',
    'testSimulateEvents.R:79:1', 'testSimulateEvents.R:112:1',
    'testSimulateEventsLogic.R:8:1', 'testSimulateEventsLogic.R:41:1',
    'testSimulateEventsLogic.R:79:1', 'testSimulateEventsLogic.R:125:1',
    'testSimulateEventsLogic.R:148:1', 'testSimulateIOV.R:8:1',
    'testSimulateIOV.R:47:1', 'testSimulateIOV.R:79:1',
    'testSimulateIOV.R:129:1', 'testSimulateIOV.R:193:1',
    'testSimulateIOV.R:241:1', 'testSimulateInfusion.R:8:1',
    'testSimulateNocbLocf.R:8:1', 'testSimulateNocbLocf.R:49:1',
    'testSimulateNocbLocf.R:89:1', 'testSimulateReplicates.R:9:1',
    'testSimulateReplicates.R:26:1', 'testSimulateReplicates.R:50:1',
    'testSimulateReplicates.R:71:1', 'testSimulateReplicates.R:107:1',
    'testSimulateTTEModels.R:7:1'
  
  ══ Failed tests ════════════════════════════════════════════════════════════════
  ── Failure ('testCompilationError.R:19:7'): Incorrect Campsis model does not compile and a clear error message is shown (rxode2/mrgsolve) ──
  `simulate(model = model, dataset = ds, dest = destEngine)` threw an error with unexpected message.
  Expected match: "The following parameter\\(s\\) are required for solving: BW"
  Actual message: "ℹ In index: 1.\nCaused by error:\nℹ In index: 1.\nCaused by error:\n! .onLoad failed in loadNamespace() for 'rxode2', details:\n  call: NULL\n  error: rxode2 compiled with rxode2parse with a different solving structure\ncan try: install.packages('rxode2', type='source')"
  ── Error ('testSimulateArguments.R:21:3'): Argument dest works well ────────────
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "KA", rhs = "THETA_KA*exp(ETA_KA)", 
              comment = NA_character_), new("equation", lhs = "CL", 
              rhs = "THETA_CL*exp(ETA_CL)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "V3", 
                  rhs = "THETA_V3*exp(ETA_V3)", comment = NA_character_), 
              new("equation", lhs = "Q", rhs = "THETA_Q*exp(ETA_Q)", 
                  comment = NA_character_), new("equation", lhs = "S2", 
                  rhs = "V2", comment = NA_character_)), type = "model_statement")), 
          new("ode_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("ode", lhs = "A_DEPOT", rhs = "-KA*A_DEPOT", 
                  comment = NA_character_), new("ode", lhs = "A_CENTRAL", 
                  rhs = "KA*A_DEPOT + Q*A_PERIPHERAL/V3 + (-CL/V2 - Q/V2)*A_CENTRAL", 
                  comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
                  rhs = "-Q*A_PERIPHERAL/V3 + Q*A_CENTRAL/V2", 
                  comment = NA_character_), new("ode", lhs = "A_OUTPUT", 
                  rhs = "CL*A_CENTRAL/V2", comment = NA_character_), 
                  new("equation", lhs = "F", rhs = "A_CENTRAL/S2", 
                    comment = NA_character_)), type = "model_statement")), 
          new("error_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("equation", lhs = "CP", rhs = "F", 
                  comment = NA_character_), new("equation", lhs = "OBS_CP", 
                  rhs = "CP*(EPS_PROP + 1)", comment = NA_character_), 
                  new("equation", lhs = "Y", rhs = "OBS_CP", comment = NA_character_)), 
              type = "model_statement"))), type = "code_record"), 
      parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
      0L)), list = list(new("theta", unit = NA_character_, name = "KA", 
          index = 1L, value = 1, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "CL", index = 2L, value = 5, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V2", index = 3L, value = 80, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V3", index = 4L, value = 20, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "Q", index = 5L, value = 4, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 1L, 
          type = "var", name = "KA", index = 1L, value = 0.025, 
          fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 2L, type = "var", name = "CL", 
              index = 2L, value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 3L, type = "var", name = "V2", index = 3L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 4L, type = "var", name = "V3", index = 4L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 5L, type = "var", name = "Q", index = 5L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("sigma", index2 = 1L, 
              type = "var", name = "PROP", index = 1L, value = 0.025, 
              fix = FALSE, label = NA_character_, comment = NA_character_)), 
          type = "parameter"), compartments = new("compartments", 
          properties = new("compartment_properties", list = list(), 
              type = "compartment_property"), list = list(new("compartment", 
              name = "DEPOT", index = 1L), new("compartment", name = "CENTRAL", 
              index = 2L), new("compartment", name = "PERIPHERAL", 
              index = 3L), new("compartment", name = "OUTPUT", 
              index = 4L)), type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateArguments.R:35:3'): Auto seed value vs fix seed + default engine ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "KA", rhs = "THETA_KA*exp(ETA_KA)", 
              comment = NA_character_), new("equation", lhs = "CL", 
              rhs = "THETA_CL*exp(ETA_CL)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "V3", 
                  rhs = "THETA_V3*exp(ETA_V3)", comment = NA_character_), 
              new("equation", lhs = "Q", rhs = "THETA_Q*exp(ETA_Q)", 
                  comment = NA_character_), new("equation", lhs = "S2", 
                  rhs = "V2", comment = NA_character_)), type = "model_statement")), 
          new("ode_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("ode", lhs = "A_DEPOT", rhs = "-KA*A_DEPOT", 
                  comment = NA_character_), new("ode", lhs = "A_CENTRAL", 
                  rhs = "KA*A_DEPOT + Q*A_PERIPHERAL/V3 + (-CL/V2 - Q/V2)*A_CENTRAL", 
                  comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
                  rhs = "-Q*A_PERIPHERAL/V3 + Q*A_CENTRAL/V2", 
                  comment = NA_character_), new("ode", lhs = "A_OUTPUT", 
                  rhs = "CL*A_CENTRAL/V2", comment = NA_character_), 
                  new("equation", lhs = "F", rhs = "A_CENTRAL/S2", 
                    comment = NA_character_)), type = "model_statement")), 
          new("error_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("equation", lhs = "CP", rhs = "F", 
                  comment = NA_character_), new("equation", lhs = "OBS_CP", 
                  rhs = "CP*(EPS_PROP + 1)", comment = NA_character_), 
                  new("equation", lhs = "Y", rhs = "OBS_CP", comment = NA_character_)), 
              type = "model_statement"))), type = "code_record"), 
      parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
      0L)), list = list(new("theta", unit = NA_character_, name = "KA", 
          index = 1L, value = 1, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "CL", index = 2L, value = 5, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V2", index = 3L, value = 80, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V3", index = 4L, value = 20, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "Q", index = 5L, value = 4, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 1L, 
          type = "var", name = "KA", index = 1L, value = 0.025, 
          fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 2L, type = "var", name = "CL", 
              index = 2L, value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 3L, type = "var", name = "V2", index = 3L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 4L, type = "var", name = "V3", index = 4L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 5L, type = "var", name = "Q", index = 5L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("sigma", index2 = 1L, 
              type = "var", name = "PROP", index = 1L, value = 0.025, 
              fix = FALSE, label = NA_character_, comment = NA_character_)), 
          type = "parameter"), compartments = new("compartments", 
          properties = new("compartment_properties", list = list(), 
              type = "compartment_property"), list = list(new("compartment", 
              name = "DEPOT", index = 1L), new("compartment", name = "CENTRAL", 
              index = 2L), new("compartment", name = "PERIPHERAL", 
              index = 3L), new("compartment", name = "OUTPUT", 
              index = 4L)), type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateBioavailability.R:30:3'): Simulate a bolus, 2 arms, F1 only in arm1, in dataset (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "KA", rhs = "THETA_KA*exp(ETA_KA)", 
              comment = NA_character_), new("equation", lhs = "CL", 
              rhs = "THETA_CL*exp(ETA_CL)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "V3", 
                  rhs = "THETA_V3*exp(ETA_V3)", comment = NA_character_), 
              new("equation", lhs = "Q", rhs = "THETA_Q*exp(ETA_Q)", 
                  comment = NA_character_), new("equation", lhs = "S2", 
                  rhs = "V2", comment = NA_character_)), type = "model_statement")), 
          new("ode_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("ode", lhs = "A_DEPOT", rhs = "-KA*A_DEPOT", 
                  comment = NA_character_), new("ode", lhs = "A_CENTRAL", 
                  rhs = "KA*A_DEPOT + Q*A_PERIPHERAL/V3 + (-CL/V2 - Q/V2)*A_CENTRAL", 
                  comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
                  rhs = "-Q*A_PERIPHERAL/V3 + Q*A_CENTRAL/V2", 
                  comment = NA_character_), new("ode", lhs = "A_OUTPUT", 
                  rhs = "CL*A_CENTRAL/V2", comment = NA_character_), 
                  new("equation", lhs = "F", rhs = "A_CENTRAL/S2", 
                    comment = NA_character_)), type = "model_statement")), 
          new("error_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("equation", lhs = "CP", rhs = "F", 
                  comment = NA_character_), new("equation", lhs = "OBS_CP", 
                  rhs = "CP*(EPS_PROP + 1)", comment = NA_character_), 
                  new("equation", lhs = "Y", rhs = "OBS_CP", comment = NA_character_)), 
              type = "model_statement"))), type = "code_record"), 
      parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
      0L)), list = list(new("theta", unit = NA_character_, name = "KA", 
          index = 1L, value = 1, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "CL", index = 2L, value = 5, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V2", index = 3L, value = 80, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V3", index = 4L, value = 20, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "Q", index = 5L, value = 4, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 1L, 
          type = "var", name = "KA", index = 1L, value = 0.025, 
          fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 2L, type = "var", name = "CL", 
              index = 2L, value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 3L, type = "var", name = "V2", index = 3L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 4L, type = "var", name = "V3", index = 4L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 5L, type = "var", name = "Q", index = 5L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("sigma", index2 = 1L, 
              type = "var", name = "PROP", index = 1L, value = 0.025, 
              fix = FALSE, label = NA_character_, comment = NA_character_)), 
          type = "parameter"), compartments = new("compartments", 
          properties = new("compartment_properties", list = list(), 
              type = "compartment_property"), list = list(new("compartment", 
              name = "DEPOT", index = 1L), new("compartment", name = "CENTRAL", 
              index = 2L), new("compartment", name = "PERIPHERAL", 
              index = 3L), new("compartment", name = "OUTPUT", 
              index = 4L)), type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateBioavailability.R:49:3'): Simulate a simple bolus with bioavailability, dataset versus model (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "KA", rhs = "THETA_KA*exp(ETA_KA)", 
              comment = NA_character_), new("equation", lhs = "CL", 
              rhs = "THETA_CL*exp(ETA_CL)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "V3", 
                  rhs = "THETA_V3*exp(ETA_V3)", comment = NA_character_), 
              new("equation", lhs = "Q", rhs = "THETA_Q*exp(ETA_Q)", 
                  comment = NA_character_), new("equation", lhs = "S2", 
                  rhs = "V2", comment = NA_character_)), type = "model_statement")), 
          new("ode_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("ode", lhs = "A_DEPOT", rhs = "-KA*A_DEPOT", 
                  comment = NA_character_), new("ode", lhs = "A_CENTRAL", 
                  rhs = "KA*A_DEPOT + Q*A_PERIPHERAL/V3 + (-CL/V2 - Q/V2)*A_CENTRAL", 
                  comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
                  rhs = "-Q*A_PERIPHERAL/V3 + Q*A_CENTRAL/V2", 
                  comment = NA_character_), new("ode", lhs = "A_OUTPUT", 
                  rhs = "CL*A_CENTRAL/V2", comment = NA_character_), 
                  new("equation", lhs = "F", rhs = "A_CENTRAL/S2", 
                    comment = NA_character_)), type = "model_statement")), 
          new("error_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("equation", lhs = "CP", rhs = "F", 
                  comment = NA_character_), new("equation", lhs = "OBS_CP", 
                  rhs = "CP*(EPS_PROP + 1)", comment = NA_character_), 
                  new("equation", lhs = "Y", rhs = "OBS_CP", comment = NA_character_)), 
              type = "model_statement"))), type = "code_record"), 
      parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
      0L)), list = list(new("theta", unit = NA_character_, name = "KA", 
          index = 1L, value = 1, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "CL", index = 2L, value = 5, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V2", index = 3L, value = 80, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V3", index = 4L, value = 20, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "Q", index = 5L, value = 4, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 1L, 
          type = "var", name = "KA", index = 1L, value = 0.025, 
          fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 2L, type = "var", name = "CL", 
              index = 2L, value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 3L, type = "var", name = "V2", index = 3L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 4L, type = "var", name = "V3", index = 4L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 5L, type = "var", name = "Q", index = 5L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("sigma", index2 = 1L, 
              type = "var", name = "PROP", index = 1L, value = 0.025, 
              fix = FALSE, label = NA_character_, comment = NA_character_)), 
          type = "parameter"), compartments = new("compartments", 
          properties = new("compartment_properties", list = list(), 
              type = "compartment_property"), list = list(new("compartment", 
              name = "DEPOT", index = 1L), new("compartment", name = "CENTRAL", 
              index = 2L), new("compartment", name = "PERIPHERAL", 
              index = 3L), new("compartment", name = "OUTPUT", 
              index = 4L)), type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateBioavailability.R:80:3'): Simulate several fixed F's at once (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "KA", rhs = "THETA_KA*exp(ETA_KA)", 
              comment = NA_character_), new("equation", lhs = "CL", 
              rhs = "THETA_CL*exp(ETA_CL)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "V3", 
                  rhs = "THETA_V3*exp(ETA_V3)", comment = NA_character_), 
              new("equation", lhs = "Q", rhs = "THETA_Q*exp(ETA_Q)", 
                  comment = NA_character_), new("equation", lhs = "S2", 
                  rhs = "V2", comment = NA_character_)), type = "model_statement")), 
          new("ode_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("ode", lhs = "A_DEPOT", rhs = "-KA*A_DEPOT", 
                  comment = NA_character_), new("ode", lhs = "A_CENTRAL", 
                  rhs = "KA*A_DEPOT + Q*A_PERIPHERAL/V3 + (-CL/V2 - Q/V2)*A_CENTRAL", 
                  comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
                  rhs = "-Q*A_PERIPHERAL/V3 + Q*A_CENTRAL/V2", 
                  comment = NA_character_), new("ode", lhs = "A_OUTPUT", 
                  rhs = "CL*A_CENTRAL/V2", comment = NA_character_), 
                  new("equation", lhs = "F", rhs = "A_CENTRAL/S2", 
                    comment = NA_character_)), type = "model_statement")), 
          new("error_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("equation", lhs = "CP", rhs = "F", 
                  comment = NA_character_), new("equation", lhs = "OBS_CP", 
                  rhs = "CP*(EPS_PROP + 1)", comment = NA_character_), 
                  new("equation", lhs = "Y", rhs = "OBS_CP", comment = NA_character_)), 
              type = "model_statement"))), type = "code_record"), 
      parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
      0L)), list = list(new("theta", unit = NA_character_, name = "KA", 
          index = 1L, value = 1, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "CL", index = 2L, value = 5, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V2", index = 3L, value = 80, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V3", index = 4L, value = 20, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "Q", index = 5L, value = 4, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 1L, 
          type = "var", name = "KA", index = 1L, value = 0, fix = FALSE, 
          label = NA_character_, comment = NA_character_), new("omega", 
          same = NA, index2 = 2L, type = "var", name = "CL", index = 2L, 
          value = 0, fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 3L, type = "var", name = "V2", 
              index = 3L, value = 0, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 4L, type = "var", name = "V3", index = 4L, 
              value = 0, fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 5L, type = "var", name = "Q", 
              index = 5L, value = 0, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("sigma", index2 = 1L, 
              type = "var", name = "PROP", index = 1L, value = 0.025, 
              fix = FALSE, label = NA_character_, comment = NA_character_)), 
          type = "parameter"), compartments = new("compartments", 
          properties = new("compartment_properties", list = list(), 
              type = "compartment_property"), list = list(new("compartment", 
              name = "DEPOT", index = 1L), new("compartment", name = "CENTRAL", 
              index = 2L), new("compartment", name = "PERIPHERAL", 
              index = 3L), new("compartment", name = "OUTPUT", 
              index = 4L)), type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateBolus.R:23:3'): Simulate a bolus (rxode2/mrgsolve) ──────
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "KA", rhs = "THETA_KA*exp(ETA_KA)", 
              comment = NA_character_), new("equation", lhs = "CL", 
              rhs = "THETA_CL*exp(ETA_CL)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "V3", 
                  rhs = "THETA_V3*exp(ETA_V3)", comment = NA_character_), 
              new("equation", lhs = "Q", rhs = "THETA_Q*exp(ETA_Q)", 
                  comment = NA_character_), new("equation", lhs = "S2", 
                  rhs = "V2", comment = NA_character_)), type = "model_statement")), 
          new("ode_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("ode", lhs = "A_DEPOT", rhs = "-KA*A_DEPOT", 
                  comment = NA_character_), new("ode", lhs = "A_CENTRAL", 
                  rhs = "KA*A_DEPOT + Q*A_PERIPHERAL/V3 + (-CL/V2 - Q/V2)*A_CENTRAL", 
                  comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
                  rhs = "-Q*A_PERIPHERAL/V3 + Q*A_CENTRAL/V2", 
                  comment = NA_character_), new("ode", lhs = "A_OUTPUT", 
                  rhs = "CL*A_CENTRAL/V2", comment = NA_character_), 
                  new("equation", lhs = "F", rhs = "A_CENTRAL/S2", 
                    comment = NA_character_)), type = "model_statement")), 
          new("error_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("equation", lhs = "CP", rhs = "F", 
                  comment = NA_character_), new("equation", lhs = "OBS_CP", 
                  rhs = "CP*(EPS_PROP + 1)", comment = NA_character_), 
                  new("equation", lhs = "Y", rhs = "OBS_CP", comment = NA_character_)), 
              type = "model_statement"))), type = "code_record"), 
      parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
      0L)), list = list(new("theta", unit = NA_character_, name = "KA", 
          index = 1L, value = 1, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "CL", index = 2L, value = 5, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V2", index = 3L, value = 80, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V3", index = 4L, value = 20, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "Q", index = 5L, value = 4, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 1L, 
          type = "var", name = "KA", index = 1L, value = 0.025, 
          fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 2L, type = "var", name = "CL", 
              index = 2L, value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 3L, type = "var", name = "V2", index = 3L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 4L, type = "var", name = "V3", index = 4L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 5L, type = "var", name = "Q", index = 5L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("sigma", index2 = 1L, 
              type = "var", name = "PROP", index = 1L, value = 0.025, 
              fix = FALSE, label = NA_character_, comment = NA_character_)), 
          type = "parameter"), compartments = new("compartments", 
          properties = new("compartment_properties", list = list(), 
              type = "compartment_property"), list = list(new("compartment", 
              name = "DEPOT", index = 1L), new("compartment", name = "CENTRAL", 
              index = 2L), new("compartment", name = "PERIPHERAL", 
              index = 3L), new("compartment", name = "OUTPUT", 
              index = 4L)), type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateBolus.R:49:3'): Simulate a bolus, 2 arms (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "KA", rhs = "THETA_KA*exp(ETA_KA)", 
              comment = NA_character_), new("equation", lhs = "CL", 
              rhs = "THETA_CL*exp(ETA_CL)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "V3", 
                  rhs = "THETA_V3*exp(ETA_V3)", comment = NA_character_), 
              new("equation", lhs = "Q", rhs = "THETA_Q*exp(ETA_Q)", 
                  comment = NA_character_), new("equation", lhs = "S2", 
                  rhs = "V2", comment = NA_character_)), type = "model_statement")), 
          new("ode_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("ode", lhs = "A_DEPOT", rhs = "-KA*A_DEPOT", 
                  comment = NA_character_), new("ode", lhs = "A_CENTRAL", 
                  rhs = "KA*A_DEPOT + Q*A_PERIPHERAL/V3 + (-CL/V2 - Q/V2)*A_CENTRAL", 
                  comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
                  rhs = "-Q*A_PERIPHERAL/V3 + Q*A_CENTRAL/V2", 
                  comment = NA_character_), new("ode", lhs = "A_OUTPUT", 
                  rhs = "CL*A_CENTRAL/V2", comment = NA_character_), 
                  new("equation", lhs = "F", rhs = "A_CENTRAL/S2", 
                    comment = NA_character_)), type = "model_statement")), 
          new("error_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("equation", lhs = "CP", rhs = "F", 
                  comment = NA_character_), new("equation", lhs = "OBS_CP", 
                  rhs = "CP*(EPS_PROP + 1)", comment = NA_character_), 
                  new("equation", lhs = "Y", rhs = "OBS_CP", comment = NA_character_)), 
              type = "model_statement"))), type = "code_record"), 
      parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
      0L)), list = list(new("theta", unit = NA_character_, name = "KA", 
          index = 1L, value = 1, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "CL", index = 2L, value = 5, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V2", index = 3L, value = 80, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V3", index = 4L, value = 20, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "Q", index = 5L, value = 4, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 1L, 
          type = "var", name = "KA", index = 1L, value = 0.025, 
          fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 2L, type = "var", name = "CL", 
              index = 2L, value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 3L, type = "var", name = "V2", index = 3L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 4L, type = "var", name = "V3", index = 4L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 5L, type = "var", name = "Q", index = 5L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("sigma", index2 = 1L, 
              type = "var", name = "PROP", index = 1L, value = 0.025, 
              fix = FALSE, label = NA_character_, comment = NA_character_)), 
          type = "parameter"), compartments = new("compartments", 
          properties = new("compartment_properties", list = list(), 
              type = "compartment_property"), list = list(new("compartment", 
              name = "DEPOT", index = 1L), new("compartment", name = "CENTRAL", 
              index = 2L), new("compartment", name = "PERIPHERAL", 
              index = 3L), new("compartment", name = "OUTPUT", 
              index = 4L)), type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateBolus.R:72:3'): Simulate a bolus, 2 labelled arms (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "KA", rhs = "THETA_KA*exp(ETA_KA)", 
              comment = NA_character_), new("equation", lhs = "CL", 
              rhs = "THETA_CL*exp(ETA_CL)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "V3", 
                  rhs = "THETA_V3*exp(ETA_V3)", comment = NA_character_), 
              new("equation", lhs = "Q", rhs = "THETA_Q*exp(ETA_Q)", 
                  comment = NA_character_), new("equation", lhs = "S2", 
                  rhs = "V2", comment = NA_character_)), type = "model_statement")), 
          new("ode_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("ode", lhs = "A_DEPOT", rhs = "-KA*A_DEPOT", 
                  comment = NA_character_), new("ode", lhs = "A_CENTRAL", 
                  rhs = "KA*A_DEPOT + Q*A_PERIPHERAL/V3 + (-CL/V2 - Q/V2)*A_CENTRAL", 
                  comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
                  rhs = "-Q*A_PERIPHERAL/V3 + Q*A_CENTRAL/V2", 
                  comment = NA_character_), new("ode", lhs = "A_OUTPUT", 
                  rhs = "CL*A_CENTRAL/V2", comment = NA_character_), 
                  new("equation", lhs = "F", rhs = "A_CENTRAL/S2", 
                    comment = NA_character_)), type = "model_statement")), 
          new("error_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("equation", lhs = "CP", rhs = "F", 
                  comment = NA_character_), new("equation", lhs = "OBS_CP", 
                  rhs = "CP*(EPS_PROP + 1)", comment = NA_character_), 
                  new("equation", lhs = "Y", rhs = "OBS_CP", comment = NA_character_)), 
              type = "model_statement"))), type = "code_record"), 
      parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
      0L)), list = list(new("theta", unit = NA_character_, name = "KA", 
          index = 1L, value = 1, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "CL", index = 2L, value = 5, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V2", index = 3L, value = 80, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V3", index = 4L, value = 20, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "Q", index = 5L, value = 4, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 1L, 
          type = "var", name = "KA", index = 1L, value = 0.025, 
          fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 2L, type = "var", name = "CL", 
              index = 2L, value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 3L, type = "var", name = "V2", index = 3L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 4L, type = "var", name = "V3", index = 4L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 5L, type = "var", name = "Q", index = 5L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("sigma", index2 = 1L, 
              type = "var", name = "PROP", index = 1L, value = 0.025, 
              fix = FALSE, label = NA_character_, comment = NA_character_)), 
          type = "parameter"), compartments = new("compartments", 
          properties = new("compartment_properties", list = list(), 
              type = "compartment_property"), list = list(new("compartment", 
              name = "DEPOT", index = 1L), new("compartment", name = "CENTRAL", 
              index = 2L), new("compartment", name = "PERIPHERAL", 
              index = 3L), new("compartment", name = "OUTPUT", 
              index = 4L)), type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateBolusInfusionSameCmt.R:27:3'): Bolus and infusion in CMT 1 (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "KA", rhs = "THETA_KA*exp(ETA_KA)", 
              comment = NA_character_), new("equation", lhs = "CL", 
              rhs = "THETA_CL*exp(ETA_CL)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "V3", 
                  rhs = "THETA_V3*exp(ETA_V3)", comment = NA_character_), 
              new("equation", lhs = "Q", rhs = "THETA_Q*exp(ETA_Q)", 
                  comment = NA_character_), new("equation", lhs = "S2", 
                  rhs = "V2", comment = NA_character_)), type = "model_statement")), 
          new("ode_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("ode", lhs = "A_DEPOT", rhs = "-KA*A_DEPOT", 
                  comment = NA_character_), new("ode", lhs = "A_CENTRAL", 
                  rhs = "KA*A_DEPOT + Q*A_PERIPHERAL/V3 + (-CL/V2 - Q/V2)*A_CENTRAL", 
                  comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
                  rhs = "-Q*A_PERIPHERAL/V3 + Q*A_CENTRAL/V2", 
                  comment = NA_character_), new("ode", lhs = "A_OUTPUT", 
                  rhs = "CL*A_CENTRAL/V2", comment = NA_character_), 
                  new("equation", lhs = "F", rhs = "A_CENTRAL/S2", 
                    comment = NA_character_)), type = "model_statement")), 
          new("error_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("equation", lhs = "CP", rhs = "F", 
                  comment = NA_character_), new("equation", lhs = "OBS_CP", 
                  rhs = "CP*(EPS_PROP + 1)", comment = NA_character_), 
                  new("equation", lhs = "Y", rhs = "OBS_CP", comment = NA_character_)), 
              type = "model_statement"))), type = "code_record"), 
      parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
      0L)), list = list(new("theta", unit = NA_character_, name = "KA", 
          index = 1L, value = 1, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "CL", index = 2L, value = 5, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V2", index = 3L, value = 80, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V3", index = 4L, value = 20, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "Q", index = 5L, value = 4, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 1L, 
          type = "var", name = "KA", index = 1L, value = 0.025, 
          fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 2L, type = "var", name = "CL", 
              index = 2L, value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 3L, type = "var", name = "V2", index = 3L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 4L, type = "var", name = "V3", index = 4L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 5L, type = "var", name = "Q", index = 5L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("sigma", index2 = 1L, 
              type = "var", name = "PROP", index = 1L, value = 0.025, 
              fix = FALSE, label = NA_character_, comment = NA_character_)), 
          type = "parameter"), compartments = new("compartments", 
          properties = new("compartment_properties", list = list(
              new("compartment_infusion_duration", compartment = 1L, 
                  rhs = "5", comment = NA_character_)), type = "compartment_property"), 
          list = list(new("compartment", name = "DEPOT", index = 1L), 
              new("compartment", name = "CENTRAL", index = 2L), 
              new("compartment", name = "PERIPHERAL", index = 3L), 
              new("compartment", name = "OUTPUT", index = 4L)), 
          type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateDeclare.R:24:7'): Declare argument with mrgsolve ────────
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "KA", rhs = "THETA_KA*exp(ETA_KA + SOMETHING)", 
              comment = NA_character_), new("equation", lhs = "CL", 
              rhs = "THETA_CL*exp(ETA_CL)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "V3", 
                  rhs = "THETA_V3*exp(ETA_V3)", comment = NA_character_), 
              new("equation", lhs = "Q", rhs = "THETA_Q*exp(ETA_Q)", 
                  comment = NA_character_), new("equation", lhs = "S2", 
                  rhs = "V2", comment = NA_character_)), type = "model_statement")), 
          new("ode_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("ode", lhs = "A_DEPOT", rhs = "-KA*A_DEPOT", 
                  comment = NA_character_), new("ode", lhs = "A_CENTRAL", 
                  rhs = "KA*A_DEPOT + Q*A_PERIPHERAL/V3 + (-CL/V2 - Q/V2)*A_CENTRAL", 
                  comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
                  rhs = "-Q*A_PERIPHERAL/V3 + Q*A_CENTRAL/V2", 
                  comment = NA_character_), new("ode", lhs = "A_OUTPUT", 
                  rhs = "CL*A_CENTRAL/V2", comment = NA_character_), 
                  new("equation", lhs = "F", rhs = "A_CENTRAL/S2", 
                    comment = NA_character_)), type = "model_statement")), 
          new("error_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("equation", lhs = "CP", rhs = "F", 
                  comment = NA_character_), new("equation", lhs = "OBS_CP", 
                  rhs = "CP*(EPS_PROP + 1)", comment = NA_character_), 
                  new("equation", lhs = "Y", rhs = "OBS_CP", comment = NA_character_)), 
              type = "model_statement"))), type = "code_record"), 
      parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
      0L)), list = list(new("theta", unit = NA_character_, name = "KA", 
          index = 1L, value = 1, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "CL", index = 2L, value = 5, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V2", index = 3L, value = 80, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V3", index = 4L, value = 20, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "Q", index = 5L, value = 4, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 1L, 
          type = "var", name = "KA", index = 1L, value = 0.025, 
          fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 2L, type = "var", name = "CL", 
              index = 2L, value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 3L, type = "var", name = "V2", index = 3L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 4L, type = "var", name = "V3", index = 4L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 5L, type = "var", name = "Q", index = 5L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("sigma", index2 = 1L, 
              type = "var", name = "PROP", index = 1L, value = 0.025, 
              fix = FALSE, label = NA_character_, comment = NA_character_)), 
          type = "parameter"), compartments = new("compartments", 
          properties = new("compartment_properties", list = list(), 
              type = "compartment_property"), list = list(new("compartment", 
              name = "DEPOT", index = 1L), new("compartment", name = "CENTRAL", 
              index = 2L), new("compartment", name = "PERIPHERAL", 
              index = 3L), new("compartment", name = "OUTPUT", 
              index = 4L)), type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateDoseAdaptation.R:24:3'): Dose adaptations based on weight work well (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "KA", rhs = "THETA_KA*exp(ETA_KA)", 
              comment = NA_character_), new("equation", lhs = "CL", 
              rhs = "THETA_CL*exp(ETA_CL)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "V3", 
                  rhs = "THETA_V3*exp(ETA_V3)", comment = NA_character_), 
              new("equation", lhs = "Q", rhs = "THETA_Q*exp(ETA_Q)", 
                  comment = NA_character_), new("equation", lhs = "S2", 
                  rhs = "V2", comment = NA_character_)), type = "model_statement")), 
          new("ode_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("ode", lhs = "A_DEPOT", rhs = "-KA*A_DEPOT", 
                  comment = NA_character_), new("ode", lhs = "A_CENTRAL", 
                  rhs = "KA*A_DEPOT + Q*A_PERIPHERAL/V3 + (-CL/V2 - Q/V2)*A_CENTRAL", 
                  comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
                  rhs = "-Q*A_PERIPHERAL/V3 + Q*A_CENTRAL/V2", 
                  comment = NA_character_), new("ode", lhs = "A_OUTPUT", 
                  rhs = "CL*A_CENTRAL/V2", comment = NA_character_), 
                  new("equation", lhs = "F", rhs = "A_CENTRAL/S2", 
                    comment = NA_character_)), type = "model_statement")), 
          new("error_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("equation", lhs = "CP", rhs = "F", 
                  comment = NA_character_), new("equation", lhs = "OBS_CP", 
                  rhs = "CP*(EPS_PROP + 1)", comment = NA_character_), 
                  new("equation", lhs = "Y", rhs = "OBS_CP", comment = NA_character_)), 
              type = "model_statement"))), type = "code_record"), 
      parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
      0L)), list = list(new("theta", unit = NA_character_, name = "KA", 
          index = 1L, value = 1, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "CL", index = 2L, value = 5, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V2", index = 3L, value = 80, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V3", index = 4L, value = 20, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "Q", index = 5L, value = 4, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 1L, 
          type = "var", name = "KA", index = 1L, value = 0, fix = FALSE, 
          label = NA_character_, comment = NA_character_), new("omega", 
          same = NA, index2 = 2L, type = "var", name = "CL", index = 2L, 
          value = 0, fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 3L, type = "var", name = "V2", 
              index = 3L, value = 0, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 4L, type = "var", name = "V3", index = 4L, 
              value = 0, fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 5L, type = "var", name = "Q", 
              index = 5L, value = 0, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("sigma", index2 = 1L, 
              type = "var", name = "PROP", index = 1L, value = 0.025, 
              fix = FALSE, label = NA_character_, comment = NA_character_)), 
          type = "parameter"), compartments = new("compartments", 
          properties = new("compartment_properties", list = list(), 
              type = "compartment_property"), list = list(new("compartment", 
              name = "DEPOT", index = 1L), new("compartment", name = "CENTRAL", 
              index = 2L), new("compartment", name = "PERIPHERAL", 
              index = 3L), new("compartment", name = "OUTPUT", 
              index = 4L)), type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateDoseAdaptation.R:43:3'): Dose adaptations preserve specified infusion duration (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "CL", rhs = "THETA_CL*exp(ETA_CL)", 
              comment = NA_character_), new("equation", lhs = "V1", 
              rhs = "THETA_V1*exp(ETA_V1)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "Q", 
                  rhs = "THETA_Q*exp(ETA_Q)", comment = NA_character_), 
              new("equation", lhs = "S1", rhs = "V1", comment = NA_character_)), 
          type = "model_statement")), new("ode_record", comment = NA_character_, 
          statements = new("model_statements", list = list(new("ode", 
              lhs = "A_CENTRAL", rhs = "Q*A_PERIPHERAL/V2 + (-CL/V1 - Q/V1)*A_CENTRAL", 
              comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
              rhs = "-Q*A_PERIPHERAL/V2 + Q*A_CENTRAL/V1", comment = NA_character_), 
              new("ode", lhs = "A_OUTPUT", rhs = "CL*A_CENTRAL/V1", 
                  comment = NA_character_), new("equation", lhs = "F", 
                  rhs = "A_CENTRAL/S1", comment = NA_character_)), 
              type = "model_statement")), new("error_record", comment = NA_character_, 
          statements = new("model_statements", list = list(new("equation", 
              lhs = "CP", rhs = "F", comment = NA_character_), 
              new("equation", lhs = "OBS_CP", rhs = "CP*(EPS_PROP + 1)", 
                  comment = NA_character_), new("equation", lhs = "Y", 
                  rhs = "OBS_CP", comment = NA_character_)), type = "model_statement"))), 
      type = "code_record"), parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
  0L)), list = list(new("theta", unit = NA_character_, name = "CL", 
      index = 1L, value = 5, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("theta", unit = NA_character_, 
      name = "V1", index = 2L, value = 80, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("theta", unit = NA_character_, 
      name = "V2", index = 3L, value = 20, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("theta", unit = NA_character_, 
      name = "Q", index = 4L, value = 4, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("omega", same = NA, index2 = 1L, 
      type = "var", name = "CL", index = 1L, value = 0, fix = FALSE, 
      label = NA_character_, comment = NA_character_), new("omega", 
      same = NA, index2 = 2L, type = "var", name = "V1", index = 2L, 
      value = 0, fix = FALSE, label = NA_character_, comment = NA_character_), 
      new("omega", same = NA, index2 = 3L, type = "var", name = "V2", 
          index = 3L, value = 0, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 4L, 
          type = "var", name = "Q", index = 4L, value = 0, fix = FALSE, 
          label = NA_character_, comment = NA_character_), new("sigma", 
          index2 = 1L, type = "var", name = "PROP", index = 1L, 
          value = 0.025, fix = FALSE, label = NA_character_, comment = NA_character_)), 
      type = "parameter"), compartments = new("compartments", properties = new("compartment_properties", 
      list = list(), type = "compartment_property"), list = list(
      new("compartment", name = "CENTRAL", index = 1L), new("compartment", 
          name = "PERIPHERAL", index = 2L), new("compartment", 
          name = "OUTPUT", index = 3L)), type = "compartment"))), 
      .y = 1L, .f = function (...) 
      {
          NULL
          NULL
          ...furrr_out <- ...furrr_fn(...)
          ...furrr_out
      })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateDoseAdaptation.R:68:3'): Dose adaptations preserve specified infusion rate (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "CL", rhs = "THETA_CL*exp(ETA_CL)", 
              comment = NA_character_), new("equation", lhs = "V1", 
              rhs = "THETA_V1*exp(ETA_V1)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "Q", 
                  rhs = "THETA_Q*exp(ETA_Q)", comment = NA_character_), 
              new("equation", lhs = "S1", rhs = "V1", comment = NA_character_)), 
          type = "model_statement")), new("ode_record", comment = NA_character_, 
          statements = new("model_statements", list = list(new("ode", 
              lhs = "A_CENTRAL", rhs = "Q*A_PERIPHERAL/V2 + (-CL/V1 - Q/V1)*A_CENTRAL", 
              comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
              rhs = "-Q*A_PERIPHERAL/V2 + Q*A_CENTRAL/V1", comment = NA_character_), 
              new("ode", lhs = "A_OUTPUT", rhs = "CL*A_CENTRAL/V1", 
                  comment = NA_character_), new("equation", lhs = "F", 
                  rhs = "A_CENTRAL/S1", comment = NA_character_)), 
              type = "model_statement")), new("error_record", comment = NA_character_, 
          statements = new("model_statements", list = list(new("equation", 
              lhs = "CP", rhs = "F", comment = NA_character_), 
              new("equation", lhs = "OBS_CP", rhs = "CP*(EPS_PROP + 1)", 
                  comment = NA_character_), new("equation", lhs = "Y", 
                  rhs = "OBS_CP", comment = NA_character_)), type = "model_statement"))), 
      type = "code_record"), parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
  0L)), list = list(new("theta", unit = NA_character_, name = "CL", 
      index = 1L, value = 5, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("theta", unit = NA_character_, 
      name = "V1", index = 2L, value = 80, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("theta", unit = NA_character_, 
      name = "V2", index = 3L, value = 20, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("theta", unit = NA_character_, 
      name = "Q", index = 4L, value = 4, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("omega", same = NA, index2 = 1L, 
      type = "var", name = "CL", index = 1L, value = 0, fix = FALSE, 
      label = NA_character_, comment = NA_character_), new("omega", 
      same = NA, index2 = 2L, type = "var", name = "V1", index = 2L, 
      value = 0, fix = FALSE, label = NA_character_, comment = NA_character_), 
      new("omega", same = NA, index2 = 3L, type = "var", name = "V2", 
          index = 3L, value = 0, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 4L, 
          type = "var", name = "Q", index = 4L, value = 0, fix = FALSE, 
          label = NA_character_, comment = NA_character_), new("sigma", 
          index2 = 1L, type = "var", name = "PROP", index = 1L, 
          value = 0.025, fix = FALSE, label = NA_character_, comment = NA_character_)), 
      type = "parameter"), compartments = new("compartments", properties = new("compartment_properties", 
      list = list(), type = "compartment_property"), list = list(
      new("compartment", name = "CENTRAL", index = 1L), new("compartment", 
          name = "PERIPHERAL", index = 2L), new("compartment", 
          name = "OUTPUT", index = 3L)), type = "compartment"))), 
      .y = 1L, .f = function (...) 
      {
          NULL
          NULL
          ...furrr_out <- ...furrr_fn(...)
          ...furrr_out
      })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateDoseAdaptation.R:95:3'): Dose adaptations based on weight work well, check argument compartments works as expected (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "KA", rhs = "THETA_KA*exp(ETA_KA)", 
              comment = NA_character_), new("equation", lhs = "CL", 
              rhs = "THETA_CL*exp(ETA_CL)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "V3", 
                  rhs = "THETA_V3*exp(ETA_V3)", comment = NA_character_), 
              new("equation", lhs = "Q", rhs = "THETA_Q*exp(ETA_Q)", 
                  comment = NA_character_), new("equation", lhs = "S2", 
                  rhs = "V2", comment = NA_character_)), type = "model_statement")), 
          new("ode_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("ode", lhs = "A_DEPOT", rhs = "-KA*A_DEPOT", 
                  comment = NA_character_), new("ode", lhs = "A_CENTRAL", 
                  rhs = "KA*A_DEPOT + Q*A_PERIPHERAL/V3 + (-CL/V2 - Q/V2)*A_CENTRAL", 
                  comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
                  rhs = "-Q*A_PERIPHERAL/V3 + Q*A_CENTRAL/V2", 
                  comment = NA_character_), new("ode", lhs = "A_OUTPUT", 
                  rhs = "CL*A_CENTRAL/V2", comment = NA_character_), 
                  new("equation", lhs = "F", rhs = "A_CENTRAL/S2", 
                    comment = NA_character_), new("ode", lhs = "A_TEST", 
                    rhs = "-log(2)/100*A_TEST", comment = NA_character_)), 
              type = "model_statement")), new("error_record", comment = NA_character_, 
              statements = new("model_statements", list = list(
                  new("equation", lhs = "CP", rhs = "F", comment = NA_character_), 
                  new("equation", lhs = "OBS_CP", rhs = "CP*(EPS_PROP + 1)", 
                    comment = NA_character_), new("equation", lhs = "Y", 
                    rhs = "OBS_CP", comment = NA_character_)), 
                  type = "model_statement"))), type = "code_record"), 
      parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
      0L)), list = list(new("theta", unit = NA_character_, name = "KA", 
          index = 1L, value = 1, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "CL", index = 2L, value = 5, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V2", index = 3L, value = 80, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V3", index = 4L, value = 20, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "Q", index = 5L, value = 4, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 1L, 
          type = "var", name = "KA", index = 1L, value = 0, fix = FALSE, 
          label = NA_character_, comment = NA_character_), new("omega", 
          same = NA, index2 = 2L, type = "var", name = "CL", index = 2L, 
          value = 0, fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 3L, type = "var", name = "V2", 
              index = 3L, value = 0, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 4L, type = "var", name = "V3", index = 4L, 
              value = 0, fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 5L, type = "var", name = "Q", 
              index = 5L, value = 0, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("sigma", index2 = 1L, 
              type = "var", name = "PROP", index = 1L, value = 0.025, 
              fix = FALSE, label = NA_character_, comment = NA_character_)), 
          type = "parameter"), compartments = new("compartments", 
          properties = new("compartment_properties", list = list(), 
              type = "compartment_property"), list = list(new("compartment", 
              name = "DEPOT", index = 1L), new("compartment", name = "CENTRAL", 
              index = 2L), new("compartment", name = "PERIPHERAL", 
              index = 3L), new("compartment", name = "OUTPUT", 
              index = 4L), new("compartment", name = "TEST", index = 5L)), 
          type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateDosingInfo.R:24:3'): Dose adaptations can be checked in CAMPSIS output if dosing is TRUE (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "KA", rhs = "THETA_KA*exp(ETA_KA)", 
              comment = NA_character_), new("equation", lhs = "CL", 
              rhs = "THETA_CL*exp(ETA_CL)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "V3", 
                  rhs = "THETA_V3*exp(ETA_V3)", comment = NA_character_), 
              new("equation", lhs = "Q", rhs = "THETA_Q*exp(ETA_Q)", 
                  comment = NA_character_), new("equation", lhs = "S2", 
                  rhs = "V2", comment = NA_character_)), type = "model_statement")), 
          new("ode_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("ode", lhs = "A_DEPOT", rhs = "-KA*A_DEPOT", 
                  comment = NA_character_), new("ode", lhs = "A_CENTRAL", 
                  rhs = "KA*A_DEPOT + Q*A_PERIPHERAL/V3 + (-CL/V2 - Q/V2)*A_CENTRAL", 
                  comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
                  rhs = "-Q*A_PERIPHERAL/V3 + Q*A_CENTRAL/V2", 
                  comment = NA_character_), new("ode", lhs = "A_OUTPUT", 
                  rhs = "CL*A_CENTRAL/V2", comment = NA_character_), 
                  new("equation", lhs = "F", rhs = "A_CENTRAL/S2", 
                    comment = NA_character_)), type = "model_statement")), 
          new("error_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("equation", lhs = "CP", rhs = "F", 
                  comment = NA_character_), new("equation", lhs = "OBS_CP", 
                  rhs = "CP*(EPS_PROP + 1)", comment = NA_character_), 
                  new("equation", lhs = "Y", rhs = "OBS_CP", comment = NA_character_)), 
              type = "model_statement"))), type = "code_record"), 
      parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
      0L)), list = list(new("theta", unit = NA_character_, name = "KA", 
          index = 1L, value = 1, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "CL", index = 2L, value = 5, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V2", index = 3L, value = 80, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V3", index = 4L, value = 20, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "Q", index = 5L, value = 4, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 1L, 
          type = "var", name = "KA", index = 1L, value = 0.025, 
          fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 2L, type = "var", name = "CL", 
              index = 2L, value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 3L, type = "var", name = "V2", index = 3L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 4L, type = "var", name = "V3", index = 4L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 5L, type = "var", name = "Q", index = 5L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("sigma", index2 = 1L, 
              type = "var", name = "PROP", index = 1L, value = 0.025, 
              fix = FALSE, label = NA_character_, comment = NA_character_)), 
          type = "parameter"), compartments = new("compartments", 
          properties = new("compartment_properties", list = list(), 
              type = "compartment_property"), list = list(new("compartment", 
              name = "DEPOT", index = 1L), new("compartment", name = "CENTRAL", 
              index = 2L), new("compartment", name = "PERIPHERAL", 
              index = 3L), new("compartment", name = "OUTPUT", 
              index = 4L)), type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateInfusion.R:59:3'): Simulate infusion using rate in dataset (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "CL", rhs = "THETA_CL*exp(ETA_CL)", 
              comment = NA_character_), new("equation", lhs = "V1", 
              rhs = "THETA_V1*exp(ETA_V1)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "Q", 
                  rhs = "THETA_Q*exp(ETA_Q)", comment = NA_character_), 
              new("equation", lhs = "S1", rhs = "V1", comment = NA_character_)), 
          type = "model_statement")), new("ode_record", comment = NA_character_, 
          statements = new("model_statements", list = list(new("ode", 
              lhs = "A_CENTRAL", rhs = "Q*A_PERIPHERAL/V2 + (-CL/V1 - Q/V1)*A_CENTRAL", 
              comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
              rhs = "-Q*A_PERIPHERAL/V2 + Q*A_CENTRAL/V1", comment = NA_character_), 
              new("ode", lhs = "A_OUTPUT", rhs = "CL*A_CENTRAL/V1", 
                  comment = NA_character_), new("equation", lhs = "F", 
                  rhs = "A_CENTRAL/S1", comment = NA_character_)), 
              type = "model_statement")), new("error_record", comment = NA_character_, 
          statements = new("model_statements", list = list(new("equation", 
              lhs = "CP", rhs = "F", comment = NA_character_), 
              new("equation", lhs = "OBS_CP", rhs = "CP*(EPS_PROP + 1)", 
                  comment = NA_character_), new("equation", lhs = "Y", 
                  rhs = "OBS_CP", comment = NA_character_)), type = "model_statement"))), 
      type = "code_record"), parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
  0L)), list = list(new("theta", unit = NA_character_, name = "CL", 
      index = 1L, value = 5, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("theta", unit = NA_character_, 
      name = "V1", index = 2L, value = 80, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("theta", unit = NA_character_, 
      name = "V2", index = 3L, value = 20, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("theta", unit = NA_character_, 
      name = "Q", index = 4L, value = 4, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("omega", same = NA, index2 = 1L, 
      type = "var", name = "CL", index = 1L, value = 0.025, fix = FALSE, 
      label = NA_character_, comment = NA_character_), new("omega", 
      same = NA, index2 = 2L, type = "var", name = "V1", index = 2L, 
      value = 0.025, fix = FALSE, label = NA_character_, comment = NA_character_), 
      new("omega", same = NA, index2 = 3L, type = "var", name = "V2", 
          index = 3L, value = 0.025, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 4L, 
          type = "var", name = "Q", index = 4L, value = 0.025, 
          fix = FALSE, label = NA_character_, comment = NA_character_), 
      new("sigma", index2 = 1L, type = "var", name = "PROP", index = 1L, 
          value = 0.025, fix = FALSE, label = NA_character_, comment = NA_character_)), 
      type = "parameter"), compartments = new("compartments", properties = new("compartment_properties", 
      list = list(), type = "compartment_property"), list = list(
      new("compartment", name = "CENTRAL", index = 1L), new("compartment", 
          name = "PERIPHERAL", index = 2L), new("compartment", 
          name = "OUTPUT", index = 3L)), type = "compartment"))), 
      .y = 1L, .f = function (...) 
      {
          NULL
          NULL
          ...furrr_out <- ...furrr_fn(...)
          ...furrr_out
      })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateInfusion.R:96:3'): Simulate infusion using rate and lag time in dataset (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "CL", rhs = "THETA_CL*exp(ETA_CL)", 
              comment = NA_character_), new("equation", lhs = "V1", 
              rhs = "THETA_V1*exp(ETA_V1)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "Q", 
                  rhs = "THETA_Q*exp(ETA_Q)", comment = NA_character_), 
              new("equation", lhs = "S1", rhs = "V1", comment = NA_character_)), 
          type = "model_statement")), new("ode_record", comment = NA_character_, 
          statements = new("model_statements", list = list(new("ode", 
              lhs = "A_CENTRAL", rhs = "Q*A_PERIPHERAL/V2 + (-CL/V1 - Q/V1)*A_CENTRAL", 
              comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
              rhs = "-Q*A_PERIPHERAL/V2 + Q*A_CENTRAL/V1", comment = NA_character_), 
              new("ode", lhs = "A_OUTPUT", rhs = "CL*A_CENTRAL/V1", 
                  comment = NA_character_), new("equation", lhs = "F", 
                  rhs = "A_CENTRAL/S1", comment = NA_character_)), 
              type = "model_statement")), new("error_record", comment = NA_character_, 
          statements = new("model_statements", list = list(new("equation", 
              lhs = "CP", rhs = "F", comment = NA_character_), 
              new("equation", lhs = "OBS_CP", rhs = "CP*(EPS_PROP + 1)", 
                  comment = NA_character_), new("equation", lhs = "Y", 
                  rhs = "OBS_CP", comment = NA_character_)), type = "model_statement"))), 
      type = "code_record"), parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
  0L)), list = list(new("theta", unit = NA_character_, name = "CL", 
      index = 1L, value = 5, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("theta", unit = NA_character_, 
      name = "V1", index = 2L, value = 80, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("theta", unit = NA_character_, 
      name = "V2", index = 3L, value = 20, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("theta", unit = NA_character_, 
      name = "Q", index = 4L, value = 4, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("omega", same = NA, index2 = 1L, 
      type = "var", name = "CL", index = 1L, value = 0.025, fix = FALSE, 
      label = NA_character_, comment = NA_character_), new("omega", 
      same = NA, index2 = 2L, type = "var", name = "V1", index = 2L, 
      value = 0.025, fix = FALSE, label = NA_character_, comment = NA_character_), 
      new("omega", same = NA, index2 = 3L, type = "var", name = "V2", 
          index = 3L, value = 0.025, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 4L, 
          type = "var", name = "Q", index = 4L, value = 0.025, 
          fix = FALSE, label = NA_character_, comment = NA_character_), 
      new("sigma", index2 = 1L, type = "var", name = "PROP", index = 1L, 
          value = 0.025, fix = FALSE, label = NA_character_, comment = NA_character_)), 
      type = "parameter"), compartments = new("compartments", properties = new("compartment_properties", 
      list = list(), type = "compartment_property"), list = list(
      new("compartment", name = "CENTRAL", index = 1L), new("compartment", 
          name = "PERIPHERAL", index = 2L), new("compartment", 
          name = "OUTPUT", index = 3L)), type = "compartment"))), 
      .y = 1L, .f = function (...) 
      {
          NULL
          NULL
          ...furrr_out <- ...furrr_fn(...)
          ...furrr_out
      })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateInfusion.R:117:3'): Simulate infusion using rate and lag time (parameter distribution) in dataset (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "CL", rhs = "THETA_CL*exp(ETA_CL)", 
              comment = NA_character_), new("equation", lhs = "V1", 
              rhs = "THETA_V1*exp(ETA_V1)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "Q", 
                  rhs = "THETA_Q*exp(ETA_Q)", comment = NA_character_), 
              new("equation", lhs = "S1", rhs = "V1", comment = NA_character_)), 
          type = "model_statement")), new("ode_record", comment = NA_character_, 
          statements = new("model_statements", list = list(new("ode", 
              lhs = "A_CENTRAL", rhs = "Q*A_PERIPHERAL/V2 + (-CL/V1 - Q/V1)*A_CENTRAL", 
              comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
              rhs = "-Q*A_PERIPHERAL/V2 + Q*A_CENTRAL/V1", comment = NA_character_), 
              new("ode", lhs = "A_OUTPUT", rhs = "CL*A_CENTRAL/V1", 
                  comment = NA_character_), new("equation", lhs = "F", 
                  rhs = "A_CENTRAL/S1", comment = NA_character_)), 
              type = "model_statement")), new("error_record", comment = NA_character_, 
          statements = new("model_statements", list = list(new("equation", 
              lhs = "CP", rhs = "F", comment = NA_character_), 
              new("equation", lhs = "OBS_CP", rhs = "CP*(EPS_PROP + 1)", 
                  comment = NA_character_), new("equation", lhs = "Y", 
                  rhs = "OBS_CP", comment = NA_character_)), type = "model_statement"))), 
      type = "code_record"), parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
  0L)), list = list(new("theta", unit = NA_character_, name = "CL", 
      index = 1L, value = 5, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("theta", unit = NA_character_, 
      name = "V1", index = 2L, value = 80, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("theta", unit = NA_character_, 
      name = "V2", index = 3L, value = 20, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("theta", unit = NA_character_, 
      name = "Q", index = 4L, value = 4, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("omega", same = NA, index2 = 1L, 
      type = "var", name = "CL", index = 1L, value = 0.025, fix = FALSE, 
      label = NA_character_, comment = NA_character_), new("omega", 
      same = NA, index2 = 2L, type = "var", name = "V1", index = 2L, 
      value = 0.025, fix = FALSE, label = NA_character_, comment = NA_character_), 
      new("omega", same = NA, index2 = 3L, type = "var", name = "V2", 
          index = 3L, value = 0.025, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 4L, 
          type = "var", name = "Q", index = 4L, value = 0.025, 
          fix = FALSE, label = NA_character_, comment = NA_character_), 
      new("sigma", index2 = 1L, type = "var", name = "PROP", index = 1L, 
          value = 0.025, fix = FALSE, label = NA_character_, comment = NA_character_), 
      new("theta", unit = NA_character_, name = "ALAG1", index = 5L, 
          value = 2, fix = FALSE, label = NA_character_, comment = NA_character_), 
      new("omega", same = NA, index2 = 5L, type = "var", name = "ALAG1", 
          index = 5L, value = 0.04, fix = FALSE, label = NA_character_, 
          comment = NA_character_)), type = "parameter"), compartments = new("compartments", 
      properties = new("compartment_properties", list = list(), 
          type = "compartment_property"), list = list(new("compartment", 
          name = "CENTRAL", index = 1L), new("compartment", name = "PERIPHERAL", 
          index = 2L), new("compartment", name = "OUTPUT", index = 3L)), 
      type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateInit.R:21:3'): Simulate initial conditions, observations starting at 0 (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "CL", rhs = "THETA_CL*exp(ETA_CL)", 
              comment = NA_character_), new("equation", lhs = "V1", 
              rhs = "THETA_V1*exp(ETA_V1)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "Q", 
                  rhs = "THETA_Q*exp(ETA_Q)", comment = NA_character_), 
              new("equation", lhs = "S1", rhs = "V1", comment = NA_character_)), 
          type = "model_statement")), new("ode_record", comment = NA_character_, 
          statements = new("model_statements", list = list(new("ode", 
              lhs = "A_CENTRAL", rhs = "Q*A_PERIPHERAL/V2 + (-CL/V1 - Q/V1)*A_CENTRAL", 
              comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
              rhs = "-Q*A_PERIPHERAL/V2 + Q*A_CENTRAL/V1", comment = NA_character_), 
              new("ode", lhs = "A_OUTPUT", rhs = "CL*A_CENTRAL/V1", 
                  comment = NA_character_), new("equation", lhs = "F", 
                  rhs = "A_CENTRAL/S1", comment = NA_character_)), 
              type = "model_statement")), new("error_record", comment = NA_character_, 
          statements = new("model_statements", list = list(new("equation", 
              lhs = "CP", rhs = "F", comment = NA_character_), 
              new("equation", lhs = "OBS_CP", rhs = "CP*(EPS_PROP + 1)", 
                  comment = NA_character_), new("equation", lhs = "Y", 
                  rhs = "OBS_CP", comment = NA_character_)), type = "model_statement"))), 
      type = "code_record"), parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
  0L)), list = list(new("theta", unit = NA_character_, name = "CL", 
      index = 1L, value = 5, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("theta", unit = NA_character_, 
      name = "V1", index = 2L, value = 80, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("theta", unit = NA_character_, 
      name = "V2", index = 3L, value = 20, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("theta", unit = NA_character_, 
      name = "Q", index = 4L, value = 4, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("omega", same = NA, index2 = 1L, 
      type = "var", name = "CL", index = 1L, value = 0.025, fix = FALSE, 
      label = NA_character_, comment = NA_character_), new("omega", 
      same = NA, index2 = 2L, type = "var", name = "V1", index = 2L, 
      value = 0.025, fix = FALSE, label = NA_character_, comment = NA_character_), 
      new("omega", same = NA, index2 = 3L, type = "var", name = "V2", 
          index = 3L, value = 0.025, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 4L, 
          type = "var", name = "Q", index = 4L, value = 0.025, 
          fix = FALSE, label = NA_character_, comment = NA_character_), 
      new("sigma", index2 = 1L, type = "var", name = "PROP", index = 1L, 
          value = 0.025, fix = FALSE, label = NA_character_, comment = NA_character_)), 
      type = "parameter"), compartments = new("compartments", properties = new("compartment_properties", 
      list = list(new("compartment_initial_condition", compartment = 1L, 
          rhs = "1000", comment = NA_character_)), type = "compartment_property"), 
      list = list(new("compartment", name = "CENTRAL", index = 1L), 
          new("compartment", name = "PERIPHERAL", index = 2L), 
          new("compartment", name = "OUTPUT", index = 3L)), type = "compartment"))), 
      .y = 1L, .f = function (...) 
      {
          NULL
          NULL
          ...furrr_out <- ...furrr_fn(...)
          ...furrr_out
      })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateInit.R:37:3'): Simulate initial conditions, observations starting at 5 (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "CL", rhs = "THETA_CL*exp(ETA_CL)", 
              comment = NA_character_), new("equation", lhs = "V1", 
              rhs = "THETA_V1*exp(ETA_V1)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "Q", 
                  rhs = "THETA_Q*exp(ETA_Q)", comment = NA_character_), 
              new("equation", lhs = "S1", rhs = "V1", comment = NA_character_)), 
          type = "model_statement")), new("ode_record", comment = NA_character_, 
          statements = new("model_statements", list = list(new("ode", 
              lhs = "A_CENTRAL", rhs = "Q*A_PERIPHERAL/V2 + (-CL/V1 - Q/V1)*A_CENTRAL", 
              comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
              rhs = "-Q*A_PERIPHERAL/V2 + Q*A_CENTRAL/V1", comment = NA_character_), 
              new("ode", lhs = "A_OUTPUT", rhs = "CL*A_CENTRAL/V1", 
                  comment = NA_character_), new("equation", lhs = "F", 
                  rhs = "A_CENTRAL/S1", comment = NA_character_)), 
              type = "model_statement")), new("error_record", comment = NA_character_, 
          statements = new("model_statements", list = list(new("equation", 
              lhs = "CP", rhs = "F", comment = NA_character_), 
              new("equation", lhs = "OBS_CP", rhs = "CP*(EPS_PROP + 1)", 
                  comment = NA_character_), new("equation", lhs = "Y", 
                  rhs = "OBS_CP", comment = NA_character_)), type = "model_statement"))), 
      type = "code_record"), parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
  0L)), list = list(new("theta", unit = NA_character_, name = "CL", 
      index = 1L, value = 5, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("theta", unit = NA_character_, 
      name = "V1", index = 2L, value = 80, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("theta", unit = NA_character_, 
      name = "V2", index = 3L, value = 20, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("theta", unit = NA_character_, 
      name = "Q", index = 4L, value = 4, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("omega", same = NA, index2 = 1L, 
      type = "var", name = "CL", index = 1L, value = 0.025, fix = FALSE, 
      label = NA_character_, comment = NA_character_), new("omega", 
      same = NA, index2 = 2L, type = "var", name = "V1", index = 2L, 
      value = 0.025, fix = FALSE, label = NA_character_, comment = NA_character_), 
      new("omega", same = NA, index2 = 3L, type = "var", name = "V2", 
          index = 3L, value = 0.025, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 4L, 
          type = "var", name = "Q", index = 4L, value = 0.025, 
          fix = FALSE, label = NA_character_, comment = NA_character_), 
      new("sigma", index2 = 1L, type = "var", name = "PROP", index = 1L, 
          value = 0.025, fix = FALSE, label = NA_character_, comment = NA_character_)), 
      type = "parameter"), compartments = new("compartments", properties = new("compartment_properties", 
      list = list(new("compartment_initial_condition", compartment = 1L, 
          rhs = "1000", comment = NA_character_)), type = "compartment_property"), 
      list = list(new("compartment", name = "CENTRAL", index = 1L), 
          new("compartment", name = "PERIPHERAL", index = 2L), 
          new("compartment", name = "OUTPUT", index = 3L)), type = "compartment"))), 
      .y = 1L, .f = function (...) 
      {
          NULL
          NULL
          ...furrr_out <- ...furrr_fn(...)
          ...furrr_out
      })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateLagTime.R:23:3'): Simulate a bolus with fixed lag time in dataset (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "KA", rhs = "THETA_KA*exp(ETA_KA)", 
              comment = NA_character_), new("equation", lhs = "CL", 
              rhs = "THETA_CL*exp(ETA_CL)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "V3", 
                  rhs = "THETA_V3*exp(ETA_V3)", comment = NA_character_), 
              new("equation", lhs = "Q", rhs = "THETA_Q*exp(ETA_Q)", 
                  comment = NA_character_), new("equation", lhs = "S2", 
                  rhs = "V2", comment = NA_character_)), type = "model_statement")), 
          new("ode_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("ode", lhs = "A_DEPOT", rhs = "-KA*A_DEPOT", 
                  comment = NA_character_), new("ode", lhs = "A_CENTRAL", 
                  rhs = "KA*A_DEPOT + Q*A_PERIPHERAL/V3 + (-CL/V2 - Q/V2)*A_CENTRAL", 
                  comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
                  rhs = "-Q*A_PERIPHERAL/V3 + Q*A_CENTRAL/V2", 
                  comment = NA_character_), new("ode", lhs = "A_OUTPUT", 
                  rhs = "CL*A_CENTRAL/V2", comment = NA_character_), 
                  new("equation", lhs = "F", rhs = "A_CENTRAL/S2", 
                    comment = NA_character_)), type = "model_statement")), 
          new("error_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("equation", lhs = "CP", rhs = "F", 
                  comment = NA_character_), new("equation", lhs = "OBS_CP", 
                  rhs = "CP*(EPS_PROP + 1)", comment = NA_character_), 
                  new("equation", lhs = "Y", rhs = "OBS_CP", comment = NA_character_)), 
              type = "model_statement"))), type = "code_record"), 
      parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
      0L)), list = list(new("theta", unit = NA_character_, name = "KA", 
          index = 1L, value = 1, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "CL", index = 2L, value = 5, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V2", index = 3L, value = 80, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V3", index = 4L, value = 20, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "Q", index = 5L, value = 4, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 1L, 
          type = "var", name = "KA", index = 1L, value = 0.025, 
          fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 2L, type = "var", name = "CL", 
              index = 2L, value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 3L, type = "var", name = "V2", index = 3L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 4L, type = "var", name = "V3", index = 4L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 5L, type = "var", name = "Q", index = 5L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("sigma", index2 = 1L, 
              type = "var", name = "PROP", index = 1L, value = 0.025, 
              fix = FALSE, label = NA_character_, comment = NA_character_)), 
          type = "parameter"), compartments = new("compartments", 
          properties = new("compartment_properties", list = list(), 
              type = "compartment_property"), list = list(new("compartment", 
              name = "DEPOT", index = 1L), new("compartment", name = "CENTRAL", 
              index = 2L), new("compartment", name = "PERIPHERAL", 
              index = 3L), new("compartment", name = "OUTPUT", 
              index = 4L)), type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateLagTime.R:44:3'): Simulate a bolus with fixed lag time in model (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "KA", rhs = "THETA_KA*exp(ETA_KA)", 
              comment = NA_character_), new("equation", lhs = "CL", 
              rhs = "THETA_CL*exp(ETA_CL)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "V3", 
                  rhs = "THETA_V3*exp(ETA_V3)", comment = NA_character_), 
              new("equation", lhs = "Q", rhs = "THETA_Q*exp(ETA_Q)", 
                  comment = NA_character_), new("equation", lhs = "S2", 
                  rhs = "V2", comment = NA_character_)), type = "model_statement")), 
          new("ode_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("ode", lhs = "A_DEPOT", rhs = "-KA*A_DEPOT", 
                  comment = NA_character_), new("ode", lhs = "A_CENTRAL", 
                  rhs = "KA*A_DEPOT + Q*A_PERIPHERAL/V3 + (-CL/V2 - Q/V2)*A_CENTRAL", 
                  comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
                  rhs = "-Q*A_PERIPHERAL/V3 + Q*A_CENTRAL/V2", 
                  comment = NA_character_), new("ode", lhs = "A_OUTPUT", 
                  rhs = "CL*A_CENTRAL/V2", comment = NA_character_), 
                  new("equation", lhs = "F", rhs = "A_CENTRAL/S2", 
                    comment = NA_character_)), type = "model_statement")), 
          new("error_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("equation", lhs = "CP", rhs = "F", 
                  comment = NA_character_), new("equation", lhs = "OBS_CP", 
                  rhs = "CP*(EPS_PROP + 1)", comment = NA_character_), 
                  new("equation", lhs = "Y", rhs = "OBS_CP", comment = NA_character_)), 
              type = "model_statement"))), type = "code_record"), 
      parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
      0L)), list = list(new("theta", unit = NA_character_, name = "KA", 
          index = 1L, value = 1, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "CL", index = 2L, value = 5, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V2", index = 3L, value = 80, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V3", index = 4L, value = 20, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "Q", index = 5L, value = 4, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 1L, 
          type = "var", name = "KA", index = 1L, value = 0.025, 
          fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 2L, type = "var", name = "CL", 
              index = 2L, value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 3L, type = "var", name = "V2", index = 3L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 4L, type = "var", name = "V3", index = 4L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 5L, type = "var", name = "Q", index = 5L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("sigma", index2 = 1L, 
              type = "var", name = "PROP", index = 1L, value = 0.025, 
              fix = FALSE, label = NA_character_, comment = NA_character_)), 
          type = "parameter"), compartments = new("compartments", 
          properties = new("compartment_properties", list = list(
              new("compartment_lag_time", compartment = 1L, rhs = "2", 
                  comment = NA_character_)), type = "compartment_property"), 
          list = list(new("compartment", name = "DEPOT", index = 1L), 
              new("compartment", name = "CENTRAL", index = 2L), 
              new("compartment", name = "PERIPHERAL", index = 3L), 
              new("compartment", name = "OUTPUT", index = 4L)), 
          type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateMinimalistModel.R:20:3'): Simulate a minimalist model, no parameters/IIV/error model (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "THALF", rhs = "12", 
              comment = NA_character_), new("equation", lhs = "K", 
              rhs = "log(2)/THALF", comment = NA_character_)), 
          type = "model_statement")), new("ode_record", comment = NA_character_, 
          statements = new("model_statements", list = list(new("ode", 
              lhs = "A_CENTRAL", rhs = "-K*A_CENTRAL", comment = NA_character_)), 
              type = "model_statement"))), type = "code_record"), 
      parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
      0L)), list = list(), type = "parameter"), compartments = new("compartments", 
          properties = new("compartment_properties", list = list(
              new("compartment_initial_condition", compartment = 1L, 
                  rhs = "1000", comment = NA_character_)), type = "compartment_property"), 
          list = list(new("compartment", name = "CENTRAL", index = 1L)), 
          type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateNoOmega.R:30:3'): Simulate a 1-cpt model without OMEGA's/SIGMA's - Github issue #8 (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "KA", rhs = "THETA_KA*exp(ETA_KA)", 
              comment = NA_character_), new("equation", lhs = "CL", 
              rhs = "THETA_CL*exp(ETA_CL)", comment = NA_character_), 
              new("equation", lhs = "V", rhs = "THETA_V*exp(ETA_V)", 
                  comment = NA_character_), new("equation", lhs = "S2", 
                  rhs = "V", comment = NA_character_)), type = "model_statement")), 
          new("ode_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("ode", lhs = "A_DEPOT", rhs = "-KA*A_DEPOT", 
                  comment = NA_character_), new("ode", lhs = "A_CENTRAL", 
                  rhs = "-CL*A_CENTRAL/V + KA*A_DEPOT", comment = NA_character_), 
                  new("ode", lhs = "A_OUTPUT", rhs = "CL*A_CENTRAL/V", 
                    comment = NA_character_), new("equation", lhs = "F", 
                    rhs = "A_CENTRAL/S2", comment = NA_character_)), 
              type = "model_statement")), new("error_record", comment = NA_character_, 
              statements = new("model_statements", list = list(
                  new("equation", lhs = "CP", rhs = "F", comment = NA_character_), 
                  new("equation", lhs = "OBS_CP", rhs = "CP*(EPS_PROP + 1)", 
                    comment = NA_character_), new("equation", lhs = "Y", 
                    rhs = "OBS_CP", comment = NA_character_)), 
                  type = "model_statement"))), type = "code_record"), 
      parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
      0L)), list = list(new("theta", unit = NA_character_, name = "KA", 
          index = 1L, value = 1, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "CL", index = 2L, value = 5, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V", index = 3L, value = 80, fix = FALSE, label = NA_character_, 
          comment = NA_character_)), type = "parameter"), compartments = new("compartments", 
          properties = new("compartment_properties", list = list(), 
              type = "compartment_property"), list = list(new("compartment", 
              name = "DEPOT", index = 1L), new("compartment", name = "CENTRAL", 
              index = 2L), new("compartment", name = "OUTPUT", 
              index = 3L)), type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateOutvars.R:23:3'): NULL outvars (rxode2/mrgsolve) ────────
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "KA", rhs = "THETA_KA*exp(ETA_KA)", 
              comment = NA_character_), new("equation", lhs = "CL", 
              rhs = "THETA_CL*exp(ETA_CL)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "V3", 
                  rhs = "THETA_V3*exp(ETA_V3)", comment = NA_character_), 
              new("equation", lhs = "Q", rhs = "THETA_Q*exp(ETA_Q)", 
                  comment = NA_character_), new("equation", lhs = "S2", 
                  rhs = "V2", comment = NA_character_)), type = "model_statement")), 
          new("ode_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("ode", lhs = "A_DEPOT", rhs = "-KA*A_DEPOT", 
                  comment = NA_character_), new("ode", lhs = "A_CENTRAL", 
                  rhs = "KA*A_DEPOT + Q*A_PERIPHERAL/V3 + (-CL/V2 - Q/V2)*A_CENTRAL", 
                  comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
                  rhs = "-Q*A_PERIPHERAL/V3 + Q*A_CENTRAL/V2", 
                  comment = NA_character_), new("ode", lhs = "A_OUTPUT", 
                  rhs = "CL*A_CENTRAL/V2", comment = NA_character_), 
                  new("equation", lhs = "F", rhs = "A_CENTRAL/S2", 
                    comment = NA_character_)), type = "model_statement")), 
          new("error_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("equation", lhs = "CP", rhs = "F", 
                  comment = NA_character_), new("equation", lhs = "OBS_CP", 
                  rhs = "CP*(EPS_PROP + 1)", comment = NA_character_), 
                  new("equation", lhs = "Y", rhs = "OBS_CP", comment = NA_character_)), 
              type = "model_statement"))), type = "code_record"), 
      parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
      0L)), list = list(new("theta", unit = NA_character_, name = "KA", 
          index = 1L, value = 1, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "CL", index = 2L, value = 5, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V2", index = 3L, value = 80, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V3", index = 4L, value = 20, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "Q", index = 5L, value = 4, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 1L, 
          type = "var", name = "KA", index = 1L, value = 0.025, 
          fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 2L, type = "var", name = "CL", 
              index = 2L, value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 3L, type = "var", name = "V2", index = 3L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 4L, type = "var", name = "V3", index = 4L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 5L, type = "var", name = "Q", index = 5L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("sigma", index2 = 1L, 
              type = "var", name = "PROP", index = 1L, value = 0.025, 
              fix = FALSE, label = NA_character_, comment = NA_character_)), 
          type = "parameter"), compartments = new("compartments", 
          properties = new("compartment_properties", list = list(), 
              type = "compartment_property"), list = list(new("compartment", 
              name = "DEPOT", index = 1L), new("compartment", name = "CENTRAL", 
              index = 2L), new("compartment", name = "PERIPHERAL", 
              index = 3L), new("compartment", name = "OUTPUT", 
              index = 4L)), type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateOutvars.R:42:3'): Not NULL outvars (rxode2/mrgsolve) ────
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "KA", rhs = "THETA_KA*exp(ETA_KA)", 
              comment = NA_character_), new("equation", lhs = "CL", 
              rhs = "THETA_CL*exp(ETA_CL)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "V3", 
                  rhs = "THETA_V3*exp(ETA_V3)", comment = NA_character_), 
              new("equation", lhs = "Q", rhs = "THETA_Q*exp(ETA_Q)", 
                  comment = NA_character_), new("equation", lhs = "S2", 
                  rhs = "V2", comment = NA_character_)), type = "model_statement")), 
          new("ode_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("ode", lhs = "A_DEPOT", rhs = "-KA*A_DEPOT", 
                  comment = NA_character_), new("ode", lhs = "A_CENTRAL", 
                  rhs = "KA*A_DEPOT + Q*A_PERIPHERAL/V3 + (-CL/V2 - Q/V2)*A_CENTRAL", 
                  comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
                  rhs = "-Q*A_PERIPHERAL/V3 + Q*A_CENTRAL/V2", 
                  comment = NA_character_), new("ode", lhs = "A_OUTPUT", 
                  rhs = "CL*A_CENTRAL/V2", comment = NA_character_), 
                  new("equation", lhs = "F", rhs = "A_CENTRAL/S2", 
                    comment = NA_character_)), type = "model_statement")), 
          new("error_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("equation", lhs = "CP", rhs = "F", 
                  comment = NA_character_), new("equation", lhs = "OBS_CP", 
                  rhs = "CP*(EPS_PROP + 1)", comment = NA_character_), 
                  new("equation", lhs = "Y", rhs = "OBS_CP", comment = NA_character_)), 
              type = "model_statement"))), type = "code_record"), 
      parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
      0L)), list = list(new("theta", unit = NA_character_, name = "KA", 
          index = 1L, value = 1, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "CL", index = 2L, value = 5, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V2", index = 3L, value = 80, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V3", index = 4L, value = 20, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "Q", index = 5L, value = 4, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 1L, 
          type = "var", name = "KA", index = 1L, value = 0.025, 
          fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 2L, type = "var", name = "CL", 
              index = 2L, value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 3L, type = "var", name = "V2", index = 3L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 4L, type = "var", name = "V3", index = 4L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 5L, type = "var", name = "Q", index = 5L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("sigma", index2 = 1L, 
              type = "var", name = "PROP", index = 1L, value = 0.025, 
              fix = FALSE, label = NA_character_, comment = NA_character_)), 
          type = "parameter"), compartments = new("compartments", 
          properties = new("compartment_properties", list = list(), 
              type = "compartment_property"), list = list(new("compartment", 
              name = "DEPOT", index = 1L), new("compartment", name = "CENTRAL", 
              index = 2L), new("compartment", name = "PERIPHERAL", 
              index = 3L), new("compartment", name = "OUTPUT", 
              index = 4L)), type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateOutvars.R:57:3'): Not NULL outvars + DROP_OTHERS (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "KA", rhs = "THETA_KA*exp(ETA_KA)", 
              comment = NA_character_), new("equation", lhs = "CL", 
              rhs = "THETA_CL*exp(ETA_CL)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "V3", 
                  rhs = "THETA_V3*exp(ETA_V3)", comment = NA_character_), 
              new("equation", lhs = "Q", rhs = "THETA_Q*exp(ETA_Q)", 
                  comment = NA_character_), new("equation", lhs = "S2", 
                  rhs = "V2", comment = NA_character_)), type = "model_statement")), 
          new("ode_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("ode", lhs = "A_DEPOT", rhs = "-KA*A_DEPOT", 
                  comment = NA_character_), new("ode", lhs = "A_CENTRAL", 
                  rhs = "KA*A_DEPOT + Q*A_PERIPHERAL/V3 + (-CL/V2 - Q/V2)*A_CENTRAL", 
                  comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
                  rhs = "-Q*A_PERIPHERAL/V3 + Q*A_CENTRAL/V2", 
                  comment = NA_character_), new("ode", lhs = "A_OUTPUT", 
                  rhs = "CL*A_CENTRAL/V2", comment = NA_character_), 
                  new("equation", lhs = "F", rhs = "A_CENTRAL/S2", 
                    comment = NA_character_)), type = "model_statement")), 
          new("error_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("equation", lhs = "CP", rhs = "F", 
                  comment = NA_character_), new("equation", lhs = "OBS_CP", 
                  rhs = "CP*(EPS_PROP + 1)", comment = NA_character_), 
                  new("equation", lhs = "Y", rhs = "OBS_CP", comment = NA_character_)), 
              type = "model_statement"))), type = "code_record"), 
      parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
      0L)), list = list(new("theta", unit = NA_character_, name = "KA", 
          index = 1L, value = 1, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "CL", index = 2L, value = 5, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V2", index = 3L, value = 80, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V3", index = 4L, value = 20, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "Q", index = 5L, value = 4, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 1L, 
          type = "var", name = "KA", index = 1L, value = 0.025, 
          fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 2L, type = "var", name = "CL", 
              index = 2L, value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 3L, type = "var", name = "V2", index = 3L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 4L, type = "var", name = "V3", index = 4L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 5L, type = "var", name = "Q", index = 5L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("sigma", index2 = 1L, 
              type = "var", name = "PROP", index = 1L, value = 0.025, 
              fix = FALSE, label = NA_character_, comment = NA_character_)), 
          type = "parameter"), compartments = new("compartments", 
          properties = new("compartment_properties", list = list(), 
              type = "compartment_property"), list = list(new("compartment", 
              name = "DEPOT", index = 1L), new("compartment", name = "CENTRAL", 
              index = 2L), new("compartment", name = "PERIPHERAL", 
              index = 3L), new("compartment", name = "OUTPUT", 
              index = 4L)), type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateOutvars.R:72:3'): Not NULL outvars from ERROR block + DROP_OTHERS (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "KA", rhs = "THETA_KA*exp(ETA_KA)", 
              comment = NA_character_), new("equation", lhs = "CL", 
              rhs = "THETA_CL*exp(ETA_CL)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "V3", 
                  rhs = "THETA_V3*exp(ETA_V3)", comment = NA_character_), 
              new("equation", lhs = "Q", rhs = "THETA_Q*exp(ETA_Q)", 
                  comment = NA_character_), new("equation", lhs = "S2", 
                  rhs = "V2", comment = NA_character_)), type = "model_statement")), 
          new("ode_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("ode", lhs = "A_DEPOT", rhs = "-KA*A_DEPOT", 
                  comment = NA_character_), new("ode", lhs = "A_CENTRAL", 
                  rhs = "KA*A_DEPOT + Q*A_PERIPHERAL/V3 + (-CL/V2 - Q/V2)*A_CENTRAL", 
                  comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
                  rhs = "-Q*A_PERIPHERAL/V3 + Q*A_CENTRAL/V2", 
                  comment = NA_character_), new("ode", lhs = "A_OUTPUT", 
                  rhs = "CL*A_CENTRAL/V2", comment = NA_character_), 
                  new("equation", lhs = "F", rhs = "A_CENTRAL/S2", 
                    comment = NA_character_)), type = "model_statement")), 
          new("error_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("equation", lhs = "CP", rhs = "F", 
                  comment = NA_character_), new("equation", lhs = "OBS_CP", 
                  rhs = "CP*(EPS_PROP + 1)", comment = NA_character_), 
                  new("equation", lhs = "Y", rhs = "OBS_CP", comment = NA_character_)), 
              type = "model_statement"))), type = "code_record"), 
      parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
      0L)), list = list(new("theta", unit = NA_character_, name = "KA", 
          index = 1L, value = 1, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "CL", index = 2L, value = 5, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V2", index = 3L, value = 80, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V3", index = 4L, value = 20, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "Q", index = 5L, value = 4, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 1L, 
          type = "var", name = "KA", index = 1L, value = 0.025, 
          fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 2L, type = "var", name = "CL", 
              index = 2L, value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 3L, type = "var", name = "V2", index = 3L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 4L, type = "var", name = "V3", index = 4L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 5L, type = "var", name = "Q", index = 5L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("sigma", index2 = 1L, 
              type = "var", name = "PROP", index = 1L, value = 0.025, 
              fix = FALSE, label = NA_character_, comment = NA_character_)), 
          type = "parameter"), compartments = new("compartments", 
          properties = new("compartment_properties", list = list(), 
              type = "compartment_property"), list = list(new("compartment", 
              name = "DEPOT", index = 1L), new("compartment", name = "CENTRAL", 
              index = 2L), new("compartment", name = "PERIPHERAL", 
              index = 3L), new("compartment", name = "OUTPUT", 
              index = 4L)), type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateOutvars.R:89:3'): Covariates in outvars can be output well (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "KA", rhs = "THETA_KA*exp(ETA_KA)", 
              comment = NA_character_), new("equation", lhs = "CL", 
              rhs = "THETA_CL*exp(ETA_CL)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "V3", 
                  rhs = "THETA_V3*exp(ETA_V3)", comment = NA_character_), 
              new("equation", lhs = "Q", rhs = "THETA_Q*exp(ETA_Q)", 
                  comment = NA_character_), new("equation", lhs = "S2", 
                  rhs = "V2", comment = NA_character_)), type = "model_statement")), 
          new("ode_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("ode", lhs = "A_DEPOT", rhs = "-KA*A_DEPOT", 
                  comment = NA_character_), new("ode", lhs = "A_CENTRAL", 
                  rhs = "KA*A_DEPOT + Q*A_PERIPHERAL/V3 + (-CL/V2 - Q/V2)*A_CENTRAL", 
                  comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
                  rhs = "-Q*A_PERIPHERAL/V3 + Q*A_CENTRAL/V2", 
                  comment = NA_character_), new("ode", lhs = "A_OUTPUT", 
                  rhs = "CL*A_CENTRAL/V2", comment = NA_character_), 
                  new("equation", lhs = "F", rhs = "A_CENTRAL/S2", 
                    comment = NA_character_)), type = "model_statement")), 
          new("error_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("equation", lhs = "CP", rhs = "F", 
                  comment = NA_character_), new("equation", lhs = "OBS_CP", 
                  rhs = "CP*(EPS_PROP + 1)", comment = NA_character_), 
                  new("equation", lhs = "Y", rhs = "OBS_CP", comment = NA_character_)), 
              type = "model_statement"))), type = "code_record"), 
      parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
      0L)), list = list(new("theta", unit = NA_character_, name = "KA", 
          index = 1L, value = 1, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "CL", index = 2L, value = 5, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V2", index = 3L, value = 80, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V3", index = 4L, value = 20, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "Q", index = 5L, value = 4, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 1L, 
          type = "var", name = "KA", index = 1L, value = 0.025, 
          fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 2L, type = "var", name = "CL", 
              index = 2L, value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 3L, type = "var", name = "V2", index = 3L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 4L, type = "var", name = "V3", index = 4L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 5L, type = "var", name = "Q", index = 5L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("sigma", index2 = 1L, 
              type = "var", name = "PROP", index = 1L, value = 0.025, 
              fix = FALSE, label = NA_character_, comment = NA_character_)), 
          type = "parameter"), compartments = new("compartments", 
          properties = new("compartment_properties", list = list(), 
              type = "compartment_property"), list = list(new("compartment", 
              name = "DEPOT", index = 1L), new("compartment", name = "CENTRAL", 
              index = 2L), new("compartment", name = "PERIPHERAL", 
              index = 3L), new("compartment", name = "OUTPUT", 
              index = 4L)), type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateOutvars.R:105:3'): ETAs in outvars can be output well (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "KA", rhs = "THETA_KA*exp(ETA_KA)", 
              comment = NA_character_), new("equation", lhs = "CL", 
              rhs = "THETA_CL*exp(ETA_CL)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "V3", 
                  rhs = "THETA_V3*exp(ETA_V3)", comment = NA_character_), 
              new("equation", lhs = "Q", rhs = "THETA_Q*exp(ETA_Q)", 
                  comment = NA_character_), new("equation", lhs = "S2", 
                  rhs = "V2", comment = NA_character_)), type = "model_statement")), 
          new("ode_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("ode", lhs = "A_DEPOT", rhs = "-KA*A_DEPOT", 
                  comment = NA_character_), new("ode", lhs = "A_CENTRAL", 
                  rhs = "KA*A_DEPOT + Q*A_PERIPHERAL/V3 + (-CL/V2 - Q/V2)*A_CENTRAL", 
                  comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
                  rhs = "-Q*A_PERIPHERAL/V3 + Q*A_CENTRAL/V2", 
                  comment = NA_character_), new("ode", lhs = "A_OUTPUT", 
                  rhs = "CL*A_CENTRAL/V2", comment = NA_character_), 
                  new("equation", lhs = "F", rhs = "A_CENTRAL/S2", 
                    comment = NA_character_)), type = "model_statement")), 
          new("error_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("equation", lhs = "CP", rhs = "F", 
                  comment = NA_character_), new("equation", lhs = "OBS_CP", 
                  rhs = "CP*(EPS_PROP + 1)", comment = NA_character_), 
                  new("equation", lhs = "Y", rhs = "OBS_CP", comment = NA_character_)), 
              type = "model_statement"))), type = "code_record"), 
      parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
      0L)), list = list(new("theta", unit = NA_character_, name = "KA", 
          index = 1L, value = 1, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "CL", index = 2L, value = 5, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V2", index = 3L, value = 80, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V3", index = 4L, value = 20, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "Q", index = 5L, value = 4, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 1L, 
          type = "var", name = "KA", index = 1L, value = 0.025, 
          fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 2L, type = "var", name = "CL", 
              index = 2L, value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 3L, type = "var", name = "V2", index = 3L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 4L, type = "var", name = "V3", index = 4L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 5L, type = "var", name = "Q", index = 5L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("sigma", index2 = 1L, 
              type = "var", name = "PROP", index = 1L, value = 0.025, 
              fix = FALSE, label = NA_character_, comment = NA_character_)), 
          type = "parameter"), compartments = new("compartments", 
          properties = new("compartment_properties", list = list(), 
              type = "compartment_property"), list = list(new("compartment", 
              name = "DEPOT", index = 1L), new("compartment", name = "CENTRAL", 
              index = 2L), new("compartment", name = "PERIPHERAL", 
              index = 3L), new("compartment", name = "OUTPUT", 
              index = 4L)), type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateRUV.R:25:3'): Test generated RUV is correct (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "KA", rhs = "THETA_KA*exp(ETA_KA)", 
              comment = NA_character_), new("equation", lhs = "CL", 
              rhs = "THETA_CL*exp(ETA_CL)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "V3", 
                  rhs = "THETA_V3*exp(ETA_V3)", comment = NA_character_), 
              new("equation", lhs = "Q", rhs = "THETA_Q*exp(ETA_Q)", 
                  comment = NA_character_), new("equation", lhs = "S2", 
                  rhs = "V2", comment = NA_character_)), type = "model_statement")), 
          new("ode_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("ode", lhs = "A_DEPOT", rhs = "-KA*A_DEPOT", 
                  comment = NA_character_), new("ode", lhs = "A_CENTRAL", 
                  rhs = "KA*A_DEPOT + Q*A_PERIPHERAL/V3 + (-CL/V2 - Q/V2)*A_CENTRAL", 
                  comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
                  rhs = "-Q*A_PERIPHERAL/V3 + Q*A_CENTRAL/V2", 
                  comment = NA_character_), new("ode", lhs = "A_OUTPUT", 
                  rhs = "CL*A_CENTRAL/V2", comment = NA_character_), 
                  new("equation", lhs = "F", rhs = "A_CENTRAL/S2", 
                    comment = NA_character_)), type = "model_statement")), 
          new("error_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("equation", lhs = "CP", rhs = "F", 
                  comment = NA_character_), new("equation", lhs = "OBS_CP", 
                  rhs = "CP*(EPS_PROP + 1)", comment = NA_character_), 
                  new("equation", lhs = "Y", rhs = "OBS_CP", comment = NA_character_)), 
              type = "model_statement"))), type = "code_record"), 
      parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
      0L)), list = list(new("theta", unit = NA_character_, name = "KA", 
          index = 1L, value = 1, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "CL", index = 2L, value = 5, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V2", index = 3L, value = 80, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V3", index = 4L, value = 20, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "Q", index = 5L, value = 4, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 1L, 
          type = "var", name = "KA", index = 1L, value = 0.025, 
          fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 2L, type = "var", name = "CL", 
              index = 2L, value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 3L, type = "var", name = "V2", index = 3L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 4L, type = "var", name = "V3", index = 4L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 5L, type = "var", name = "Q", index = 5L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("sigma", index2 = 1L, 
              type = "var", name = "PROP", index = 1L, value = 0.025, 
              fix = FALSE, label = NA_character_, comment = NA_character_)), 
          type = "parameter"), compartments = new("compartments", 
          properties = new("compartment_properties", list = list(), 
              type = "compartment_property"), list = list(new("compartment", 
              name = "DEPOT", index = 1L), new("compartment", name = "CENTRAL", 
              index = 2L), new("compartment", name = "PERIPHERAL", 
              index = 3L), new("compartment", name = "OUTPUT", 
              index = 4L)), type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateScenarios.R:44:3'): Simulate scenarios - make few changes on dataset (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "KA", rhs = "THETA_KA*exp(ETA_KA)", 
              comment = NA_character_), new("equation", lhs = "CL", 
              rhs = "THETA_CL*exp(ETA_CL)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "V3", 
                  rhs = "THETA_V3*exp(ETA_V3)", comment = NA_character_), 
              new("equation", lhs = "Q", rhs = "THETA_Q*exp(ETA_Q)", 
                  comment = NA_character_), new("equation", lhs = "S2", 
                  rhs = "V2", comment = NA_character_)), type = "model_statement")), 
          new("ode_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("ode", lhs = "A_DEPOT", rhs = "-KA*A_DEPOT", 
                  comment = NA_character_), new("ode", lhs = "A_CENTRAL", 
                  rhs = "KA*A_DEPOT + Q*A_PERIPHERAL/V3 + (-CL/V2 - Q/V2)*A_CENTRAL", 
                  comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
                  rhs = "-Q*A_PERIPHERAL/V3 + Q*A_CENTRAL/V2", 
                  comment = NA_character_), new("ode", lhs = "A_OUTPUT", 
                  rhs = "CL*A_CENTRAL/V2", comment = NA_character_), 
                  new("equation", lhs = "F", rhs = "A_CENTRAL/S2", 
                    comment = NA_character_)), type = "model_statement")), 
          new("error_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("equation", lhs = "CP", rhs = "F", 
                  comment = NA_character_), new("equation", lhs = "OBS_CP", 
                  rhs = "CP*(EPS_PROP + 1)", comment = NA_character_), 
                  new("equation", lhs = "Y", rhs = "OBS_CP", comment = NA_character_)), 
              type = "model_statement"))), type = "code_record"), 
      parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
      0L)), list = list(new("theta", unit = NA_character_, name = "KA", 
          index = 1L, value = 1, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "CL", index = 2L, value = 5, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V2", index = 3L, value = 80, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V3", index = 4L, value = 20, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "Q", index = 5L, value = 4, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 1L, 
          type = "var", name = "KA", index = 1L, value = 0.025, 
          fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 2L, type = "var", name = "CL", 
              index = 2L, value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 3L, type = "var", name = "V2", index = 3L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 4L, type = "var", name = "V3", index = 4L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 5L, type = "var", name = "Q", index = 5L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("sigma", index2 = 1L, 
              type = "var", name = "PROP", index = 1L, value = 0.025, 
              fix = FALSE, label = NA_character_, comment = NA_character_)), 
          type = "parameter"), compartments = new("compartments", 
          properties = new("compartment_properties", list = list(), 
              type = "compartment_property"), list = list(new("compartment", 
              name = "DEPOT", index = 1L), new("compartment", name = "CENTRAL", 
              index = 2L), new("compartment", name = "PERIPHERAL", 
              index = 3L), new("compartment", name = "OUTPUT", 
              index = 4L)), type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateScenarios.R:69:3'): Simulate scenarios - make few changes on model (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "KA", rhs = "THETA_KA*exp(ETA_KA)", 
              comment = NA_character_), new("equation", lhs = "CL", 
              rhs = "THETA_CL*exp(ETA_CL)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "V3", 
                  rhs = "THETA_V3*exp(ETA_V3)", comment = NA_character_), 
              new("equation", lhs = "Q", rhs = "THETA_Q*exp(ETA_Q)", 
                  comment = NA_character_), new("equation", lhs = "S2", 
                  rhs = "V2", comment = NA_character_)), type = "model_statement")), 
          new("ode_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("ode", lhs = "A_DEPOT", rhs = "-KA*A_DEPOT", 
                  comment = NA_character_), new("ode", lhs = "A_CENTRAL", 
                  rhs = "KA*A_DEPOT + Q*A_PERIPHERAL/V3 + (-CL/V2 - Q/V2)*A_CENTRAL", 
                  comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
                  rhs = "-Q*A_PERIPHERAL/V3 + Q*A_CENTRAL/V2", 
                  comment = NA_character_), new("ode", lhs = "A_OUTPUT", 
                  rhs = "CL*A_CENTRAL/V2", comment = NA_character_), 
                  new("equation", lhs = "F", rhs = "A_CENTRAL/S2", 
                    comment = NA_character_)), type = "model_statement")), 
          new("error_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("equation", lhs = "CP", rhs = "F", 
                  comment = NA_character_), new("equation", lhs = "OBS_CP", 
                  rhs = "CP*(EPS_PROP + 1)", comment = NA_character_), 
                  new("equation", lhs = "Y", rhs = "OBS_CP", comment = NA_character_)), 
              type = "model_statement"))), type = "code_record"), 
      parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
      0L)), list = list(new("theta", unit = NA_character_, name = "KA", 
          index = 1L, value = 1, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "CL", index = 2L, value = 5, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V2", index = 3L, value = 80, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V3", index = 4L, value = 20, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "Q", index = 5L, value = 4, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 1L, 
          type = "var", name = "KA", index = 1L, value = 0.025, 
          fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 2L, type = "var", name = "CL", 
              index = 2L, value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 3L, type = "var", name = "V2", index = 3L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 4L, type = "var", name = "V3", index = 4L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 5L, type = "var", name = "Q", index = 5L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("sigma", index2 = 1L, 
              type = "var", name = "PROP", index = 1L, value = 0.025, 
              fix = FALSE, label = NA_character_, comment = NA_character_)), 
          type = "parameter"), compartments = new("compartments", 
          properties = new("compartment_properties", list = list(), 
              type = "compartment_property"), list = list(new("compartment", 
              name = "DEPOT", index = 1L), new("compartment", name = "CENTRAL", 
              index = 2L), new("compartment", name = "PERIPHERAL", 
              index = 3L), new("compartment", name = "OUTPUT", 
              index = 4L)), type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateTSLD.R:27:3'): Weibull model simulation works as expected (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "TVTAU", rhs = "THETA_TAU", 
              comment = NA_character_), new("equation", lhs = "TVSH", 
              rhs = "THETA_SH", comment = NA_character_), new("equation", 
              lhs = "TVCL", rhs = "THETA_CL", comment = NA_character_), 
              new("equation", lhs = "TVVC", rhs = "THETA_VC", comment = NA_character_), 
              new("equation", lhs = "TVPROP_RUV", rhs = "THETA_PROP_RUV", 
                  comment = NA_character_), new("line_break", comment = NA_character_), 
              new("equation", lhs = "TAU", rhs = "TVTAU * exp(ETA_TAU)", 
                  comment = NA_character_), new("equation", lhs = "SH", 
                  rhs = "TVSH * exp(ETA_SH)", comment = NA_character_), 
              new("equation", lhs = "CL", rhs = "TVCL * exp(ETA_CL)", 
                  comment = NA_character_), new("equation", lhs = "VC", 
                  rhs = "TVVC * exp(ETA_VC)", comment = NA_character_), 
              new("equation", lhs = "PROP_RUV", rhs = "TVPROP_RUV", 
                  comment = NA_character_)), type = "model_statement")), 
          new("ode_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("equation", lhs = "WF", rhs = "DOSE * exp(-pow((t-TDOS)/TAU, SH)) * (SH/TAU) * pow((t-TDOS)/TAU, SH-1)", 
                  comment = NA_character_), new("ode", lhs = "A_ABS", 
                  rhs = "-WF*A_ABS", comment = NA_character_), 
                  new("ode", lhs = "A_CENTRAL", rhs = "WF*A_ABS - CL/VC*A_CENTRAL", 
                    comment = NA_character_)), type = "model_statement")), 
          new("error_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("equation", lhs = "CONC", rhs = "A_CENTRAL/VC", 
                  comment = NA_character_), new("equation", lhs = "IPRED", 
                  rhs = "log(0.001)", comment = NA_character_), 
                  new("if_statement", condition = "CONC > 0", equation = new("equation", 
                    lhs = "IPRED", rhs = "log(CONC)", comment = NA_character_), 
                    comment = NA_character_), new("equation", lhs = "W", 
                    rhs = "PROP_RUV", comment = NA_character_), 
                  new("equation", lhs = "Y", rhs = "IPRED + W*EPS_RUV_FIX", 
                    comment = NA_character_)), type = "model_statement"))), 
      type = "code_record"), parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
  0L)), list = list(new("theta", unit = NA_character_, name = "TAU", 
      index = 1L, value = 50, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("theta", unit = NA_character_, 
      name = "SH", index = 2L, value = 3, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("theta", unit = NA_character_, 
      name = "CL", index = 3L, value = 3, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("theta", unit = NA_character_, 
      name = "VC", index = 4L, value = 60, fix = FALSE, label = NA_character_, 
      comment = NA_character_), new("theta", unit = NA_character_, 
      name = "PROP_RUV", index = 5L, value = 0.15, fix = FALSE, 
      label = NA_character_, comment = NA_character_), new("omega", 
      same = NA, index2 = 1L, type = "cv%", name = "TAU", index = 1L, 
      value = 20, fix = TRUE, label = NA_character_, comment = NA_character_), 
      new("omega", same = NA, index2 = 2L, type = "cv%", name = "SH", 
          index = 2L, value = 5, fix = TRUE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 3L, 
          type = "cv%", name = "CL", index = 3L, value = 20, fix = TRUE, 
          label = NA_character_, comment = NA_character_), new("omega", 
          same = NA, index2 = 4L, type = "cv%", name = "VC", index = 4L, 
          value = 20, fix = TRUE, label = NA_character_, comment = NA_character_), 
      new("sigma", index2 = 1L, type = "var", name = "RUV_FIX", 
          index = 1L, value = 1, fix = TRUE, label = NA_character_, 
          comment = NA_character_)), type = "parameter"), compartments = new("compartments", 
      properties = new("compartment_properties", list = list(), 
          type = "compartment_property"), list = list(new("compartment", 
          name = "ABS", index = 1L), new("compartment", name = "CENTRAL", 
          index = 2L)), type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateTable.R:23:3'): Simulate a bolus using the tabular dataset (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "KA", rhs = "THETA_KA*exp(ETA_KA)", 
              comment = NA_character_), new("equation", lhs = "CL", 
              rhs = "THETA_CL*exp(ETA_CL)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "V3", 
                  rhs = "THETA_V3*exp(ETA_V3)", comment = NA_character_), 
              new("equation", lhs = "Q", rhs = "THETA_Q*exp(ETA_Q)", 
                  comment = NA_character_), new("equation", lhs = "S2", 
                  rhs = "V2", comment = NA_character_)), type = "model_statement")), 
          new("ode_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("ode", lhs = "A_DEPOT", rhs = "-KA*A_DEPOT", 
                  comment = NA_character_), new("ode", lhs = "A_CENTRAL", 
                  rhs = "KA*A_DEPOT + Q*A_PERIPHERAL/V3 + (-CL/V2 - Q/V2)*A_CENTRAL", 
                  comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
                  rhs = "-Q*A_PERIPHERAL/V3 + Q*A_CENTRAL/V2", 
                  comment = NA_character_), new("ode", lhs = "A_OUTPUT", 
                  rhs = "CL*A_CENTRAL/V2", comment = NA_character_), 
                  new("equation", lhs = "F", rhs = "A_CENTRAL/S2", 
                    comment = NA_character_)), type = "model_statement")), 
          new("error_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("equation", lhs = "CP", rhs = "F", 
                  comment = NA_character_), new("equation", lhs = "OBS_CP", 
                  rhs = "CP*(EPS_PROP + 1)", comment = NA_character_), 
                  new("equation", lhs = "Y", rhs = "OBS_CP", comment = NA_character_)), 
              type = "model_statement"))), type = "code_record"), 
      parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
      0L)), list = list(new("theta", unit = NA_character_, name = "KA", 
          index = 1L, value = 1, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "CL", index = 2L, value = 5, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V2", index = 3L, value = 80, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V3", index = 4L, value = 20, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "Q", index = 5L, value = 4, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 1L, 
          type = "var", name = "KA", index = 1L, value = 0.025, 
          fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 2L, type = "var", name = "CL", 
              index = 2L, value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 3L, type = "var", name = "V2", index = 3L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 4L, type = "var", name = "V3", index = 4L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 5L, type = "var", name = "Q", index = 5L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("sigma", index2 = 1L, 
              type = "var", name = "PROP", index = 1L, value = 0.025, 
              fix = FALSE, label = NA_character_, comment = NA_character_)), 
          type = "parameter"), compartments = new("compartments", 
          properties = new("compartment_properties", list = list(), 
              type = "compartment_property"), list = list(new("compartment", 
              name = "DEPOT", index = 1L), new("compartment", name = "CENTRAL", 
              index = 2L), new("compartment", name = "PERIPHERAL", 
              index = 3L), new("compartment", name = "OUTPUT", 
              index = 4L)), type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateTimeVaryingCovariates.R:24:3'): Body weight as a true time varying covariate (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "KA", rhs = "THETA_KA*exp(ETA_KA)", 
              comment = NA_character_), new("equation", lhs = "CL", 
              rhs = "THETA_CL*exp(ETA_CL)*pow(BW/70, 0.75)", comment = NA_character_), 
              new("equation", lhs = "V", rhs = "THETA_V*exp(ETA_V)", 
                  comment = NA_character_), new("equation", lhs = "S2", 
                  rhs = "V", comment = NA_character_)), type = "model_statement")), 
          new("ode_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("ode", lhs = "A_DEPOT", rhs = "-KA*A_DEPOT", 
                  comment = NA_character_), new("ode", lhs = "A_CENTRAL", 
                  rhs = "-CL*A_CENTRAL/V + KA*A_DEPOT", comment = NA_character_), 
                  new("ode", lhs = "A_OUTPUT", rhs = "CL*A_CENTRAL/V", 
                    comment = NA_character_), new("equation", lhs = "F", 
                    rhs = "A_CENTRAL/S2", comment = NA_character_)), 
              type = "model_statement")), new("error_record", comment = NA_character_, 
              statements = new("model_statements", list = list(
                  new("equation", lhs = "CP", rhs = "F", comment = NA_character_), 
                  new("equation", lhs = "OBS_CP", rhs = "CP*(EPS_PROP + 1)", 
                    comment = NA_character_), new("equation", lhs = "Y", 
                    rhs = "OBS_CP", comment = NA_character_)), 
                  type = "model_statement"))), type = "code_record"), 
      parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
      0L)), list = list(new("theta", unit = NA_character_, name = "KA", 
          index = 1L, value = 1, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "CL", index = 2L, value = 5, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V", index = 3L, value = 80, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 1L, 
          type = "var", name = "KA", index = 1L, value = 0.025, 
          fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 2L, type = "var", name = "CL", 
              index = 2L, value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 3L, type = "var", name = "V", index = 3L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("sigma", index2 = 1L, 
              type = "var", name = "PROP", index = 1L, value = 0.025, 
              fix = FALSE, label = NA_character_, comment = NA_character_)), 
          type = "parameter"), compartments = new("compartments", 
          properties = new("compartment_properties", list = list(), 
              type = "compartment_property"), list = list(new("compartment", 
              name = "DEPOT", index = 1L), new("compartment", name = "CENTRAL", 
              index = 2L), new("compartment", name = "OUTPUT", 
              index = 3L)), type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateTimeVaryingCovariates.R:57:3'): Body weight as a true time varying covariate, 2 arms, individual body weights (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "KA", rhs = "THETA_KA*exp(ETA_KA)", 
              comment = NA_character_), new("equation", lhs = "CL", 
              rhs = "THETA_CL*exp(ETA_CL)*pow(BW/70, 0.75)", comment = NA_character_), 
              new("equation", lhs = "V", rhs = "THETA_V*exp(ETA_V)", 
                  comment = NA_character_), new("equation", lhs = "S2", 
                  rhs = "V", comment = NA_character_)), type = "model_statement")), 
          new("ode_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("ode", lhs = "A_DEPOT", rhs = "-KA*A_DEPOT", 
                  comment = NA_character_), new("ode", lhs = "A_CENTRAL", 
                  rhs = "-CL*A_CENTRAL/V + KA*A_DEPOT", comment = NA_character_), 
                  new("ode", lhs = "A_OUTPUT", rhs = "CL*A_CENTRAL/V", 
                    comment = NA_character_), new("equation", lhs = "F", 
                    rhs = "A_CENTRAL/S2", comment = NA_character_)), 
              type = "model_statement")), new("error_record", comment = NA_character_, 
              statements = new("model_statements", list = list(
                  new("equation", lhs = "CP", rhs = "F", comment = NA_character_), 
                  new("equation", lhs = "OBS_CP", rhs = "CP*(EPS_PROP + 1)", 
                    comment = NA_character_), new("equation", lhs = "Y", 
                    rhs = "OBS_CP", comment = NA_character_)), 
                  type = "model_statement"))), type = "code_record"), 
      parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
      0L)), list = list(new("theta", unit = NA_character_, name = "KA", 
          index = 1L, value = 1, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "CL", index = 2L, value = 5, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V", index = 3L, value = 80, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 1L, 
          type = "var", name = "KA", index = 1L, value = 0, fix = FALSE, 
          label = NA_character_, comment = NA_character_), new("omega", 
          same = NA, index2 = 2L, type = "var", name = "CL", index = 2L, 
          value = 0, fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 3L, type = "var", name = "V", 
              index = 3L, value = 0, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("sigma", index2 = 1L, 
              type = "var", name = "PROP", index = 1L, value = 0.025, 
              fix = FALSE, label = NA_character_, comment = NA_character_)), 
          type = "parameter"), compartments = new("compartments", 
          properties = new("compartment_properties", list = list(), 
              type = "compartment_property"), list = list(new("compartment", 
              name = "DEPOT", index = 1L), new("compartment", name = "CENTRAL", 
              index = 2L), new("compartment", name = "OUTPUT", 
              index = 3L)), type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  ── Error ('testSimulateWeirdCases.R:34:3'): Simulate a bolus with single observation at time 0 (rxode2/mrgsolve) ──
  <purrr_error_indexed/rlang_error/error/condition>
  Error in `(function (.x, .y, .f, ..., .progress = FALSE) 
  {
      map2_("list", .x, .y, .f, ..., .progress = .progress)
  })(.x = list(new("campsis_model", model = new("code_records", 
      list = list(new("main_record", comment = NA_character_, statements = new("model_statements", 
          list = list(new("equation", lhs = "KA", rhs = "THETA_KA*exp(ETA_KA)", 
              comment = NA_character_), new("equation", lhs = "CL", 
              rhs = "THETA_CL*exp(ETA_CL)", comment = NA_character_), 
              new("equation", lhs = "V2", rhs = "THETA_V2*exp(ETA_V2)", 
                  comment = NA_character_), new("equation", lhs = "V3", 
                  rhs = "THETA_V3*exp(ETA_V3)", comment = NA_character_), 
              new("equation", lhs = "Q", rhs = "THETA_Q*exp(ETA_Q)", 
                  comment = NA_character_), new("equation", lhs = "S2", 
                  rhs = "V2", comment = NA_character_)), type = "model_statement")), 
          new("ode_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("ode", lhs = "A_DEPOT", rhs = "-KA*A_DEPOT", 
                  comment = NA_character_), new("ode", lhs = "A_CENTRAL", 
                  rhs = "KA*A_DEPOT + Q*A_PERIPHERAL/V3 + (-CL/V2 - Q/V2)*A_CENTRAL", 
                  comment = NA_character_), new("ode", lhs = "A_PERIPHERAL", 
                  rhs = "-Q*A_PERIPHERAL/V3 + Q*A_CENTRAL/V2", 
                  comment = NA_character_), new("ode", lhs = "A_OUTPUT", 
                  rhs = "CL*A_CENTRAL/V2", comment = NA_character_), 
                  new("equation", lhs = "F", rhs = "A_CENTRAL/S2", 
                    comment = NA_character_)), type = "model_statement")), 
          new("error_record", comment = NA_character_, statements = new("model_statements", 
              list = list(new("equation", lhs = "CP", rhs = "F", 
                  comment = NA_character_), new("equation", lhs = "OBS_CP", 
                  rhs = "CP*(EPS_PROP + 1)", comment = NA_character_), 
                  new("equation", lhs = "Y", rhs = "OBS_CP", comment = NA_character_)), 
              type = "model_statement"))), type = "code_record"), 
      parameters = new("parameters", varcov = structure(numeric(0), dim = c(0L, 
      0L)), list = list(new("theta", unit = NA_character_, name = "KA", 
          index = 1L, value = 1, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "CL", index = 2L, value = 5, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V2", index = 3L, value = 80, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "V3", index = 4L, value = 20, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("theta", unit = NA_character_, 
          name = "Q", index = 5L, value = 4, fix = FALSE, label = NA_character_, 
          comment = NA_character_), new("omega", same = NA, index2 = 1L, 
          type = "var", name = "KA", index = 1L, value = 0.025, 
          fix = FALSE, label = NA_character_, comment = NA_character_), 
          new("omega", same = NA, index2 = 2L, type = "var", name = "CL", 
              index = 2L, value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 3L, type = "var", name = "V2", index = 3L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 4L, type = "var", name = "V3", index = 4L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("omega", same = NA, 
              index2 = 5L, type = "var", name = "Q", index = 5L, 
              value = 0.025, fix = FALSE, label = NA_character_, 
              comment = NA_character_), new("sigma", index2 = 1L, 
              type = "var", name = "PROP", index = 1L, value = 0.025, 
              fix = FALSE, label = NA_character_, comment = NA_character_)), 
          type = "parameter"), compartments = new("compartments", 
          properties = new("compartment_properties", list = list(), 
              type = "compartment_property"), list = list(new("compartment", 
              name = "DEPOT", index = 1L), new("compartment", name = "CENTRAL", 
              index = 2L), new("compartment", name = "PERIPHERAL", 
              index = 3L), new("compartment", name = "OUTPUT", 
              index = 4L)), type = "compartment"))), .y = 1L, .f = function (...) 
  {
      NULL
      NULL
      ...furrr_out <- ...furrr_fn(...)
      ...furrr_out
  })`: ℹ In index: 1.
  Caused by error:
  ℹ In index: 1.
  Caused by error:
  ! .onLoad failed in loadNamespace() for 'rxode2', details:
    call: NULL
    error: rxode2 compiled with rxode2parse with a different solving structure
  can try: install.packages('rxode2', type='source')
  
  [ FAIL 39 | WARN 0 | SKIP 29 | PASS 389 ]
  Deleting unused snapshots:
  • testDefaultPlots/scatterplot-colour-scenario.svg
  • testDefaultPlots/scatterplot-no-colour.svg
  • testDefaultPlots/shadedplot-colour-arm-strat-wt.svg
  • testDefaultPlots/shadedplot-colour-arm-wt-strat-wt.svg
  • testDefaultPlots/shadedplot-colour-scenario.svg
  • testDefaultPlots/spaghettiplot-colour-arm-strat-wt.svg
  • testDefaultPlots/spaghettiplot-colour-scenario.svg
  Error: Test failures
  Execution halted