python - 期望的概率分布

Question

我有一个遵循特定分布的数据“zcosmo” 。它有很多波动（即它在某个点略微增加然后变平等，但总体而言它从 0 增加到 0.5）。所以我使用一个遵循数据分布的函数来拟合数据。

This function is the probability distribution lets say of the data. 在这张图片中，拟合给了我一条蓝线，我想将其用作概率分布函数。

def fit(x,p1,p2):
    return (p1*x)+(p2*(x**2))

现在我理论上知道如何创建遵循这个函数的数据了！

• 1）根据簇红移的分布（拟合函数）计算红移中的选择函数并将其归一化 以具有概率分布，例如p（y）

• 2) 计算概率分布的不定积分 F(y)= sum (from 0 to y) p(y)dy

• 3) 抛出一个在 0 和 1 之间均匀偏离 x 的随机数

• 4) 计算 y = F-1(x)

然后“y”将具有所需的概率分布 p(y)

I am struggling to find a way to do it in Python. Is there a package that does the above??

Answer 1

你应该看看SciPy堆栈。它包括用于概率和统计、数值积分和绘图的优秀库。