np.log1p(df['population']).plot.hist(bins=30) plt.axvline(np.log1p(df['population']).mean(),c='C1') plt.axvline(np.log1p(df['population']).median(),c='C1',linestyle='--') plt.axvline(np.log1p(df['population']).quantile(0.25),c='C2',linestyle=':') plt.axvline(np.log1p(df['population']).quantile(0.75),c='C2',linestyle=':')