Added Stata precedure