ivtob_1_lf.ado

是一个经济学管理应用软件 很难找的 但是经济学学生又必须用到

*! version 1.0.2  28jun2004

program ivtob_1_lf

	version 8.0
	args lnf xb1 xb2 alpha lns lnv
	
	tempvar s v
	qui gen double `s' = exp(`lns')
	qui gen double `v' = exp(`lnv')
	
	tempvar eq2resid w
	qui gen double `eq2resid' = $ML_y2 - `xb2'
	qui gen double `w' = `xb1' + `eq2resid'*`alpha'
	
	qui replace `lnf' = -0.5*ln(2*_pi) - 0.5*ln(`v'^2) - ///
		`eq2resid'^2 / (2*`v'^2)
	qui replace `lnf' = `lnf' - 0.5*ln(2*_pi) - 0.5*ln(`s'^2) - ///
		(($ML_y1-`w')^2 / (2*`s'^2)) ///
		if $ML_y1 > $IVT_ll & $ML_y1 < $IVT_ul
	qui replace `lnf' = `lnf' + ln(1 - norm(`w' / `s')) ///
		if $ML_y1 <= $IVT_ll
	qui replace `lnf' = `lnf' + ln(norm(`w'/`s')) if $ML_y1 >= $IVT_ul
	
end