2015-10-10 45 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

回答