# Simulate data
clear
set obs 100
set seed 321
matrix C = (1, .7 \ .7, 1)
corr2data x y, means(50 70) corr(C)
gen z=x
sum x y z
#Modify z variable
replace z=z+5 if y>70
replace z=z-5 if y<70
replace z=z+1 if x>51
replace z=z-1 if x<49
replace z=z-2 if x<51 & x>49 & z<50
replace z=z+2 if x<51 & x>49 & z>50
replace z=z+1 if x<49
replace z=z-3 if x>=50.4 & x<=50.6 & z>50
replace z=z-2 if x<=49.4 & z>50
replace z=z+3 if x>=50.4 & x<=50.6 & z<50
replace z=z+2 if x>=48.5 & x<=49 & z<50
replace z=z+2 if x>=48.1 & x<=48.5 & z<50