具有给定最大值的离散概率分布

Question

我必须生成具有给定最大值的n元素的随机离散概率分布。

我知道如何创建经典的但我没有线索如何生成尊重n的{​​{1}}元素的随机离散概率分布。我的意思是分布的最高概率必须是max(distribution) = p。

有什么想法吗？

Answer 1

您可以使用命中和未命中方法（通常用于概率模拟，以随机选择满足特定约束的元素）。除非np太接近1 ，否则性能是可接受的。这是一个Python实现，应该很容易翻译成其他语言：

from random import random

def makeDist(n,p):
    #assumes p*n > 1
    miss = True
    while miss:
        x = [p]
        rest = [random() for i in range(1,n)]
        total = sum(rest)
        rest = [x*(1-p)/total for x in rest]
        if all(x < p for x in rest):
            x.extend(rest)
            miss = False #we have a hit!
    return x

典型输出：

>>> dist = makeDist(6,0.2)
>>> dist
[0.2, 0.08986510724051082, 0.18690143846768711, 0.19758176720598397, 0.19299989610231708, 0.13265179098350102]
>>> max(dist)
0.2
>>> sum(dist)
1.0
>>> makeDist(3,0.35)
[0.35, 0.31154704906869274, 0.33845295093130723]
>>>