缩放PDF以显示100%的峰值

时间:2019-05-21 03:22:54

标签: python-3.x matplotlib scipy probability-density

我正在显示数据的直方图和覆盖的PDF。我的情节看起来都像这样:

enter image description here

我正在尝试缩放红色曲线以在峰值处显示100%。

以下两个玩具代码与我实际使用的玩具代码相同,除了两个%之间的行:

%
import pandas as pd
import matplotlib.pyplot as plt
import scipy.stats as stats
import numpy as np
my_randoms = np.random.normal(0.5, 1, 50000)
dictOne = {"delta z":my_randoms}
df = pd.DataFrame(dictOne)
df = df[df['delta z'] > -999]
%

fig, ax = plt.subplots()
h, edges, _ = ax.hist(df['delta z'], alpha = 1, density = False, bins = 100)

param = stats.norm.fit(df['delta z'].dropna())   # Fit a normal distribution to the data
pdf_fitted = stats.norm.pdf(df['delta z'], *param)
x = np.linspace(*df['delta z'].agg([min, max]), 100) # x-values
binwidth = np.diff(edges).mean()
ax.plot(x, stats.norm.pdf(x, *param)*h.sum()*binwidth, color = 'r')

#   Decorations
graph_title = 'U-B'
plt.grid(which = 'both')
plt.title(r'$\Delta z$ distribution for %s'%graph_title, fontsize = 25)
plt.xlabel(r'$\Delta z = z_{spec} - z_{photo}$', fontsize = 25)
plt.ylabel('Number', fontsize = 25)
plt.xticks(fontsize = 25)
plt.yticks(fontsize = 25)
xmin, xmax = min(df['delta z']), max(df['delta z'])
plt.xlim(xmin, xmax)
plt.annotate(
r'''$\mu_{\Delta z}$ = %.3f
$\sigma_{\Delta z}$ = %.3f'''%(param[0], param[1]),
         fontsize = 25, color = 'r', xy=(0.85, 0.85), xycoords='axes fraction')

我该如何在右侧定义另一个axes对象(从0到100)并将PDF映射到该对象? 还是有更好的方法呢?

这是对my previous question的后续行动。

2 个答案:

答案 0 :(得分:0)

您使用.twinx()

fig = plt.figure(figsize=(10, 8), dpi=72.0)

n_rows = 2
n_cols = 2

ax1 = fig.add_subplot(n_rows, n_cols, 1)
ax2 = fig.add_subplot(n_rows, n_cols, 2)

ax3 = ax1.twinx()

https://matplotlib.org/gallery/api/two_scales.html

答案 1 :(得分:0)

您可以使用density=True绘制直方图。