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的帮助回答