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()
Квантиль-квантильный график с использованием SciPy