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 。