Implementation of the bunching estimator for kinks and notches.
Allows for flexible estimation of counterfactual (e.g. controlling for round number bunching, accounting for other bunching masses within bunching window, fixing bunching point to be minimum, maximum or median value in its bin, etc.).
It produces publication-ready plots in the style followed since Chetty et al. (2011)