*The data files referenced in the first three lines of the following stata program may be
*downloaded from https://www.gesis.org/issp/home/

use "C:\Users\ecagj\ISSP\2014\ZA6670_v2-0-0.dta" 
append using "C:\Users\ecagj\ISSP\2013\ZA5950_v2-0-0.dta"
append using "C:\Users\ecagj\ISSP\2012\ZA5900_v2-0-0.dta"
keep year V1 C_ALPHAN AGE SEX EDUCYRS WRKHRS TW_RINC HR_RINC CZ_RINC DE_RINC DK_RINC ES_RINC FI_RINC FR_RINC GB_RINC IL_RINC IN_RINC IS_RINC JP_RINC KR_RINC LT_RINC NO_RINC PH_RINC RU_RINC SE_RINC SI_RINC SK_RINC TR_RINC US_RINC ZA_RINC CH_RINC
drop if WRKHRS==0|WRKHRS>80
drop if EDUCYRS>25|EDUCYRS<5|EDUCYRS==.
drop if WRKHRS==.
g AGE2=AGE*AGE
g y1=0
g y2=0
replace y1=1 if year==2013
replace y2=1 if year==2014
save "C:\Users\ecagj\ISSP\2012\pooled1012-14.dta", replace
g inc=TW_RINC if C_ALPHAN=="TW"
drop if inc>500000
replace inc=ln(inc/(WRKHRS*4)) if C_ALPHAN=="TW"
regr inc AGE AGE2 EDUCYRS SEX y1 y2
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.1)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.25)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.5)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.75)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.9)
use "C:\Users\ecagj\ISSP\2012\pooled1012-14.dta", replace
g inc=HR_RINC if C_ALPHAN=="HR"
drop if inc>20000
replace inc=ln(inc/(WRKHRS*4)) if C_ALPHAN=="HR"
regr inc AGE AGE2 EDUCYRS SEX y1 y2
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.1)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.25)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.5)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.75)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.9)
use "C:\Users\ecagj\ISSP\2012\pooled1012-14.dta", replace
g inc=CZ_RINC if C_ALPHAN=="CZ"
drop if inc>200000
replace inc=ln(inc/(WRKHRS*4)) if C_ALPHAN=="CZ"
regr inc AGE AGE2 EDUCYRS SEX y1 y2
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.1)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.25)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.5)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.75)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.9)
use "C:\Users\ecagj\ISSP\2012\pooled1012-14.dta", replace
drop if inc>800000
g inc=DK_RINC if C_ALPHAN=="DK"
replace inc=ln(inc/(WRKHRS*4)) if C_ALPHAN=="DK"
regr inc AGE AGE2 EDUCYRS SEX y1 y2
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.1)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.25)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.5)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.75)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.9)
use "C:\Users\ecagj\ISSP\2012\pooled1012-14.dta", replace
g inc=FI_RINC if C_ALPHAN=="FI"
drop if inc>100000
replace inc=ln(inc/(WRKHRS*4)) if C_ALPHAN=="FI"
regr inc AGE AGE2 EDUCYRS SEX y1 y2
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.1)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.25)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.5)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.75)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.9)
use "C:\Users\ecagj\ISSP\2012\pooled1012-14.dta", replace
g inc=FR_RINC if C_ALPHAN=="FR"
drop if inc>100000
replace inc=ln(inc/(WRKHRS*4)) if C_ALPHAN=="FR"
regr inc AGE AGE2 EDUCYRS SEX y1 y2
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.1)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.25)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.5)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.75)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.9)
use "C:\Users\ecagj\ISSP\2012\pooled1012-14.dta", replace
g inc=IS_RINC if C_ALPHAN=="IS"
drop if inc>5000000
replace inc=ln(inc/(WRKHRS*4)) if C_ALPHAN=="IS"
regr inc AGE AGE2 EDUCYRS SEX y1 y2
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.1)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.25)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.5)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.75)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.9)
use "C:\Users\ecagj\ISSP\2012\pooled1012-14.dta", replace
g inc=IN_RINC if C_ALPHAN=="IN"
drop if inc>200000
replace inc=ln(inc/(WRKHRS*4)) if C_ALPHAN=="IN"
regr inc AGE AGE2 EDUCYRS SEX y1 y2
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.1)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.25)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.5)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.75)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.9)
use "C:\Users\ecagj\ISSP\2012\pooled1012-14.dta", replace
g inc=JP_RINC if C_ALPHAN=="JP"
drop if inc>25000000
replace inc=ln(inc/(WRKHRS*4)) if C_ALPHAN=="JP"
regr inc AGE AGE2 EDUCYRS SEX y1 y2
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.1)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.25)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.5)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.75)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.9)
use "C:\Users\ecagj\ISSP\2012\pooled1012-14.dta", replace
g inc=KR_RINC if C_ALPHAN=="KR"
drop if inc>25000000
replace inc=ln(inc/(WRKHRS*4)) if C_ALPHAN=="KR"
regr inc AGE AGE2 EDUCYRS SEX y1 y2
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.1)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.25)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.5)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.75)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.9)
use "C:\Users\ecagj\ISSP\2012\pooled1012-14.dta", replace
g inc=LT_RINC if C_ALPHAN=="LT"
drop if inc>5000
replace inc=ln(inc/(WRKHRS*4)) if C_ALPHAN=="LT"
regr inc AGE AGE2 EDUCYRS SEX y1 y2
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.1)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.25)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.5)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.75)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.9)
use "C:\Users\ecagj\ISSP\2012\pooled1012-14.dta", replace
g inc=NO_RINC if C_ALPHAN=="NO"
drop if inc>10000000
replace inc=ln(inc/(WRKHRS*4)) if C_ALPHAN=="NO"
regr inc AGE AGE2 EDUCYRS SEX y1 y2
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.1)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.25)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.5)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.75)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.9)
use "C:\Users\ecagj\ISSP\2012\pooled1012-14.dta", replace
g inc=PH_RINC if C_ALPHAN=="PH"
drop if inc>200000
replace inc=ln(inc/(WRKHRS*4)) if C_ALPHAN=="PH"
regr inc AGE AGE2 EDUCYRS SEX y1 y2
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.1)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.25)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.5)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.75)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.9)
use "C:\Users\ecagj\ISSP\2012\pooled1012-14.dta", replace
g inc=RU_RINC if C_ALPHAN=="RU"
drop if inc>200000
replace inc=ln(inc/(WRKHRS*4)) if C_ALPHAN=="RU"
regr inc AGE AGE2 EDUCYRS SEX y1 y2
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.1)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.25)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.5)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.75)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.9)
use "C:\Users\ecagj\ISSP\2012\pooled1012-14.dta", replace
g inc=SK_RINC if C_ALPHAN=="SK"
drop if inc>5000
replace inc=ln(inc/(WRKHRS*4)) if C_ALPHAN=="SK"
regr inc AGE AGE2 EDUCYRS SEX y1 y2
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.1)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.25)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.5)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.75)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.9)
use "C:\Users\ecagj\ISSP\2012\pooled1012-14.dta", replace
g inc=SI_RINC if C_ALPHAN=="SI"
drop if inc>10000
replace inc=ln(inc/(WRKHRS*4)) if C_ALPHAN=="SI"
regr inc AGE AGE2 EDUCYRS SEX y1 y2
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.1)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.25)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.5)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.75)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.9)
use "C:\Users\ecagj\ISSP\2012\pooled1012-14.dta", replace
g inc=ZA_RINC if C_ALPHAN=="ZA"
drop if inc>100000
replace inc=ln(inc/(WRKHRS*4)) if C_ALPHAN=="ZA"
regr inc AGE AGE2 EDUCYRS SEX y1 y2
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.1)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.25)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.5)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.75)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.9)
use "C:\Users\ecagj\ISSP\2012\pooled1012-14.dta", replace
g inc=ES_RINC if C_ALPHAN=="ES"
drop if inc>10000
replace inc=ln(inc/(WRKHRS*4)) if C_ALPHAN=="ES"
regr inc AGE AGE2 EDUCYRS SEX y1 y2
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.1)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.25)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.5)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.75)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.9)
use "C:\Users\ecagj\ISSP\2012\pooled1012-14.dta", replace
g inc=SE_RINC if C_ALPHAN=="SE"
drop if inc>200000
replace inc=ln(inc/(WRKHRS*4)) if C_ALPHAN=="SE"
regr inc AGE AGE2 EDUCYRS SEX y1 y2
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.1)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.25)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.5)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.75)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.9)
use "C:\Users\ecagj\ISSP\2012\pooled1012-14.dta", replace
g inc=CH_RINC if C_ALPHAN=="CH"
drop if inc>20000
replace inc=ln(inc/(WRKHRS*4)) if C_ALPHAN=="CH"
regr inc AGE AGE2 EDUCYRS SEX y1 y2
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.1)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.25)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.5)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.75)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.9)
use "C:\Users\ecagj\ISSP\2012\pooled1012-14.dta", replace
g inc=TR_RINC if C_ALPHAN=="TR"
drop if inc>10000
replace inc=ln(inc/(WRKHRS*4)) if C_ALPHAN=="TR"
regr inc AGE AGE2 EDUCYRS SEX y1 y2
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.1)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.25)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.5)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.75)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.9)
use "C:\Users\ecagj\ISSP\2012\pooled1012-14.dta", replace
g inc=US_RINC if C_ALPHAN=="US"
drop if inc>200000
replace inc=ln(inc/(WRKHRS*4)) if C_ALPHAN=="US"
regr inc AGE AGE2 EDUCYRS SEX y1 y2
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.1)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.25)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.5)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.75)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.9)
use "C:\Users\ecagj\ISSP\2012\pooled1012-14.dta", replace
g inc=DE_RINC if C_ALPHAN=="DE-W"|C_ALPHAN=="DE-E"
drop if inc>100000
replace inc=ln(inc/(WRKHRS*4)) if C_ALPHAN=="DE-W"|C_ALPHAN=="DE-E"
regr inc AGE AGE2 EDUCYRS SEX y1 y2
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.1)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.25)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.5)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.75)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.9)
use "C:\Users\ecagj\ISSP\2012\pooled1012-14.dta", replace
g inc=GB_RINC if C_ALPHAN=="GB-GBN"
drop if inc>5000
replace inc=ln(inc/(WRKHRS*4)) if C_ALPHAN=="GB-GBN"
regr inc AGE AGE2 EDUCYRS SEX y1 y2
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.1)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.25)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.5)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.75)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.9)
use "C:\Users\ecagj\ISSP\2012\pooled1012-14.dta", replace
g inc=IL_RINC if C_ALPHAN=="IL"
drop if inc>20000
replace inc=ln(inc/(WRKHRS*4)) if C_ALPHAN=="IL"
regr inc AGE AGE2 EDUCYRS SEX y1 y2
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.1)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.25)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.5)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.75)
qreg inc AGE AGE2 EDUCYRS SEX y1 y2, q(.9)