Estimates average treatment effects at the cutoff based on sharp
regression discontinuity designs (RDD) and multiple imputation regression
discontinuity designs (MIRDD). It provides diagnostic tools for RDD by
comparing results with those from MIRDD, as proposed in Takahashi (2023)