我刚才用一个笨方法,您给看一下。
use 2010_total_2.dta, clear
keep if year== 2009 | year== 1989 // 只保留需要用的两年的样本
forvalues i= 1/99 {
by year, sort: egen w`i'= pctile(wage), p(`i')
gen lnw`i'= ln(w`i')
}
forvalues i= 1/99 {
duplicates drop year w`i', force
}
forvalues i= 1/99 {
drop w`i'
}
xpose, clear
gen ratio = v2 - v1
drop if v1 == 1989
gen n= _n
twoway line ratio n, sort lw(*1.1) msize(large) ///
xlabel(0(10)100) ylabel(-1(1)3) ///
xtitle("分位") ytitle(" ") lw(*2)