mixed_ll_reparm.ado

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

*! version 1.0.5  1dec2004
program mixed_ll_reparm
	version 8.2
	args todo b lnf
	
	tempname d b2 err1 err2 s2
	local p = colsof(`b')
	scalar `s2' = exp(2*`b'[1,`p'])
	mata: _xtm_gamma_to_theta("`b2'","`b'","`err1'")
	if scalar(`err1') == 1 {
		scalar `lnf' = .
	}
	else {
		mata: _xtm_theta_to_delta(`"`d'"',`"`b2'"',"`err2'")
		if scalar(`err2') == 1 {
			scalar `lnf' = .
		}
		else {
			mata: _xtm_mixed_ll(`"`d'"',"`s2'",0)
			scalar `lnf' = r(ll)
		}
	}
end