clear
set obs 25
generate e=invnorm(uniform())
generate x0=10*uniform()
generate u1=invnorm(uniform())
generate u2=invnorm(uniform())
generate x1=x0+u1
generate x2=x0+u2
generate y=x1+5*e
regress y x1
regress y x2
regress y x1 x2