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

R 如何确定最大似然法的初始点,r,mle,R,Mle,我目前正在做配电装置。我使用了fitdistr函数,但在确定MLE的初始点时遇到了问题。例如,我想用伽马分布拟合我的数据(降雨量-13149乘以1矩阵) fit.gamma = fitdistr(rainfall,dgamma,start=list(shape = ?, scale = ?),method="Nelder-Mead") 库FitDistripPlus对此非常有用。如果没有起始值,它将为您猜测gamma参数。此外,如果你的猜测失败,你可以使用矩量法 x <- rgamm

我目前正在做配电装置。我使用了
fitdistr
函数,但在确定MLE的初始点时遇到了问题。例如,我想用伽马分布拟合我的数据(降雨量-13149乘以1矩阵)

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

FitDistripPlus
对此非常有用。如果没有起始值,它将为您猜测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

x也许这会有帮助:谢谢@Pascal的建议谢谢你。。我找到了几个用于此目的的软件包,包括FitDistripPlus,但不知道哪一个是最好的。你的意思是使用矩量法来寻找起始值吗?@lagizelt我必须对起始值进行一些研究,我不知道是否有一些方法可以找到最佳的起始点,从而给出最佳的估计结果。好的,非常感谢@lagizelt的帮助回答