clear
drawnorm error, n(1000)
gen x1=uniform()
gen x2=uniform()
gen x3=uniform()
histogram error
gen y=10*x1-5*x2+7*x3^2+error
reg y x1 x2 x3
gen z2=x2+1
reg y x1 x2 z2
replace z2=x2*x1
corr x1 x2 z2
reg y x1 x2 z2
g sex=1
replace sex=0 if x1>.3
replace y=10*x1-5*z2+7*x3^2+2*sex+error
reg y x1 x2 z2 sex
scatter y z2, by(sex)