如何在sas环境下做分位数回归

各位大神们！
        本人最近被导师要求写一篇分位数回归模型的东东，之前在网上搜到一些分位数回归的SAS过程步，但是还想多多了解一些细节的过程步，比如说如何做因变量和自变量的线性关系图（在各个分位数点）上，情况紧急，望各位大神赐教，不胜感激！！！

参考查一下sas的帮助文档，具体可以查proc quantreg的用法。以下是帮助文档里面的一个示例,
希望对你有帮助。
data a (drop=i);
   do i=1 to 1000;
      x1=rannor(1234);
      x2=rannor(1234);
      e=rannor(1234);
      if i > 950 then y=100 + 10*e;
      else y=10 + 5*x1 + 3*x2 + 0.5 * e;
      output;
   end;
run;
proc quantreg data=a;
   model y = x1 x2;
run;
proc quantreg algorithm=interior(tolerance=1e-6)
              ci=none data=a;
   model y = x1 x2 / itprint nosummary;
run;