2015-10-10 46 views
0

对于Python的numpy或scipy,特别是Python 3,是否支持将概率表示为各种分布的对数概率?如果是这样,在哪里?有没有对log概率的支持?

当概率非常小时,这些特别重要。那么,例如,在numpy.random.choice中,可以指定一个概率列表?对于我正在寻找的情况,将它们全部指定为日志概率将会很好。

例如,R中相当的能力是参数log.p,其在许多统计分布中出现或具有对应物。这意味着,在功能qnormpnorm的特定情况下,作为布尔开关,概率应该作为概率返回为,而不是概率本身。

我想,例如,而不是说

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

能够这样说

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

有什么阻止你只是转换到非对数表示?即使下溢并且您的概率变为0,也不会影响样本。 – user2357112

回答