2012-09-10 56 views
4

我是Mathematica的初學者,只能理解到目前爲止的基本命令。在編寫最大似然估計函數時需要幫助。我想用2個參數找到二項分佈的MLE:n和p。我知道這是一個非常簡單的問題,實際上我可以手動獲得答案,這是我做的。但是我需要知道如何在Mathematica中編寫這個以獲得相同的答案。謝謝!如何在Mathematica中編寫最大似然估計?

回答

6
data = RandomVariate[BinomialDistribution[2, .3], 1000]; 
EstimatedDistribution[data, BinomialDistribution[n, p], 
         ParameterEstimator -> {"MaximumLikelihood"}] 

(* -> BinomialDistribution[2, 0.3155] *) 
+0

非常感謝!我用這個代碼,它的工作原理。然而,我得到這個消息爲我的實際數據:FindRoot :: cvmit:無法在100次迭代內收斂到要求的精度或精度。二項式分佈[4,0.0252702]。你知道那裏發生了什麼嗎? – user1525225

+1

@ user1525225你可以嘗試這裏描述的其他方法http://reference.wolfram.com/mathematica/ref/ParameterEstimator.html。表格示例'ParameterEstimator - > {「MaximumLikelihood」,Method - >「NMaximize」}]' –