Я пытаюсь закодировать это гамма-распределение на Python; однако мне трудно понять, как включить интеграл в функцию.
Распределение: gamma.cdf
Будет ли первая половина уравнения прежде чем интеграл будет определен как отдельная функция, а затем умножен на функцию, определенную интегралом?
Для большего контекста: «t» — это моя переменная времени, которая состоит из моих точки x-данных в диапазоне от 0 до 24.
Это моя попытка написать это:
Код: Выделить всё
#Define Gamma Function
def func(t,g,n):
y = (1/gamma(g)*n*g)
return y
def func1(t,g,s,n,d):
y1 = (s**(g-1))*(np.exp(1)**(-s/n))*d*s
return y1
I = integrate.quad(func1,0,t)
funcdata = func(xdata,2,1)*func1(xdata,1,1,1,1)
plt.plot(xdata,funcdata,label="Model")
Источник: https://stackoverflow.com/questions/689 ... -in-python
