use dlsy_khb.dta
khb logit univ fses || abil

capture drop fses2
generate fses2 = fses + rnormal(0, 0.15)
pwcorr fses fses2
khb logit univ fses || fses2, summary

khb logit univ fses2 || fses, summary

drawnorm x1 x2, corr(1, .5, 1) cstorage(l)