python - 使用 SciPy 的分位数-分位数图

Question

你会如何使用 Python 创建一个 qq-plot？

假设您有大量测量值，并且正在使用一些将 XY 值作为输入的绘图函数。该函数应根据某些分布(正态、均匀...)的相应分位数绘制测量值的分位数。

结果图让我们可以评估我们的测量是否遵循假设的分布。

http://en.wikipedia.org/wiki/Quantile-quantile_plot

R 和 Matlab 都为此提供了现成的函数，但我想知道在 Python 中实现的最简洁的方法是什么。

Answer 1

更新:正如人们指出的那样，这个答案不正确。概率图不同于分位数-分位数图。在您在解释或传达您的分布关系时出错之前，请先查看这些评论和其他答案。

我认为 scipy.stats.probplot 会做你想做的事。见 documentation了解更多详情。

import numpy as np 
import pylab 
import scipy.stats as stats

measurements = np.random.normal(loc = 20, scale = 5, size=100)   
stats.probplot(measurements, dist="norm", plot=pylab)
pylab.show()

结果