Skewnorm不合适

Question

这是我先前的问题here的后续内容。我正在尝试使用scipy.stats.skewnorm使this csv文件中的数据适合我，但是我无法正常工作：

import matplotlib.pyplot as plt
import pandas as pd
from scipy.stats import skewnorm

df = pd.read_csv('astro_data.csv')
x = df['delta z']
number_bins = 50

fig, ax = plt.subplots()
h, edges, _ = ax.hist(x, alpha = 0.5,
                      density = False,
                      bins = number_bins)

a_est, loc_est, scale_est = skewnorm.fit(x)
ax.plot(x, skewnorm.pdf(x, a_est, loc_est, scale_est), 'r-', lw=5, alpha=0.6, label='skewnorm pdf')

谁能看到我该如何解决这个问题？

编辑：当我更改为density=True时，结果是这样的：