Scipy：是否可以使用自定义概率密度函数制作qq图？

Question

我有一个2d pdf函数，它是一个双变量正态分布，

我绘制它angle-wisely，以便与我得到的原始数据进行比较。也就是说，我将(x,y)转换为(r, theta)，然后进行一些整合。使用Scipy.integrate

可以非常轻松地完成此操作

# Let f be the pdf
x_vals = np.arange(0, max_speed, 0.5)
y_vals_ =[integrate.nquad(f, 
                          [[x_val-0.5, x_val+0.5],[angle_radian-incre_radian/2, angle_radian+incre_radian/2]]) 
          for x_val in x_vals]

我可以将pdf与原始数据集进行比较

问题在于，density plot还不够，我还想制作一个q-q plot。但是，由于pdf是集成的结果，是否可以使用Scipy来实现？

我不知道。也许可能的解决方案是Sympy。由于bivariate pdf的参数都是已知的，我们可以生成角度方向的pdf。