clear
set seed 13579
set obs 1000
gen x1=invnorm(uniform())
gen x2=2*invnorm(uniform())
gen e=invnorm(uniform())
gen y=3+x1+3*x2+e
gen z=y if y<10
replace z=10 if y>=10
*drop in 1/100
reg y x1 x2
tobit z x1 x2, ul(10) //eq1
tobit y x1 x2, ul(10) //eq2
sum y
/*
Variable | Obs Mean Std. Dev. Min Max
-------------+--------------------------------------------------------
y | 1000 2.792235 5.905514 -14.7648 24.75783
*/