R/simulation.r

In Thibaudpmx/VirtualTumor: What the Package Does (One Line, Title Case)

# model_equations <- function(){
#
#   test <-  exists("model_RxODE")
#
#   if(! test){
#
#
#
#
#   }
#
#   return(model_RxODE)
#
# }


# source("D:/these/Second_project/QSP/VirtualTumor/inst/model.r")


# simulations(add_events = add_events)

#' Title
#'
#' @param ind_param
#' @param add_events
#'
#' @return
#' @export
#'
#' @examples
#'
# simulations(ind_param = celltheque %>% slice(1), add_events = add_events) %>%
#   ggplot()+
#   geom_line(aes(time, Bcl2_I))
simulations <- function(ind_param = double(), add_events = tibble(), returnSim = T){



events <-  tibble(cmt = names(initial_cmt_values)[[1]], time = 0, amt = 0) %>%
    bind_rows(add_events) %>%
    mutate(evid = 1) %>%
    bind_rows(crossing(time = times, evid = 0, cmt = c(NA)))



parameter <- parameters_default_values
# colnames(parameter) <- names(parameters_default_values)

# ind_param <- c(ke_Venetoclax = 0.1)
# ind_param[]

for(a in names(ind_param)[map_lgl(ind_param, ~ is.numeric(.x))]){

  parameter[a] <- ind_param[[a]]

}

# initial_cmt_values
states <- initial_cmt_values


# evaluate the character
for(a in 1:length(states)){

states[a] <- with(data = as.tibble(parameter) %>%
                    mutate(rowname = names(parameter)) %>%
                    spread(key = "rowname", value = "value"), expr = eval(parse_expr(states[a])) )

}

states <- as.double(states)
names(states) <- names(initial_cmt_values)


res <- as_tibble(model_RxODE$solve(parameter, events, states))

if(returnSim == T ) return(res)

eval(criteria)

}