Stata code to produce the analysis results