*! mynormal_lf_d0_c1 v1.0.0 CFBaum 11aug2008 program mynormal_lf_d0_c1 version 10.1 args todo b lnf tempvar xb mu tempname a lnsigma mleval `a' = `b', eq(1) scalar mleval `xb' = `b', eq(2) mleval `lnsigma' = `b', eq(3) scalar quietly { generate double `mu' = `xb' - exp(`a')* $x1 mlsum `lnf' = ln( normalden($ML_y1,`mu', exp(`lnsigma')) ) } end