разность двух случайных распределений питона

Question

Добрый день! У меня есть два гамма-распределения, и я хочу найти распределение их разницы. Используйте np.random.gamma для генерации распределения по параметрам, но результирующее распределение время от времени сильно отличается. Код:

import numpy as np 
from scipy.stats import gamma

for i in range(0, 10):
    s1 = np.random.gamma(1.242619972, 0.062172619, 2000) +  0.479719122 
    s2 = np.random.gamma(456.1387112, 0.002811328, 2000) - 0.586076723
    r_a, r_loc, r_scale = gamma.fit(s1 - s2)
    print(1 - gamma.cdf(0.0, r_a, r_loc, r_scale))

Результат:

0.4795655021157602
0.07061938039031612
0.06960741675590854
0.4957568913729331
0.4889900326940878
0.07381963810128422
0.0690800784280835
0.07198551429809896
0.07659274505827551
0.06967441935502583

Я получаю два совершенно разных файла формата cdf: 0,48 и 0,07. В чем может быть проблема?

Answer 1

Вы подгоняете гамма-распределение к разнице между двумя другими гамма-распределениями. Распределение гаммы может быть только положительным, так что это не имеет смысла, и вы не можете ожидать последовательного ответа. Если вы напечатаете среднюю разницу, вы получите последовательные результаты.