Dear Stata Listers,

I'm using clustered standard errors in my RD design. For my linear split model I cannot figure out how to plot my RD with clustered standard errors displayed to read the discontinuity from the graph.
Code:
*linear split regression
reg y treatment_dummy running_variable treatmentdum_runningvar_interaction, cluster(running_variable)

This is the twoway scatter I can draw with robust standard errors, but I want to display clustered standard errors.
Code:
twoway (scatter y running_variable)                                     ///
    (lfitci y running_variable if treatment_dummy==0, ciplot(rline))    ///
    (lfitci y running_variable if treatment_dummy==1, ciplot(rline)),   ///
    xline(cutoff)

Also, cmogram and rdplot seem not to offer this option.

Would you know a way I can plot my regression discontinuity with clustered standard errors?


Thanks so much for your time and help!
Sipke Bontekoe