capt log close
log using 771ps1,replace
set more off
clear

* EC 771 S2004 PS 1

use http://fmwww.bc.edu/ec-p/data/wooldridge/discrim
describe
tabulate state
summarize wagest
ttest income, by(state)
ttest pentree==1.39

capt program drop pinar_lf1
program define pinar_lf1
   version 8.0
   args lnf mu sigma
   quietly replace `lnf'=ln(normden($ML_y1, `mu', `sigma'))
end
ml model lf pinar_lf1 (mu: pfries=income prpblck) (sigma:)
* ml check
ml maximize

ml model lf pinar_lf1 (mu: pfries=income prpblck) (sigma:) if state==1
* ml check
ml maximize

regress pfries income prpblck

bootstrap " regress pfries income prpblck" _b, reps(1000)

set more on
log close