program define spautosacnw version 11.0 args lf $spat_ARGS tempvar L1 L2 L3 Ro D0 D1 qui gen `D0'=0 qui gen `D1'=0 qui replace `D0' =1 if `e(depvar)' ==spat_llt qui replace `D1' =1 if `e(depvar)' > spat_llt qui gen double `Ro'=`Rho'*w1y_`e(depvar)' qui gen double `L1'=`Lambda'*spat_eigw qui gen double `L2'=`Lambda'*w1y_`e(depvar)' qui gen double `L3'= `Rho'*`Lambda'*w2y_`e(depvar)' qui forval i=1/$spat_kx { tempvar X`i' XX`i' qui gen double `X`i''=`beta`i''*`:word `i' of `:colnames(spat_ols)'' local XB "`XB'`X`i''-" qui gen double `XX`i''=`Lambda'*`beta`i''*spat_w1x_`i' local ZB "`ZB'`XX`i''+" } qui replace `lf' =`D1'*(((`e(depvar)'-`Ro'-`L2'+`L3'-`XB'`beta0'+`ZB'`Lambda'*`beta0')/`Sigma')-exp((`e(depvar)'-`Ro'-`L2'+`L3'-`XB'`beta0'+`ZB'`Lambda'*`beta0')/`Sigma')) /// -`D0'*exp((`Ro'+`L2'-`L3'+`XB'`beta0'-`ZB'`Lambda'*`beta0')/`Sigma')-`D1'*log(`Sigma') end