I am conducting a sharp RD, with my running variable being the age of the woman in months. I am trying to calculate whether a particular legislation affected the age a woman gets married.
My dependent variable is age of marriage (agemar) and my running variable is age of the woman in months (agePCMAfunctionYrs), which is discrete since it is in months
Currently I am running the following command:
rdrobust agemar agePCMAfunctionYrs, c(15) covs(region religion) p(2)

However, because it is discrete I know I need to cluster the standard errors at the running variable level. I cannot see anywhere online how to do this - I would be very grateful if somebody would be able to say how I do this on STATA.