我想将线叠加在此直方图上,使其两端都为零。我以为将0值添加到数据数组的开头和结尾就足够了,但这似乎行不通。
plt.style.use('seaborn-white')
data = np.array(sizes_)
n, bins, _ = plt.hist(data, bins=bins, alpha=0.2, density=True)
plt.grid(axis='y', alpha=0.5)
plt.xlim([0, 500])
bin_centers = 0.5 * (bins[1:] + bins[:-1])
## Tried this but doesn't work
np.insert(bin_centers, 0, 0)
np.append(bin_centers, 500)
np.insert(n, 0, 0)
np.append(n, 0)
##
plt.plot(bin_centers, n)
plt.xlabel('length (nts)')
plt.ylabel('frequency')
plt.title('Length Distribution')
plt.savefig(outfile)
有什么想法可以做到吗?
最佳答案
插入后需要保存数组。此外,您需要在数组末尾附加一对x-y值。由于您未提供MCVE,因此以下是使用示例数据集进行此操作的一种方法
plt.style.use('seaborn-white')
np.random.seed(121)
sizes_ = np.random.randint(2, 100, 50)
data = np.array(sizes_)
n, bins, _ = plt.hist(data, bins=20, alpha=0.2, density=True)
plt.grid(axis='y', alpha=0.5)
bin_centers = 0.5 * (bins[1:] + bins[:-1])
bin_diff = np.diff(bin_centers)[-1]
bin_centers = np.insert(bin_centers, 0, 0)
bin_centers = np.insert(bin_centers, len(bin_centers), bin_centers[-1] + bin_diff)
n = np.insert(n, 0, 0)
n = np.insert(n, len(n), 0)
plt.plot(bin_centers, n)
plt.xlabel('length (nts)')
plt.ylabel('frequency')
plt.title('Length Distribution')