*-
http://www.pinggu.org/bbs/thread-813133-1-1.html
use E:\stata答疑\graph\yellowriver_data\d_sample1.dta, clear
gen id = _n
reshape long salary, i(id) j(t)
tsset id t
label define gender_lbl 1 "female" 0 "male"
label value gender gender_lbl
label define edu_ba_lbl 1 "Yes_BA" 0 "No_BA"
label value edu_ba edu_ba_lbl
label define edu_ma_lbl 1 "Yes_MA" 0 "No_MA"
label value edu_ma edu_ma_lbl
label define edu_phd_lbl 1 "Yes_PhD" 0 "No_PhD"
label value edu_phd edu_phd_lbl
egen sex_ba = group(gender edu_ba) if edu_ba==1, label
egen sex_ma = group(gender edu_ma) if edu_ma==1, label
egen sex_phd = group(gender edu_phd) if edu_phd==1, label
gen g_score = irecode(entra_score, 80, 90)
label define g_score_lbl 0 "bad" 1 "good" 2 "excellent"
label value g_score g_score_lbl
bysort t g_score: egen avs = mean(salary)
bysort t sex_ba: egen avba = mean(salary)
bysort t sex_ma: egen avma = mean(salary)
bysort t sex_phd: egen avphd = mean(salary)
duplicates drop av* t g_score sex_*, force
twoway (connect avs t if g_score==0, sort) ///
(connect avs t if g_score==1) ///
(connect avs t if g_score==2) ///
(connect avba t if sex_ba==1) ///
(connect avba t if sex_ba==2) ///
(connect avma t if sex_ma==1) ///
(connect avma t if sex_ma==2) ///
(connect avphd t if sex_phd==1) ///
(connect avphd t if sex_phd==2) ///
, ///
legend(label(1 "Bad") label(2 "Good") label(3 "Excellent") ///
label(4 "BA-Male") label(5 "BA-FeMale") ///
label(6 "MA-Male") label(7 "MA-FeMale") ///
label(8 "PhD-Male") label(9 "PhD-FeMale") ///
rows(3)) ///
xlabel(1 "1-Year" 5 "5-Year" 10 "10-Year") ///
ylabel(,angle(0)) ytitle(Average Salary)