我想画各个分位上的两个样本的变化率，不知如何写命令？


先前尝试过：
use 2010_total_2.dta, clear
gen wage_2009= wage if year== 2009
gen wage_1989= wage if year== 1989
gen lnwage_2009 = ln(wage_2009)
gen lnwage_1989=  ln(wage_1989)
gen lnwage_gap = lnwage_2009 - lnwage_1989
quantile lnwage_gap
还尝试过
forvalues i= 1/99 {
by year, sort: egen w`i'= pctile(wage), p(`i')
gen lnw`i'= ln(w`i')
gen lnw`i'_2009 = lnw`i' if year== 2009
gen lnw`i'_1989 = lnw`i' if year== 1989

}
不过这样产生的数据还是需要贴入excel整理，不知道有没有简便的方法