Average Treatment Effects estimation by simulation

The parallelization is commented out (#)