The R code for performing the simulations described in the paper is included in "simulation_code.R" and the results are included in "simulation_400.RData" (for n=400) and
"simulation_100.RData" (for n=100). Required packages are MASS and snowfall, where the number of CPUs may be adjusted depending on computational resource. To run the
simulation, specify the working directory and run all code.

The file "functions.R" contains all functions necessary for performing the simulations. The command for sourcing the functions are included within each function in 
"simulation_code.R."



