sim.ratebystate function

Conduct simulation of state dependent rate variation