R code and Stata DO file