我有一个关于概率分布函数的问题我有一个时间序列数据，我想计算不同时间窗口内数据的概率分布。
我开发了以下代码，但我找不到该函数的概率分布值。
a = pd.DataFrame([0.0,
21.660332407421638,
20.56428943581567,
20.597329924045983,
19.313207915827956,
19.104973174542806,
18.031361568112377,
17.904747973652125,
16.705687654209264,
16.534206966165637,
16.347782724271802,
13.994284547628721,
12.870120434556945,
12.794530081249571,
10.660675400742669])
这是我的数据的直方图和密度图：
a.plot.hist()
a.plot.density()
但我不知道如何计算密度曲线下的面积值。
可以直接调用该方法scipy.stats.gaussian_kdepandas 内部也使用它。
该方法返回所需的函数。
然后您可以调用其中的方法之一scipy.integrate计算核密度估计下的面积，例如
from scipy import stats, integrate
kde = stats.gaussian_kde(a[0])
# Calculate the integral of the kde between 10 and 20:
xmin, xmax = 10, 20
integral, err = integrate.quad(kde, xmin, xmax)
x = np.linspace(-5,20,100)
x_integral = np.linspace(xmin, xmax, 100)
plt.plot(x, kde(x), label="KDE")
plt.fill_between(x_integral, 0, kde(x_integral),
alpha=0.3, color='b', label="Area: {:.3f}".format(integral))
plt.legend()
本文内容由网友自发贡献，版权归原作者所有，本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容，请联系:hwhale#tublm.com(使用前将#替换为@)
