Simulation code - main file

This code runs simulations presented in the paper. Note, this script requires a function found in simulate_func.R (accessed with source() at the top of the script).