4

问题

我计算了一个取决于两个变量的概率密度函数。我想使用这个多元分布来生成一些随机数，这些随机数的概率与 PDF 成正比。

看起来，SciPy 目前只支持单变量分布。是否有任何允许二维分布的简单方法或易于使用的软件包？

作为一种解决方法，我可能会尝试在感兴趣的域上创建随机数并将它们丢弃或保留它们与我的 PDF 相关的机会，但仍然可能有其他选择。随机数生成不必很快。

谢谢你的帮助！

这是一个可能的解决方案

根据答案（非常感谢！），我破解了一些代码，您可能会在此 gist中找到这些代码。如果您使用 sin^2*Gauss PDF 运行此示例，将在 PDF 上绘制 2000 个满足给定条件（在圆圈内）的随机随机变量。也许这对其他人也有帮助。

4

2 回答 2

2

所以你有一个PDF F(x,y)，你想生成成对x并y根据这个PDF分发？

我会说，除非您可以使用反演技术（wiki）的多元版本，否则拒绝采样是可行的方法。

于 2012-06-18T14:14:24.933 回答
0

对于变量 X 和 Y，您不能通过仅生成具有 X 的独立分布的 x 和具有给定 x 的 Y 分布的 ay 来将其分成两个单变量分布的采样吗？

于 2012-06-18T14:14:28.160 回答