I am using python to calculate Inverse of a CDF of gamma distribution (using scipy. stat.gamma.fit). But for probability value 1, it is coming infinite. If it is replaced from 1 to 0.99 it works but the values changes with the different number of significant figures. Like it is 61 for 0.99 and 130 for 0.9999.
I do not know the best way to handle these infinite values in my workflow.I need to get a ideals but valid inverse value for probabilities arbitrarily close to 1. But i don't know, how to decide how many digits are needed to round off the probability value. Also, I am not sure about, Will it be fine to round off?
Can we do something to get some meaningful information instead of inf value?
Thank you for your help