numpy.random.logistic#
- random.logistic(loc=0.0, scale=1.0, size=None)#
从逻辑斯蒂分布中抽取样本。
从具有指定参数的逻辑斯蒂分布中抽取样本,loc(位置或均值,也是中位数)和 scale(>0)。
- 参数:
- locfloat 或浮点数的类数组,可选
分布的参数。默认为 0。
- scalefloat 或浮点数的类数组,可选
分布的参数。必须是非负数。默认为 1。
- sizeint 或 int 的元组,可选
输出形状。如果给定的形状为,例如,
(m, n, k)
,则会抽取m * n * k
个样本。如果 size 为None
(默认),如果loc
和scale
都是标量,则返回单个值。否则,将抽取np.broadcast(loc, scale).size
个样本。
- 返回:
- outndarray 或标量
从参数化的逻辑斯蒂分布中抽取的样本。
参见
scipy.stats.logistic
概率密度函数、分布或累积密度函数等。
random.Generator.logistic
应该用于新代码。
注释
逻辑斯蒂分布的概率密度为
\[P(x) = P(x) = \frac{e^{-(x-\mu)/s}}{s(1+e^{-(x-\mu)/s})^2},\]其中 \(\mu\) = 位置,\(s\) = 尺度。
逻辑斯蒂分布用于极值问题,在其中它可以充当 Gumbel 分布的混合,在流行病学中,以及国际象棋联合会 (FIDE) 使用它,它在 Elo 排名系统中使用,假设每个棋手的表现是逻辑斯蒂分布的随机变量。
参考文献
[1]Reiss, R.-D. 和 Thomas M. (2001),“极值的统计分析,来自保险、金融、水文和其他领域”,Birkhauser Verlag,巴塞尔,第 132-133 页。
[2]Weisstein,Eric W.“逻辑斯蒂分布。”来自 MathWorld——Wolfram 网页资源。https://mathworld.wolfram.com/LogisticDistribution.html
[3]维基百科,“逻辑斯蒂分布”,https://en.wikipedia.org/wiki/Logistic_distribution
示例
从分布中抽取样本
>>> loc, scale = 10, 1 >>> s = np.random.logistic(loc, scale, 10000) >>> import matplotlib.pyplot as plt >>> count, bins, ignored = plt.hist(s, bins=50)
# 对比分布作图
>>> def logist(x, loc, scale): ... return np.exp((loc-x)/scale)/(scale*(1+np.exp((loc-x)/scale))**2) >>> lgst_val = logist(bins, loc, scale) >>> plt.plot(bins, lgst_val * count.max() / lgst_val.max()) >>> plt.show()