R/run_markov_simulation.R

In twig: For Streamlining Decision and Economic Evaluation Models using Grammar of Modeling

run_markov_simulation <- function(sim, twig_list, verbose = FALSE, offset_trace_cycle = 1){

    with(twig_list, {

        eval_funs <- evaluate_functions(sim, fun_core_df, fun_sim_args, prob_payoff_funs, params, arg_value_sizes, fun_args)

        F_mat <- evaluate_fun_sim(F0, IDX, prob_funs, eval_funs)

        E <- get_E(E0, F_mat, non_compl_id, event_prob_link, hash_id, compl_id)

        A <- get_A(A0_idx, E, A_idx, paths, n_paths)

        P_array <- calculate_transition_probs(P0_mat, A, dest_paths, unique_non_current_dest, dim_P, dimnames_P, is_cycle_dep, unique_dest_names, p_stay)

        T_array <- create_trace_array(arg_value_sizes, arg_values, p0_array[,,sim], 
            P_array, sim, is_cycle_dep, n_decisions, n_cycles)

    R_array <- calculate_payoffs(sim, R0_array, event_indep_payoffs, eval_funs, R_non_event_dep_idx, 
        IDX_path_dep, event_dep_payoffs, A, payoff_funs, dimnames_R0, size_core_non_event_args,
    n_cycles, is_cycle_dep)

    R_array_cycle <- return_R_array_cycle(R_array, payoff_funs, T_array, array_discount, n_cycles, offset_trace_cycle = offset_trace_cycle)

  R_sim <- apply(R_array_cycle, c(3,4), sum)
  if (verbose){

    evaluated_funs <- get_eval_funs_list(eval_funs, fun_core_df, twig_funs) 

    sim_results <- list(sim = sim,
                        
                        evaluated_funs = evaluated_funs,
                        evaluated_prob_funs_combined = F_mat, 
                        
                        path_event_options = NULL, 
                        path_probs = A, 
                        event_probs = E,
                        
                        markov_trans_probs = P_array, 
                        markov_trace = T_array, 

                        cycle_payoffs = R_array,
                        cycle_ev = R_array_cycle, 
                        
                        sim_ev = R_sim, 
                        mean_ev = R_sim
                        )
    
    
    return(sim_results)
  } else {
    return(R_sim)
  }

  gc()
  }) 

  }