如何确定最大似然法的初始点
How to determine the initial points of the maximum likelihood method
我目前正在研究分布拟合。我使用了 fitdistr 函数，但在确定 MLE 的初始点时遇到问题。例如，我想用伽马分布拟合我的数据（降雨量 - 13149 x 1 矩阵）。
fit.gamma = fitdistr(rainfall,dgamma,start=list(shape = ?, scale = ?),method="Nelder-Mead")
库 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
我目前正在研究分布拟合。我使用了 fitdistr 函数，但在确定 MLE 的初始点时遇到问题。例如，我想用伽马分布拟合我的数据（降雨量 - 13149 x 1 矩阵）。
fit.gamma = fitdistr(rainfall,dgamma,start=list(shape = ?, scale = ?),method="Nelder-Mead")
库 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