A supremely faster version of a function from dtpcrm v0.1.1 — stop_for_excess_toxicity

Originally, the sampling in stats::rnorm() (see inline comments in code) consumed 53% of run-time in a 6-cohort VIOLA DTP. After this change, it doesn't even show up! More importantly, the consumption is now dominated by (at 75%) by the objective function 'f' in integrate().

stop_for_excess_toxicity_empiric(
  x,
  tox_lim,
  prob_cert,
  dose = 1,
  suppress_dose = TRUE
)

Arguments

x: A object of class mtd
tox_lim: Scalar upper threshold on estimated toxicity rate
prob_cert: Confidence level for threshold exceedance
dose: Integer scalar, the dose being considered
suppress_dose: Logical; if TRUE the MTD is set to NA when trial stop is recommended.

Value

The mtd object x, with stop decision annotated

Author

Adapted by David C. Norris from original dtpcrm::stop_for_excess_toxicity_empiric