Инверсия scipy.stats.lognorm.interval

Question

from scipy.stats import lognorm
posterior_fb = lognorm(s=np.log(1.14), scale=0.007)
intervals = post_fb.interval(0.99)

post_fb.interval(0.99) дает мне конечные точки диапазона, который содержит 99% распределения, т. Е. (0,0049, 0,0098)

Мне нужна функция, которая выполняет обратное.То есть вы можете указать две точки, и функция вычислит процент распределения, которое лежит в этих двух точках.

Например, inverse_of_interval(0.0049, 0.0098) даст мне 0,99.

Answer 1

Вы можете попробовать использовать cdf логнорма:

 post_fb = lognorm(s=np.log(1.14), scale=0.007)
 post_fb.cdf(0.0098) - post_fb.cdf(0.0049)