【问题标题】:Scaling a PDF to show 100% at peak缩放 PDF 以在峰值显示 100%
【发布时间】:2019-05-21 03:22:54
【问题描述】:

我正在显示我的数据的直方图,并带有一个重叠的 PDF。我的情节看起来都是这样的:

我正在尝试缩放红色曲线以在峰值处显示 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')

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

这是my previous question 的后续行动。

【问题讨论】:

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


    【解决方案1】:

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

    【讨论】:

      【解决方案2】:

      你使用.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

      【讨论】:

        猜你喜欢
        • 1970-01-01
        • 1970-01-01
        • 2015-12-27
        • 1970-01-01
        • 1970-01-01
        • 1970-01-01
        • 1970-01-01
        • 1970-01-01
        • 2021-03-20
        相关资源
        最近更新 更多