Np.corrcoefが行列を返すのは奇妙に思えます。
correlation1 = corrcoef(Strategy1Returns,Strategy2Returns)
[[ 1. -0.99598935]
[-0.99598935 1. ]]
なぜこれが事実であり、古典的な意味で1つの値のみを返すことができるかどうかを知っていますか?
これにより、2つ以上のデータセットの相関係数を計算できます。
>>> from numpy import *
>>> a = array([1,2,3,4,6,7,8,9])
>>> b = array([2,4,6,8,10,12,13,15])
>>> c = array([-1,-2,-2,-3,-4,-6,-7,-8])
>>> corrcoef([a,b,c])
array([[ 1. , 0.99535001, -0.9805214 ],
[ 0.99535001, 1. , -0.97172394],
[-0.9805214 , -0.97172394, 1. ]])
ここで、a、b(0.995)、a、c(-0.981)およびb、c(-0.972)の相関係数を一度に取得できます。 2データセットの場合は、Nデータセットクラスの特別な場合にすぎません。そして、おそらく同じ戻り値の型を保持する方が良いでしょう。 「1つの値」は単純に取得できるため、
>>> corrcoef(a,b)[1,0]
0.99535001355530017
特別なケースを作成する大きな理由はありません。
corrcoef
は、正規化された共分散行列を返します。
共分散行列は行列です
Cov( X, X ) Cov( X, Y )
Cov( Y, X ) Cov( Y, Y )
正規化すると、次のマトリックスが生成されます。
Corr( X, X ) Corr( X, Y )
Corr( Y, X ) Corr( Y, Y )
correlation1[0, 0 ]
はStrategy1Returns
およびそれ自体、1でなければなりません。ただcorrelation1[ 0, 1 ]
。
相関行列は、任意の有限数の変数間の相関を表す標準的な方法です。 [〜#〜] n [〜#〜]データベクトルの相関行列は対称です[〜#〜] n [〜#〜]×[〜#〜] n [〜#〜]対角線が1の行列。 [〜#〜] n [〜#〜] = 2の場合のみ、このマトリックスには1つの自由パラメーターがあります。
関数numpyの相関は、相関させたい2つの1D配列で機能し、1つの相関値を返します。
Matplotlib.cbookの使用を検討してください
例えば:
import matplotlib.cbook as cbook
segments = cbook.pieces(np.arange(20), 3)
for s in segments:
print s
次の関数を使用して、相関係数のみを返すことができます。
def pearson_r(x, y):
"""Compute Pearson correlation coefficient between two arrays."""
# Compute correlation matrix
corr_mat = np.corrcoef(x, y)
# Return entry [0,1]
return corr_mat[0,1]