同じ形状の配列の統計を見つける

Question

こんにちは私はピアソンのR統計をいくつかのデータにマップしようとしています。しかし、データは配列の上にあります。私は配列を作っていますが、私が何をしても、同じ形と大きさではなく、エラーが発生します。
'複数の要素を持つ配列の真理値はあいまいです。 a.any（）またはa.all（） 'を使用してください

ご協力いただければ幸いです！

f = netCDF4.Dataset('C:/data/Feb.nc', 
         'r') 
    # Import variables 
    v = f.variables['store_Bio'][0:28, 0:24, 0, 0:60, 0:60])/81 
    t = f.variables['temp_a'][0:28, 0:24, 0, :, :] 
    g = f.variables['skin_temp'][0:28, 0:24, :, :] 
    cc = f.variables['cloud_cover'][0:28, 0:24, 0, :, :] 


    # Mean over array to make them the same shape 
    v1 = v.mean(axis=(0, 1)) 
    t1 = t.mean(axis=(0, 1)) 
    g1 = g.mean(axis=(0, 1)) 
    cc1 = cc.mean(axis=(0, 1)) 

    # Pearson R stats 
    p11 = pearsonr(v1, t1) 
    p21 = pearsonr(v1, g1) 
    p31 = pearsonr(v1, cc1) 

Answer 1

pearsonrの入力タイプは、配列のようにする必要があります。

>>> a = np.array([[1,2],[3,4]]) 
>>> b = np.array([[1,2],[3,4]]) 
>>> pearsonr(a,b) 
Traceback (most recent call last): 
    File "<stdin>", line 1, in <module> 
    File "/usr/lib/python2.7/dist-packages/scipy/stats/stats.py", line 2518, in pearsonr 
    r = max(min(r, 1.0), -1.0) 
ValueError: The truth value of an array with more than one element is ambiguous. Use a.any() or a.all() 

>>> pearsonr(a.ravel(),b.ravel()) 
(1.0, 0.0) 

対応する各値について、次のコードを試すことができます。あなたの入力のサイズに等しいresult配列を与えるでしょう。

>>> result = [pearsonr(x, y) for x, y in zip(a,b)]