ECDFをPythonでどのように実装しますか？

Question

私はpaperから経験的な分布関数を実装しようとしていますが、これは3ページのMATLAB実装です。ここに私のPythonバージョンがあります。私は、アカウントにどのようにstatsmodelsECDF

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 

をしながら、それはNumPy for MATLAB users documentationに従って変換

私はエラーを取得する：

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]は正しく動作するはずですか？何が起きてる？

Answer 1

D.shape [0]列、行ではなく返されます。

What does .shape[] do in "for i in range(Y.shape[0])"?

D.shapeは、[1] OPでエラーを解決し、列