[BUG] - parameters alpha, beta, gamma returned in Fit_Gamma_3P are not the same of Gamma_Distribution? about reliability HOT 1 CLOSED

I'm not sure why you think the parameters are switched. The histogram plot with fitted PDF just peaks really high. If you let reliability plot it using results.best_distribution.PDF() then the ylim will be set automatically.

Just because the PDF of the Gamma 3P doesn't follow the histogram does not mean that the fit is incorrect. Mathematically is is better than all the others so it should be chosen, but it is only the best of the 12 distributions in this package. It may just be that none of the available distributions are suited to the data. "But Gamma_3P is the best fitting distribution" I hear you say. That is true from a log likelihood perspective, but that doesn't mean that there isn't some other distribution out there that's better than Gamma_3P.

Out of interest, look at how well the Weibull Mixture models the histogram:

The dataset is interesting because of the non-zero start point, the few large peaks in the middle and the small peak in the tail.
Not many software packages implement Gamma 3P so I understand benchmarking it against other software is difficult. If you want to verify the accuracy of the Gamma_3P fit, then try drawing samples from a known distribution and refitting it like is done here.

Lastly, before I close this issue for not being a bug, I would like to point out that there is a much better way to access the parameters than using iloc and rebuilding the distribution. Do this instead:

results = Fit_Everything(data, show_histogram_plot=False,show_probability_plot=False,show_best_distribution_probability_plot=False,show_PP_plot=False)
results.best_distribution.PDF()
print("Alpha =",results.best_distribution.alpha)

The best fitting distribution is already built for you and you can extract parameters by name. Have a read of the API reference for Fit_Everything.

If you'd like to know more, please feel free to email [email protected].

from reliability.