timesteps = generate_timesteps(10, 0.8, [0.001, 0.5])
model = MCMC(timesteps)
model.sample(10000)
[np.mean(model.trace(t).gettrace()) for t in timesteps]
How can I make a discrete state Markov model with pymc?