- Bug - the values are incorrect for cond2 in ans1. Wrong order.
- 1:nrow() in apply
- Global variables
- Could add the cond1 and cond2 when each data.frame is created, i.e. in sim()
- Separate names of cond1 and cond2 and the associated mu and rate. Combine.
Turn this into a function