期望的概率分布

Question

我有一个数据＆＃34; zcosmo＆＃34; 遵循某种分布。它有很多波动（即它在某个点稍微增加然后变平等，但总体上从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）=总和（从0到y）p（y）dy

• 3）抛出一个随机数统一偏差x在0和1之间

• 4）计算y = F-1（x）

＆＃34; y＆＃34; 将具有所需的概率分布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堆栈。它包含probability and statistics，numerical integration和plotting的优秀图书馆。