0

我正在尝试从具有第 3 页上的 MATLAB 实现的论文中实现经验分布函数。这是我的 Python 版本。

我根据NumPy for MATLAB users 文档对其进行了转换，同时考虑了如何statsmodels 实现 ECDF

from statsmodels.distributions.empirical_distribution import ECDF     

def ecdf_representation(D, n):
     """calculate ECDF from D at n points"""
     m = np.mean(D)
     X = []
     for d in xrange(D.shape[0] + 1):
         func = ECDF([D[:, d] + np.random.randn(np.shape(D[:, d])) * 0.01])
         ll = func(np.linspace(0, 1, n))
         X = [X, ll]
     X = [X, m]
     plt.plot(X)
     plt.show()
     return X

我得到错误：

line 25, in ecdf_representation
func = ECDF([D[:, d] + np.random.randn(np.shape(D[:, d]))])
IndexError: too many indices for array

不D.shape[0]返回列数？那么，D[:, d]应该正常工作吗？这里发生了什么？

4

1 回答 1

0

D.shape[0] 将返回行，而不是列：

.shape[] 在“for i in range(Y.shape[0])”中做了什么？

D.shape[1] 将返回列

于 2017-06-30T13:34:15.590 回答