R/pk_2cmt_bolus_dose_from_cmin.R

In ronkeizer/PKmisc: Functions for clinical pharmacokinetics

#' Calculate dose to achieve steady state trough
#' for 2-compartmental PK model bolus dosing at steady state
#'
#' @param cmin desired trough concentration
#' @param tau dosing interval
#' @param CL clearance
#' @param V volume of distribution
#' @param Q inter-compartimental clearance
#' @param V2 volume of peripheral compartment
#' @export
pk_2cmt_bolus_dose_from_cmin <- function(cmin = 1,
                                         tau = 12,
                                         CL = 3, V = 30,
                                         Q = 2, V2 = 20) {

  # reparametrization:
  terms <- (Q/V) + (Q/V2) + (CL/V)
  beta <- 0.5 * (terms - sqrt(terms^2 - 4*(Q/V2)*(CL/V)))
  alpha <- ((Q/V2)*(CL/V))/beta
  A <- (1/V) * (alpha - (Q/V2))/(alpha-beta)
  B <- (1/V) * ((beta - Q/V2)/(beta-alpha))

  ## equation:
  dose <- cmin / (A*exp(-alpha*tau)/(1-exp(-alpha*tau)) + B*exp(-beta*tau)/(1-exp(-beta*tau)))
  return(dose)
}