This script is designed to test the number of walkers and number of steps needed for the emcee sampler to converge to a result. yerr = 0.1 * y_true + verr * np.random.rand(N) y = a_true * x**2 + ...