将逆高斯分布拟合到R中的数据

Question

我试图在R中使用fitdist函数以最大可能性将数据拟合为三个不同的分布，以将它们进行比较。对数正态和Weibull可以很好地工作，但是我在逆高斯方面挣扎。 我需要指定起始值，但是当我收到错误消息时。

fw<-fitdist(claims,"weibull") WORKS
fln<-fitdist(claims,"lnorm") WORKS
fig<-fitdist(claims,"invgauss",start=list(mu=0,lambda=1)) DOES NOT WORK

错误：“当输入的长度为零且未引发错误时，pinvgauss函数应返回零长度的矢量”

我的代码有什么问题？

Answer 1

我正在处理一个类似的问题，发现问题在于我如何标记我的起始值。我正在使用的执行器库在值上需要标签“ mean”和“ shape”。以下代码为我提供了解决方案：

library(actuar)
library(fitdistrplus)
fig <- fitdist(claims, "invgauss", start = list(mean = 5, shape = 1))