请问哪位大侠那里有stata里做分位数分解的程序？-经管之家官网！

分位数回归分解第二步：边际工资分布，第三步反事实分析怎么用stata实现啊？
2.从分位数回归系数得到边际工资分布的具体步骤
( a) 从均匀分布 U［0，1］中随机抽取一个样本容量为 l 的样本，u1，u2，…ul;
( b) 在女性样本中，分别以 θ = u1，u2…ul 做分位数回归，
Qθ( y| xf) = xf' β( θ)
得到 l 个分位数回归的系数向量 β
∧
f( ui) ，i = 1，2…l;
( c) 从女性被解释变量中有放回地重复随机抽样，抽取一个样本容量为 l 的样本，表示为
{xi* ( f) }，i =1，2…l。
( d) 把( b) 中得到的系数和( c) 中得到的样本相结合，得到一个新的样本
{yi* ( f) }≡{xi* ( f) β
∧
f( ui) } il = 1
它便是我们需要的边际分布样本。
3． 反事实分析具体步骤:
( a) 利用前面的步骤( a) —( d) 得到女性对数工资边际分布的样本{y*f };
( b) 利用男性样本计算小学人群所占的比例 rm( 1) ;
( c) 从{yf* }中分离出小学人群样本，并从中进行有放回地重复随机抽样，抽出一个样本容量为 l × rm( 1) 的样本。
( d) 对初中、高中和大学人群，重复步骤( b) 和步骤( c) ，得到其他一些类似样本。
( e) 合并所有得到的样本，就可以得到女性教育水平同男性相同时女性工资的边际分布。
只剩两个币了，能解决问题的大侠我在追加悬赏哦~~