python - 制作自定义概率分布以从 SciPy 中抽取随机样本

Question

我正在寻找使用蒙特卡罗类型模拟来总结任意数量的事物概率分布。我想随机抽取某些东西的连续分布并将它们添加到其他连续分布的其他随机样本中，最终得到它们组合的概率分布。分布本身是经验性的 - 它们不是函数，而是 P99 = 2.4、P90 = 7.12、P50 = 24.53、P10 = 82.14 等形式（实际上有很多这些点）。这些分布或多或少是对数正态分布，因此如果有必要，将它们近似为对数正态分布可能会很好。但是我怎么能把它输入到 SciPy 的lognorm 函数中呢？或者在 SciPy 或一般的 python 中以其他方式做它？

我希望很清楚我想要做什么。非常感谢，亚历克斯

Answer 1

看起来你所拥有的本质上是概率密度的直方图。然后您可以做的一件事是将逆变换采样与您的经验分布一起使用。

作为替代方案，如果您期望某个分布的函数形式（对数范数或其他某种形式），您可以尝试使用相应的函数形式拟合数据。