from scipy.stats import probplot for i in X.columns: probplot(x=X[i],dist='norm',plot=plt) plt.title(i) plt.show()