program define spautolagnw
version 10.0
args lnf XB rho sigma
tempvar grho spwy D0 D1
qui gen double `grho'=`rho'*spat_eigw
qui gen double `spwy'=`rho'*w1y_$ML_y1
tempvar D0 D1
qui gen `D0'=0
qui gen `D1'=0
qui replace `D0' =1 if $ML_y1 ==0
qui replace `D1' =1 if $ML_y1 >0
qui replace `lnf'= `D1'*((($ML_y1-`spwy'-`XB')/`sigma')-exp(($ML_y1-`spwy'-`XB')/`sigma')) ///
  -`D0'*exp(($ML_y1-`spwy'-`XB')/`sigma')-`D1'*ln(`sigma')
end