// Open Allbus data
use V388 V151 using "ZA4602_v1-0-0.dta", clear
// Income
recode V388 (99997 = .a "Refused") (99999 = .b "No answer"), gen(personal_income)
label var personal_income "Personal net income"
// Sex
recode V151 (1 = 0 "Men") (2 = 1 "Women"), gen(sex)
label var sex "Gender"
vioplot personal_income, over(sex) ytitle("Personal net income in Euro per month")
