0
对于Python的numpy或scipy,特别是Python 3,是否支持将概率表示为各种分布的对数概率?如果是这样,在哪里?有没有对log概率的支持?
当概率非常小时,这些特别重要。那么,例如,在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)
有什么阻止你只是转换到非对数表示?即使下溢并且您的概率变为0,也不会影响样本。 – user2357112