Случайное гамма-распределение с использованием python

Question

Я пытаюсь подогнать к своим данным случайное гамма-распределение, используя python. До сих пор я написал следующие коды

a = 10
b = 1
n=25
n_classes = 16       
true_val = []
for i in np.arange(n_classes):
    gamma_process= np.cumsum(np.random.gamma(n,a,b))   # to get cummaltive sum of the fit
    true_val = gamma_process   

Когда я запускал этот код, он не выдает никаких ошибок, но значение, которое я получаю в true_val, только одно. Полагаю получить список ценностей. Я не понимаю, где я делаю ошибку. Буду признателен, если вы мне поможете в этом. Заранее спасибо

Answer 1

Попробуйте добавить в список результатов:

true_val = true_val + [gamma_process[0]]

или лучше:

 true_val.append(float(gamma_process))

Теперь результирующий список имеет целочисленный тип.

Answer 2

Это связано с размером гамма-распределения. Пожалуйста, проверьте, если b> 2.

In numpy гамма-функция,

numpy.random.gamma(shape, scale=1.0, size=None)

Пожалуйста, проверьте еще раз.