difference_in_means function

Design-based difference-in-means estimator