sysuse nlsw88, clear gen lnwage = ln(wage) asl_norm lnwage if union < ., dh reps(19999) nodots