R/simtime_from_backbone.R

Defines functions simtime_from_backbone

Documented in simtime_from_backbone

#' Determine simulation time from backbone
#' 
#' @param backbone A valid dyngen backbone object
#' @param burn Whether or not to compute the simtime for only the burn phase
#' 
#' @export
#' 
#' @return An estimation of the required simulation time
#' 
#' @examples 
#' backbone <- backbone_linear()
#' 
#' simtime_from_backbone(backbone)
#' 
#' model <- initialise_model(
#'   backbone = backbone,
#'   simulation_params = simulation_default(
#'     burn_time = simtime_from_backbone(backbone, burn = TRUE),
#'     total_time = simtime_from_backbone(backbone, burn = FALSE)
#'   )
#' )
simtime_from_backbone <- function(backbone, burn = FALSE) {
  start <- from <- `.` <- NULL
  
  exp_pat <- backbone$expression_patterns
  if (burn) exp_pat <- exp_pat %>% filter(burn)
  sim_time_sum <- exp_pat %>% filter(start) %>% pull(from) 
  sim_time_sum <- set_names(rep(0, length(sim_time_sum)), sim_time_sum)
  for (i in seq_len(nrow(exp_pat))) {
    sim_time_sum[[exp_pat$to[[i]]]] <- sim_time_sum[[exp_pat$from[[i]]]] + exp_pat$time[[i]]
  }
  total_time <- max(sim_time_sum * 1.2)
  total_time
}

Try the dyngen package in your browser

Any scripts or data that you put into this service are public.

dyngen documentation built on Oct. 12, 2022, 5:06 p.m.