Ctrl+K
scipy.stats.Mixture.

ilogcdf#

Mixture.ilogcdf(p, /, *, method=None)[源代码]#

累积分布函数对数的逆函数。

累积分布函数对数的逆函数（“逆对数 CDF”）是参数 \(x\)，该参数使得累积分布函数的对数 \(\log(F(x))\) 的值为 \(\log(p)\)。

在数学上，它等价于 \(F^{-1}(\exp(y))\)，其中 \(y = \log(p)\)，但是与朴素的实现（计算 \(p = \exp(y)\)，然后 \(F^{-1}(p)\)）相比，它在数值上可能更优。

ilogcdf 接受 logp，其条件为 \(\log(p) ≤ 0\)。

参数:
logparray_like

逆对数 CDF 的参数。

method{None, ‘formula’, ‘complement’, ‘inversion’}

用于评估逆对数 CDF 的策略。默认值 (None) 时，基础设施会在以下选项中选择，按优先级顺序排列。

  • 'formula'：使用逆对数 CDF 本身的公式

  • 'complement'：在 logp 的对数补码处评估逆对数 CCDF（参见注释）

  • 'inversion'：数值求解对数 CDF 等于 logp 的参数

并非所有分布都支持所有 method 选项。如果所选的 method 不可用，则会引发 NotImplementedError。

返回:
outarray

在提供的参数处评估的逆对数 CDF。

另请参见

icdf
logcdf

注释

假设连续概率分布的支持度为 \([l, r]\)。逆对数 CDF 在 \(\log(p) = \log(0) = -\infty\) 时返回其最小值 \(l\)，在 \(\log(p) = \log(1) = 0\) 时返回其最大值 \(r\)。由于对数 CDF 的范围是 \([-\infty, 0]\)，因此逆对数 CDF 仅在负实数上定义；对于 \(\log(p) > 0\)，ilogcdf 返回 nan。

有时，需要找到 CDF 的参数，使得结果概率非常接近 0 或 1 - 太接近而无法用浮点运算精确表示。然而，在许多情况下，这个结果概率的对数可以用浮点运算表示，在这种情况下，可以使用此函数来找到 CDF 的参数，使得结果概率的对数为 \(y = \log(p)\)。

数字 \(z\) 的“对数补码”在数学上等效于 \(\log(1-\exp(z))\)，但是当 \(\exp(z)\) 接近 \(0\) 或 \(1\) 时，会进行计算以避免精度损失。

示例

使用所需的参数实例化分布

>>> import numpy as np
>>> from scipy import stats
>>> X = stats.Uniform(a=-0.5, b=0.5)

在所需的参数处评估逆对数 CDF

>>> X.ilogcdf(-0.25)
0.2788007830714034
>>> np.allclose(X.ilogcdf(-0.25), X.icdf(np.exp(-0.25)))
True
在本页