scipy.stats.Mixture.
mode#
- Mixture.mode(*, method=None)[源代码]#
众数(最可能的值)
非正式地,众数是随机变量具有最高概率(密度)的假设值。也就是说,众数是支持\(\chi\)的元素,它最大化概率密度(或质量,对于离散随机变量)函数\(f(x)\)
\[\text{mode} = \arg\max_{x \in \chi} f(x)\]- 参数:
- method{None, ‘formula’, ‘optimization’}
用于评估众数的策略。默认情况下 (
None
),基础设施在以下选项之间进行选择,并按优先级顺序排列。'formula'
: 使用公式计算中位数'optimization'
: 数值上最大化 PDF/PMF
并非所有method选项都适用于所有分布。如果所选的method不可用,则会引发
NotImplementedError
。
- 返回:
- outarray
众数
注释
对于某些分布
众数不是唯一的(例如,均匀分布);
PDF 具有一个或多个奇点,并且奇点是否被认为在域中并称为众数是有争议的(例如,形状参数小于 1 的伽马分布);和/或
概率密度函数可能具有一个或多个不是全局最大值的局部最大值(例如,混合分布)。
在这些情况下,
mode
将返回单个值,
考虑在奇点处出现众数,和/或
返回局部最大值,该值可能是也可能不是全局最大值。
如果没有为所选分布专门实现众数公式,SciPy 将尝试以数值方式计算众数,这可能不符合用户首选的众数定义。在这种情况下,鼓励用户子类化分布并覆盖
mode
。参考
[1]Mode (statistics), Wikipedia, https://en.wikipedia.org/wiki/Mode_(statistics)
示例
使用所需的参数实例化分布
>>> from scipy import stats >>> X = stats.Normal(mu=1., sigma=2.)
评估众数
>>> X.mode() 1.0
如果众数不是唯一定义的,
mode
仍然返回单个值。>>> X = stats.Uniform(a=0., b=1.) >>> X.mode() 0.5
如果此选择不满足您的要求,请子类化分布并覆盖
mode
>>> class BetterUniform(stats.Uniform): ... def mode(self): ... return self.b >>> X = BetterUniform(a=0., b=1.) >>> X.mode() 1.0