用stata10学习MLE过程中,曾发现是在写myprobit_lf.ado程序时,要注意`xb'的写法,否则会有invalid name的错误提示。
——————————————————————————————————myprobit_lf.ado
program myprobit_lf
version 10.0
args lnf xb
quietly replace `lnf' = ln(normal(`xb')) if $ML_y1==1
quietly replace `lnf'= ln(normal(-`xb')) if $ML_y1==0
end
——————————————————————————————————
ml model lf myprobit_lf (died = drug2 drug3 age )
ml maximize