Our Stata code is optimized for efficiency