How to write a maximum likelihood by yourself

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

各位好，我现在需要写一个maximum likelihood function,这是一个underreported 的poisson model,所以没有现成的程序。现在的问题是在这个ml中，遇到了一个积分，所以需要用gauss-hermite数值积分的办法解决，请问各位知道怎么实现这个数值积分吗？
likelihood function 是: 对于每个individual的prob function is
g_i=Integral{exp(-exp(X'*beta+sqrt(2)*sigma*v))*Phi((Z'*gamma+sqrt(2)*rho*v)/(sqrt(1-rho^2)))*[exp(X'*beta+sqrt(2)*sigma*v)*Phi((Z'*gamma+sqrt(2)*rho*v)/(sqrt(1-rho^2)))]^(y)exp(-v^2)dv}
我现在有一个小程序program my_truncatedpoisson_lf
version 9.1
args lnf beta gama rho sigma

tempvar a b c d v w
gen `v'=.
gen `w'=.
_GetQuad, avar(`v') wvar(`w') quad(20)

gen double `a'=(`gama'+1/(exp(`rho')+1)*`v'/exp(`sigma'))/sqrt(1-1/((exp(`rho')+1)^2))
gen double `b'=exp(`beta'+`v')*normal(`a')
gen double `c'=1/exp(lnfactorial($ML_y1))*exp(-`sigma')*normalden(`v',0,exp(`sigma'))
gen double `d'=`w'*exp(-`b')*`b'^($ML_y1)*`c'

collapse (sum)`d'
qui replace `lnf'=log(`d')
end

use auto
ml model lf my_truncatedpoisson_lf (beta: mpg= price)(gama: turn)()()
ml maximize
多谢各位了

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：Likelihood Yourself Maximum write Your

同问，也曾遇到这样的问题

看看这个命令 quadchk 的代码，是否有gauss-hermite

3# eblog
万分感谢