HistCumSum[-1]
- окончательное значение в массиве HistCumSum
.Поскольку HistCumSum
является накопленной суммой, мы можем видеть Result = HistCumSum / HistCumSum[-1]
как , нормализующий от совокупной частоты к совокупной функции распределения .
Взгляд на то, как выглядят Hist
, HistCumSum
и Result
, может помочь разобраться:
import numpy as np
from scipy.stats import norm
import matplotlib.pyplot as plt
x = np.linspace(-2, 2, num=100)
# Evaluate pdf in lieu of actual data
Hist = norm.pdf(x)
# Compute cumulative sum
HistCumSum = np.cumsum(Hist)
# Normalise cumulative sum to produce cdf
Result = HistCumSum / HistCumSum[-1]
fig, ax = plt.subplots(1, 3, figsize=[7, 3])
ax[0].plot(x, Hist)
ax[1].plot(x, HistCumSum)
ax[2].plot(x, Result)
ax[0].set_title('Hist')
ax[1].set_title('HistCumSum')
ax[2].set_title('Results (CDF)')
fig.tight_layout()