如何确定最大似然法的初始点

Question

我目前正致力于分销配件。我使用了fitdistr函数，但在确定MLE的初始点时遇到了问题。例如，我想用γ分布拟合我的数据（降雨量为13149乘1矩阵）。

fit.gamma = fitdistr(rainfall,dgamma,start=list(shape = ?, scale = ?),method="Nelder-Mead")  

Answer 1

图书馆fitdistrplus非常适合这个。如果您没有起始值，它会为您猜测伽玛参数。此外，如果您的猜测失败，您可以使用时刻方法。

x <- rgamma(100, 0.5, 0.5)

library(fitdistrplus)
(pars <- fitdist(x, "gamma"))
# Fitting of the distribution ' gamma ' by maximum likelihood 
# Parameters:
#        estimate Std. Error
# shape 0.4443304 0.05131369
# rate  0.5622472 0.10644511