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

Question

我目前正在研究配电装置。我使用fitdistr了函数，但在确定 MLE 的初始点时遇到了问题。例如，我想用伽玛分布拟合我的数据（雨量 - 13149 乘 1 矩阵）。

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

Answer 1

图书馆fitdistrplus对此非常有帮助。如果您没有起始值，它会为您猜测 gamma 参数。此外，如果您的猜测失败，您可以使用矩量法。

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