dirichlet - numpy.random.dirichlet 有什么作用？

Question

我需要一个 Dirichlet 分布并且我正在使用 numpy.random.dirichlet。当我根据 Dirichlet PDF 公式给出 alpha=[1,1,1,1] 时，它应该产生一个统一的函数。但它没有给我一个统一的向量。有人知道为什么吗？

Answer 1

据我所知，numpy.random.dirichlet功能按预期工作。我想你可能有点误解狄利克雷分布是什么。

k 维狄利克雷分布将有 k 个参数(x_1 到 x_k)，它们的总和必须为 1。如果您尝试 sum(np.random.dirichlet([1, 1, 1, 1])) ，结果将始终为 1(给出或采取一些舍入误差)。所以虽然np.random.dirichlet([1, 1, 1, 1]))在其支持上确实是统一的，这种支持受到 x_1 + x_2 + x_3 + x_4 = 1 的要求的约束，并且输出数组不会是 Unif(0, 1) 分布的。

This article给出了一个相对简单和非技术性的概念解释。