#qq plot using numpy and matplotlib p = np.linspace(0,100,100) a = np.percentile(x,p) b = np.percentile(y,p) plt.plot(a,b, ls="", marker="o") x = np.linspace(np.min((a.min(),b.min())), np.max((a.max(),b.max()))) plt.plot(x,x, color="r", ls="--")