def distributed_simulation(epsilon): y_star = np.array([[1, 282], [2, 20], [3, 13], [4, 4], [5, 2], [8, 1], [10, 1], [15, 1], [23, 1], [30, 1]], dtype=float) dim_theta = 3 if False: N_list = [50,100,200,500,1000] batches = 20 if True: N_list = [5] batches = 2 save_file = True type_sim = "mc_new" # construct the class tuberculosis_abc = loop_simulation(simulator, y_star, theta_sampler_mc, delta, epsilon, dim_theta, type_sim) # let it run ! tuberculosis_abc.loop_list_N(N_list,batches,save_file) print tuberculosis_abc.mean_theta_list print tuberculosis_abc.var_theta_list