如何用“MME”方法拟合Weibull分布并求R中的估计

Question

我试图用MME方法拟合R中的Weibull分布，并找到估计值。下面是用于拟合R中的weibull分布的代码，从适配包到实数包。

a <- rweibull(100, 10,1)
weibul_mme <- mmedist(a, "weibull", order = 1:2)

但是我得到了下面的错误

Error in mmedist(a, "weibull", order = 1:2) : 
  the empirical moment function must be defined

如果有人能告诉我我犯了什么错误，或者提供一些同样的阅读材料，那将是有帮助的。

Answer 1

下面的r代码，MicrosoftROpen3.5.3 x64，Win 10

library(fitdistrplus)
library(actuar)

x <- rweibull(1000, 10,1)

memp <- function(x, order) mean(x^order)

weibul_mme <- mmedist(x, "weibull", order = 1:2, memp=memp)
fit.weibull<- fitdist(x, "weibull", method = "mme", order=c(1, 2), memp=memp, lower = c(0, 0))

print(weibul_mme)
print(fit.weibull)

plot(fit.weibull, demp=TRUE)

用MME法产生合理的拟合

​

​