求助在STATA中怎么将两个不同数据的其中一个变量的和密度曲线做差后画图？

本人stata小白，目前在做收入差异的分位数回归分析，STATA有两个数据集，一个是URBAN.dat 一个是RURAL.dat，属于横截面数据且observation 数目不相等，应该无法merge。 我分别画出了城市和农村收入变量的核密度曲线，但是我想画一个图反映城乡收入核密度曲线在不同分位上的差异，也就是两组收入的核密度估计做差画图， 不知道该如何操作？是否只有数据都在一个数据集中才可以做核密度曲线的比较？

第二个问题，我用分位数回归命令 qui bsqreg wage x , q(.5) reps(.)  grqreg, cons ci ols olsci  做出了城市和农村两组影响收入的解释变量的分位数回归系数变化图，但是请问如何做城乡间某项同一个分位数回归系数的差异变化图？如在STATA中可以实现，可否详细打出实现代码？谢谢各位大神了。

求大神