A comparison between threshold ergodic sets and stochastic simulation of boolean networks for modelling cell differentiation