將逆高斯分布擬合到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))