问题描述：在处理一个分位数回归，Y~ X1+X2+X3 这种模型，为了进一步讨论X1对Y的影响，让X1分别取0.1/0.5/0.9分位数，分别对应低中高三种状态，余下的X2和X3取0.5分位数，对Y进行条件密度预测。想绘制成如下的图形。

这种图形，可以使用MATLAB  或者R 来绘制，但是具体怎么绘制，却不知道怎么实现。有木有懂这块儿的热心人，能帮忙处理一下呢？