*! NJC 1.0.0 9 Nov 2007 program weibullfit_lf version 8.1 args lnf b c qui replace `lnf' = /// ln(`c') - ln(`b') + (`c' - 1) * (ln($S_MLy) - ln(`b')) - ($S_MLy / `b')^`c' end