请问大家,如下是先写evaluate program,然后用ml来最大化估计logit模型的代码。为什么 loglikelihood一开始就用
replace替代,难道程序前面的部分已经建立了`lnf'这个变量?另外,即使去掉quietly,在数据里也见不到生成的变量,这是怎么回事?use http://www.ats.ucla.edu/stat/data/hsbdemo
capture program drop lflogitprogram lflogit  
version 10.1  
args lnf xb  local y "$ML_y1"  
quietly replace `lnf' = ln(  invlogit(`xb')) if `y'==1  
quietly replace `lnf' = ln(1-invlogit(`xb')) if `y'==0
end 
ml model lf lflogit (honors = read female) 
ml maximize