如何使用“ MME”方法拟合威布尔分布并在R中找到估计值

Question

我正在尝试使用MME方法拟合R中的Weibull分布并找到估计值。下面是来自fitdistrplus和actuar软件包的R中的weibull分布拟合代码。

value

但是我收到以下错误

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

如果有人可以告诉我我犯了什么错误或提供相同的阅读材料，这将很有帮助

Answer 1

下面的R代码，Microsoft R Open 3.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方法产生了合理的拟合度