program define mylogit
  args lnf theta
  quietly replace `lnf' = ln(exp(`theta')/(1 + exp(`theta'))) if $ML_y1==1
  quietly replace `lnf' = ln(1/(1 + exp(`theta'))) if $ML_y1==0
end

ml model lf mylogit (voted2004 = yrseduc)