I am using optimal bandwidth selection using STATA command rdbwselect. I am using the following command, however, this does not work when I include "fuzzy".

rdbwselect edyrtotal birthcent, c(0) p(1) kernel(triangular) bwselect(msetwo) fuzzy(edyrtotal) covs(wealthq urban) masspoints(check)

I would like to find a jump before and after of the education reform which should increase the year to complete elementary school by one year in Kenya.

edyrtotal= education year total
birthcent= Year of birth - 1971 (when the reform was implemented)
wealthq=wealth index、urban=type of residence

When I delete the fuzzy command (In other words, when trying sharp RDD), it seems work fine without any error. But the reform did not affect all people, so I think it should be Fuzzy RD.

I got the following error although I am using Fuzzy command.
-----------------------------------------------------------------------------------------------
Mass points detected in the running variable.
Try using option masspoints(adjust)
Not enough variability to compute the preliminary bandwidth. Try checking for mass points with option masspoints(check).
Not enough variability to compute the bias bandwidth (b). Try checking for masspoints with option masspoints(check).
Not enough variability to compute the loc. poly. bandwidth (h). Try checking for mass points with option masspoints(check).
-----------------------------------------------------------------------------------------------

Could someone advise?
I am using STATA 17

Thank you very much in advance.