local b1 = ln(`or') generate logit1 = log( pgty1 / ( 1 - pgty1)) generate logit2 = log( pgty2 / ( 1 - pgty2))
