Python Scipy截斷指數分布

Question

我需要構造一個截斷的指數隨機變量，其范圍在5到7之間，比率參數等於0.76。 我正在使用scipy.stats.truncexpon，loc = 5和scale = 1 / 0.76。 我不確定如何指定上限。 有任何想法嗎？

Answer 1

import scipy.stats as stats
import matplotlib.pyplot as plt

lower, upper, scale = 5, 7, 1/0.76
X = stats.truncexpon(b=(upper-lower)/scale, loc=lower, scale=scale)
data = X.rvs(10000)

fig, ax = plt.subplots()
ax.hist(data, normed=True)
plt.show()

---

stats.truncexpon是rv_continuous子類的rv_continuous 。 rv_continuous類具有a和b參數，這些參數定義了分布支持的下限和上限 。 truncexpon a參數固定為0：

truncexpon = truncexpon_gen(a=0.0, name='truncexpon')

因此，默認情況下truncexpon的支持范圍從0到b 。 根據文檔 ，

truncexpon.pdf(x, b, loc, scale)等同於truncexpon.pdf(y, b) / scale其中y = (x - loc) / scale 。

因此，當y從0變為b ， x從loc變為b*scale + loc 。 因此，要使x從lower到upper ，讓loc = lower和b = (upper-lower)/scale 。