Stata code goes here