R/split.R

In fracprolif: Fraction Proliferation via a Quiescent Growth Model

split_tot <- function(t, x0, q0, d1, d2, q1, q2, ad1, aq1, ad2, aq2, br=20)
{
  qbr <- prolif_q(br, x0, q0, d1, q1, ad1, aq1)
  dbr <- prolif_d(br, x0, q0, d1, q1, ad1, aq1)
  
  before <- prolif_tot(t,     x0,  q0, d1, q1, ad1, aq1)
  after  <- prolif_tot(t-br, dbr, qbr, d2, q2, ad2, aq2)
  
  answer       <- before
  answer[t>br] <- after[t>br]
  answer
}

split_q <- function(t, x0, q0, d1, d2, q1, q2, ad1, aq1, ad2, aq2, br=20)
{
  qbr <- prolif_q(br, x0, q0, d1, q1, ad1, aq1)
  dbr <- prolif_d(br, x0, q0, d1, q1, ad1, aq1)

  before <- prolif_q(t,     x0,  q0, d1, q1, ad1, aq1)
  after  <- prolif_q(t-br, dbr, qbr, d2, q2, ad2, aq2)

  answer       <- before
  answer[t>br] <- after[t>br]
  answer
}

split_d <- function(t, x0, q0, d1, d2, q1, q2, ad1, aq1, ad2, aq2, br=20)
{
  qbr <- prolif_q(br, x0, q0, d1, q1, ad1, aq1)
  dbr <- prolif_d(br, x0, q0, d1, q1, ad1, aq1)
  
  before <- prolif_d(t,     x0,  q0, d1, q1, ad1, aq1)
  after  <- prolif_d(t-br, dbr, qbr, d2, q2, ad2, aq2)

  answer       <- before
  answer[t>br] <- after[t>br]
  answer
}

split_fq <- function(t, x0, q0, d1, d2, q1, q2, ad1, aq1, ad2, aq2, br=20)
{
  split_q(t, x0, q0, d1, d2, q1, q2, ad1, aq1, ad2, aq2, br) /
  split_tot(t, x0, q0, d1, d2, q1, q2, ad1, aq1, ad2, aq2, br)
}

split_fd <- function(t, x0, q0, d1, d2, q1, q2, ad1, aq1, ad2, aq2, br=20)
{
  split_d(t, x0, q0, d1, d2, q1, q2, ad1, aq1, ad2, aq2, br) /
  split_tot(t, x0, q0, d1, d2, q1, q2, ad1, aq1, ad2, aq2, br)
}