是否支持对数概率？

Question

在 Python 的 numpy 或 scipy 中，特别是 Python 3，是否支持将概率表示为各种分布的对数概率？如果有，在哪里？

当概率非常小时，这些尤其重要。那么，例如，在numpy.random.choice 中，可以指定一个概率列表吗？对于我正在寻找的情况，最好将这些（全部）指定为日志概率。

例如，R 中的可比较能力是参数log.p，它在许多统计分布中出现或具有对应物。这意味着，在函数 qnorm 和 pnorm 的特定情况下，作为布尔开关，概率应该作为概率的日志返回，而不是概率本身。

例如，我想而不是说

s = numpy.random.binomial(100000, 0.0001, 1000)

能够说类似的话

s = np.random.binomial(100000, -9.2103, 1000, log_p=True)

Answer 1

Scipy stats 有一些可能适合您的需求：

scipy.stats.lognorm

• 日志正常:http://docs.scipy.org/doc/scipy-0.14.0/reference/generated/scipy.stats.lognorm.html
• 其他可以在那里找到：http://docs.scipy.org/doc/scipy/reference/stats.html