#Number of bootstrap samples
boot_samples = 5000
z_dist_boot = np.zeros(200)
z_dist_ideal_boot = np.zeros(200)
z_dist_min_boot = np.zeros(200)
Nsample = boot_samples
L0,L1 = df.get_limits(Nsample, size, rank)
#Actually run the methods and save their outputs
for obs_num in xrange(L0,L1):
boot_phi = np.random.choice(sampled_data,size=num_obs,replace=True)
z_dist_max, ideal_chisq_stop, min_chisq = lucy_replicator(50000,boot_phi)
z_dist_ideal_chisq = lucy_replicator(ideal_chisq_stop,boot_phi)[0]
z_dist_min_chisq = lucy_replicator(min_chisq,boot_phi)[0]
z_dist_boot = z_dist_boot+z_dist_max
z_dist_ideal_boot = z_dist_ideal_boot+z_dist_ideal_chisq
z_dist_min_boot = z_dist_min_boot+z_dist_min_chisq
