I'm using the ssc packages rdrobust and rdplot to estimate and plot the fuzzy RD estimates of a regression of expenditure on population, used as running variable. Basically, the estimate obtained through rdrobust is positive, while when I plot the estimate from the apparently same regression using rdplot the jump is negative. I am using the same options as much as possible as the same kernel estimation method, same bandwidth, etc. However, the results do not match, which is weird given that the two commands come from the same package.

Do you have any idea why is this the case? I am also willing to use different packages or commands if you know some that work better for my case.