cleaned up helper functions
rewrote the whole simulation from scratch, but just for the simple
symmetric setting, also added code to compare our approach to that of
sugitani et al. (2013)
sugitani provided his simulation code
implemented the REAL sugitani rule I (though sans interim selection
based on final statistic) added code to compute any power and how
often a treatment gets dropped
tuned simulation parameters, added knitr script to produce a latex table
Finalized simulations and set up table