*! version 1.6.0 MLB 12Sep2012 program propcnslogit_lf version 8.2 args lnf unconstrained constrained lambda tempvar theta quietly{ gen double `theta' = `unconstrained' + (`lambda'*`constrained') replace `lnf' = ln(invlogit( `theta')) if $ML_y1 == 1 replace `lnf' = ln(invlogit(-`theta')) if $ML_y1 == 0 } end