подгонка обратного распределения Гаусса к данным в R

Question

Я пытаюсь использовать функцию fitdist в R для подгонки данных к трем различным распределениям по максимальной вероятности их сравнения.Логнормал и Вейбулл работают нормально, но я борюсь с обратным гауссовским.Мне нужно указать начальные значения, однако, когда я получаю сообщение об ошибке.

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

Я работал с подобной проблемой и обнаружил, что проблема заключалась в том, как я пометил мои начальные значения.Библиотека актуария, с которой я работал, требовала меток «среднее» и «форма» для значений.Следующий код предоставил мне решение:

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