Hmm this is really bad, I realize this is only possible from Scipy 0.16 onwards... Before that, reproducible sampling for Scipy distributions is only possible by resetting the global Numpy random state. Yuck. 