感觉stata定义likelihood function的方式怪怪的,没有Gauss自然. Stata的语法也很
古怪,比方说下面程序里的美元符号还有奇怪的引号等。我不知道新版本有无改进。
capture program drop mle_poisson
program define mle_poisson
version 8.0
args lnf theta1
quietly replace `lnf'=-exp(`theta1')+$ML_y1*(`theta1')-lngamma($ML_y1+1)
end
ml model lf mle_poisson (y=x)
ml maximize