first stata code added!