Мне нужно создать усеченную кривую pdf гамма-распределения и гистограмму в Python 3.2 на win7.
import numpy as np
import matplotlib.pyplot as plt
import scipy.special as sps
shape, scale = 2., 2. # mean and dispersion
counter =0
s = []
upper_bound = 4
lower_bound = 0.5
while (counter <= 1000):
t = np.random.gamma(shape, scale, 1)
if (lower_bound <= t <= upper_bound) :
s.append(t)
counter += 1
count, bins, ignored = plt.hist(s, 50, normed=True)
// this part take s very long time
y = bins**(shape-1)*(np.exp(-bins/scale) /
(sps.gamma(shape)*scale**shape))
plt.plot(bins, y, linewidth=2, color='r')
plt.show()
Я обнаружил, что следующий код занимает очень много времени, и всплывающее окно рисунка перестает отвечать на запросы.
"y = bins**(shape-1)*(np.exp(-bins/scale)/(sps.gamma(shape)*scale**shape))"
Если я уберу нижнюю и верхнюю границы гамма-распределения, оно будет работать очень быстро.
Любая помощь будет оценена по достоинству.