Skip to contents

Main Modelling Functions

Modelling Core Functions

add_item()
Defining parameters that may be used in model calculations
add_reactevt()
Define the modifications to other events, costs, utilities, or other items affected by the occurrence of the event
add_tte()
Define events and the initial event time
modify_event()
Modify the time of existing events
modify_item()
Modify the value of existing items
modify_item_seq()
Modify the value of existing items
new_event()
Generate new events to be added to existing vector of events
run_sim()
Run the simulation
run_sim_parallel()
Run simulations in parallel mode (at the simulation level)

Modelling Auxiliary Functions

create_indicators()
Creates a vector of indicators (0 and 1) for sensitivity/DSA analysis
pick_psa()
Helper function to create a list with random draws or whenever a series of functions needs to be called. Can be implemented within pick_val_v.
pick_val_v()
Select which values should be applied in the corresponding loop for several values (vector or list).
replicate_profiles()
Replicate profiles data.frame

Modelling Summary Functions

ceac_des()
Calculate the cost-effectiveness acceptability curve (CEAC) for a DES model with a PSA result
evpi_des()
Calculate the Expected Value of Perfect Information (EVPI) for a DES model with a PSA result
summary_results_det()
Deterministic results for a specific treatment
summary_results_sens()
Summary of sensitivity outputs for a treatment
summary_results_sim()
Summary of PSA outputs for a treatment

Statistics Functions

Statistics Auxiliary Functions

draw_tte()
Draw a time to event from a list of parametric survival functions
luck_adj()
Perform luck adjustment

Statistics Distributions Functions

cond_dirichlet()
Calculate conditional dirichlet values
cond_mvn()
Calculate conditional multivariate normal values
rgamma_mse()
Draw from a gamma distribution based on mean and se
qgamma_mse()
Use quantiles from a gamma distribution based on mean and se
rpoisgamma()
Draw time to event (tte) from a Poisson or Poisson-Gamma (PG) Mixture/Negative Binomial (NB) Process
qbeta_mse()
Draw from a beta distribution based on mean and se (quantile)
qcond_exp()
Conditional quantile function for exponential distribution
qcond_gamma()
Conditional quantile function for gamma distribution
qcond_gompertz()
Quantile function for conditional Gompertz distribution (lower bound only)
qcond_llogis()
Conditional quantile function for loglogistic distribution
qcond_lnorm()
Conditional quantile function for lognormal distribution
qcond_norm()
Conditional quantile function for normal distribution
qcond_weibull()
Conditional quantile function for weibull distribution
rbeta_mse()
Draw from a beta distribution based on mean and se
rcond_gompertz()
Draw from a conditional Gompertz distribution (lower bound only)
rcond_gompertz_lu()
Draw from a Conditional Gompertz distribution (lower and upper bound)
rdirichlet()
Draw from a dirichlet distribution based on number of counts in transition. Adapted from brms::rdirichlet
rdirichlet_prob()
Draw from a dirichlet distribution based on mean transition probabilities and standard errors
pcond_gompertz()
Survival Probaility function for conditional Gompertz distribution (lower bound only)

Other Functions and Utilities

disc_cycle()
Cycle discounting
disc_cycle_v()
Cycle discounting for vectors
disc_instant()
Calculate instantaneous discounted costs or qalys
disc_instant_v()
Calculate instantaneous discounted costs or qalys for vectors
disc_ongoing()
Calculate discounted costs and qalys between events
disc_ongoing_v()
Calculate discounted costs and qalys between events for vectors
ast_as_list()
Transform a substituted expression to its Abstract Syntax Tree (AST) as a list
extract_elements_from_list()
Extracts items and events by looking into modify_item, modify_item_seq, modify_event and new_event
extract_from_reactions()
Extract all items and events and their interactions from the event reactions list
extract_psa_result()
Extract PSA results from a treatment
tte.df
Example TTE IPD data