Pythonによる最小二乗法 配列の最小値とそれに対応する番号を表示する。

Question

pythonを使って最小二乗法のプログラムを書こうとしています。
乱数を使って適当なデータ(1)を作成しそれを
y=a*sinx+1…(2)
でフィッティングし、最も良いaの値を決定したいです。
残差の2乗(epsilon)を計算するところまではできました。

b = np.random.normal(0,1,11)
data = 2*(1+b*0.1)*np.sin(np.radians(deg))+1  #乱数を使ったデータ(1)
deg = np.linspace(0,90,11)
for i in range(1,10):
    sin_i = (i*0.1+1.5)*np.sin(np.radians(deg))+1#a*sinx+1として1.5≦a≦2.5まで0.1刻みでaの値を変えていく…(2)
    epsilon = ((data) - (sin_i))**2#データとsin_iの差の2乗の和をとる。
for j in range(0,10):
    sum_all =np.sum(epsilon[j])#各iに対する配列epsilonの和をとる
    print(sum_all)#各iに対するepsirlonの和を表示する。
min(sum_all)#最小になるepsilonを表示する。

しかし、最後の行の

min(sum_all)

で、TypeError: 'numpy.float64' object is not iterable
とうエラーが出てしまいます。

sum_allの最小値を求める方法とその最小値が何番目のiなのか
(i,sum_all)の形で求める方法を教えていただきたいです。

Answer 1

【追記】
質問のコードですとepsilon[j]が既にnumpy.float64になっていますね。

---

下記によりsum_allはnumpy.float64になっているのではないですか？

for j in range(0,10):
    sum_all =np.sum(epsilon[j])#各iに対する配列epsilonの和をとる

Comments

[kunikuni96]

お礼が遅くなり申し訳ありません。
回答ありがとうございました。

いただいた回答をもとに、not iterable やfloatという意味を調べてどうして動かないかを理解することができました。

sum_allが反復可能な形になっておらず(epsilon[j]の時点で)、floatになっているということだと理解しました。