Масштабирование PDF, чтобы показать 100% на пике - PullRequest
0 голосов
/ 21 мая 2019

Я показываю гистограмму моих данных с наложенным 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 с этим?Или есть лучший способ сделать это?

Это своего рода продолжение моего предыдущего вопроса .

Ответы [ 2 ]

0 голосов
/ 21 мая 2019

Вы можете использовать density=True при построении гистограммы.

0 голосов
/ 21 мая 2019

Вы используете .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

...