fitdist 函数出错：函数 mle 无法估计参数，错误代码为 1（瑞利分布）答案

【问题标题】：Error in fitdist function: the function mle failed to estimate the parameters, with the error code 1 (Rayleigh distribution)fitdist 函数出错：函数 mle 无法估计参数，错误代码为 1（瑞利分布）
【发布时间】：2021-05-01 06:08:03
【问题描述】：

我正在尝试使用 fitdistrplus 包中的 fitdist 函数将我的数据拟合到瑞利分布。

x <- c(19.000000,23.000000,26.000000,45.000000,8.050000,46.900000,1.268333,30.000000,
       1.466667,3.733333,1.683333,4.000000,3.950000,1.850000,42.000000,1.333333,
       1.550000,1.000000,2.066667,1.566667,1.216667,1.850000,1.400000,8.366667,
       19.000000,29.000000,17.000000,42.000000,19.000000,10.000000,53.000000,2.550000,
       15.483333,1.533333,1.216667,1.550000,32.000000,6.583333,6.516667,5.750000,
       9.283333,46.000000,2.016667,2.133333,4.516667,46.950000,1.600000,1.433333,
       3.166667,4.416667,17.016667,2.433333,2.713333,8.633333,3.150000,1.183333,
       14.000000,10.706667,7.026944,31.000000,35.000000,21.000000,14.000000,2.200000,
       26.000000,3.316667,51.000000,13.000000,34.000000,11.650000,49.000000,12.000000,
       26.000000,20.000000,22.000000,6.483333,24.000000,5.333333,4.833333,8.750000,
       6.216667,17.000000,1.083333,19.000000,48.000000,15.000000,1.266667,54.000000,
       32.000000,3.616667,6.666667,1.600000,2.083333,6.933333,33.033333,1.883333,
       1.000000,3.072222,49.000000,1.400000)
dat <- data.frame(x)

# Generate gamma rvs
den <- density(x)
orig <- data.frame(x = den$x, y = den$y)

fit.params.2 <- fitdistrplus::fitdist(dat$x, "rayleigh", start = list(sigma = 1))

然后出现错误：

Error in fitdistrplus::fitdist(dat$x, "rayleigh", start = list(sigma = 1)) : 
   the function mle failed to estimate the parameters, 
                with the error code 1

这个问题有什么解决办法吗？感谢您的帮助。

【问题讨论】：

标签： r distribution fitdistrplus

【解决方案1】：

有两个问题。 (1) 瑞利分布似乎不太适合数据（请参见下面的绘图输出）和 (2) 您需要更好的起始值。由于 sigma 与瑞利分布的平均值成正比（请参阅wikipedia），请尝试：

library(fitdistrplus)
library(extraDistr)

fit <- fitdist(dat$x, "rayleigh", start = list(sigma = mean(dat$x)))
fit
## Fitting of the distribution ' rayleigh ' by maximum likelihood 
## Parameters:
##       estimate Std. Error
## sigma 15.00063   0.749935

plot(fit)

【讨论】：