相関係数の計算結果の保存

Question

for文で一番大きかった相関係数とそのときのiを記録したい。最後にプリント文で(最大値がでたときの)i年ずらした時の相関が最も強く、相関係数は(最大値（2年）の値）

コードimport numpy as np

data = np.loadtxt('test_data.txt')

import matplotlib.pyplot as plt

year = year.astype('int')
hares_tmp = hares[:-1]
carrots_tmp = carrots[1:]
year_tmp = np.arange(len(hares_tmp))

n = int(len(hares)/2)
carrots_tmp = carrots[:n]
for i in range(n):
    hares_tmp = hares[i:i+n]
    print(i,'年ずらした時の相関、相関係数は',np.min(np.corrcoef(hares_tmp,carrots_tmp)))    
```0 年ずらした時の相関、相関係数は -0.13383412475583126
1 年ずらした時の相関、相関係数は 0.48579256922062664
2 年ずらした時の相関、相関係数は 0.8905516204182647
3 年ずらした時の相関、相関係数は 0.7132762541264995
4 年ずらした時の相関、相関係数は 0.17636199533173733
5 年ずらした時の相関、相関係数は -0.029075254388799553
6 年ずらした時の相関、相関係数は -0.26913537526056697
7 年ずらした時の相関、相関係数は -0.5472898182646936
8 年ずらした時の相関、相関係数は -0.5611479259577979
9 年ずらした時の相関、相関係数は -0.42223540792610453







サンプルコード```# year	hare	lynx	carrot
1900	3.00E+04	4.00E+03	48300
1901	4.72E+04	6.10E+03	48200
1902	7.02E+04	9.80E+03	41500
1903	7.74E+04	3.52E+04	38200
1904	3.63E+04	5.94E+04	40600
1905	2.06E+04	4.17E+04	39800
1906	1.81E+04	1.90E+04	38600
1907	2.14E+04	1.30E+04	42300
1908	2.20E+04	8.30E+03	44500
1909	2.54E+04	9.10E+03	42100
1910	2.71E+04	7.40E+03	46000
1911	4.03E+04	8.00E+03	46800
1912	5.70E+04	1.23E+04	43800
1913	7.66E+04	1.95E+04	40900
1914	5.23E+04	4.57E+04	39400
1915	1.95E+04	5.11E+04	39000
1916	1.12E+04	2.97E+04	36700
1917	7.60E+03	1.58E+04	41800
1918	1.46E+04	9.70E+03	43300
1919	1.62E+04	1.01E+04	41300
1920	2.47E+04	8.60E+03	47300

質問はforのなかで計算した相関係数の値最大値とその時のiを保存させる方法はどうすればいいですか？というもので私はfor文のなかにif文でmaxの値を記録させその時のiも文字を置いて記録させようとしましたが失敗したので、、ここで質問させていただきました。

Comments

[toast-uz]

「質問」は何でしょうか？
また、'test_data.txtのサンプルがあると、回答しやすいです。
ここへのコメントでなく、質問そのものを修正ください。

[toast-uz]

修正ありがとうごさいます。想定して回答を作成していましたので、ご確認ください。

Answer 1

いろいろな求め方がありますが、例えば、一番下の3行（for以下）を、次の5行に書き換えてみてください。

corrs = []
for i in range(n):
    hares_tmp = hares[i:i+n]
    corrs.append(np.min(np.corrcoef(hares_tmp,carrots_tmp)))
    print(i,'年ずらした時の相関、相関係数は', corrs[i])
print(corrs.index(max(corrs)),'年ずらした時が最大の相関、相関係数は', max(corrs))