clear
set obs 10000
generate age = runiformint(18, 99)
generate income = runiform(0, 250)
generate x = rlogistic()
generate wise = x < (age-55)*0.2 + (income-55)*0.1

logit wise age income

quietly margins, at(age=40 income=(0(10)100))
marginsplot, name(m40)
quietly margins, at(age=55 income=(0(10)100))
marginsplot, name(m55)
quietly margins, at(age=65 income=(0(10)100))
marginsplot, name(m65)