How to run the difference in difference estimates for binary outcomes?