python/データの季節調整/アドバイスをお願いします

Question

データ（各月ごとに５年分）の季節調整を行い，調整結果
をすべて出力するプログラムを作成しなさい
上の問題を授業で出題され、以下のプログラムを作成しましたが、改良すべき点（こうしたらもっと見やすい、コンパクトになる）がありましたら些細なことでいいので指摘・アドバイスをお願いします。
プログラミングは大学で初めて習いましたが、ある程度は得意になりたいと考えています。Colaboratoryを用いています。
また、numpyを用いれば簡単に計算ができるそうですが、今回は使わない形でお願いします。

---

授業で与えられた情報を記載します
季節調整
季節変動（毎年のパターン）を抽出し調整する。
［季別調整法］

i年目、j月のデータを aij とする。
月平均 Mj
月平均の平均 M
調整のための指数 Sj = Mj / M
調整結果 aij / Sj

加点対象
• 調整後のデータで最大値や最大値となった年月を表示する
• 別の要素数の異なるリストが入力されても計算できる

データ
40.0 37.7 51.8 45.9 45.2 45.4 62.0 40.8 42.1 49.1 51.5 93.4
42.6 39.5 53.4 47.6 47.0 47.5 63.5 42.2 42.9 51.1 51.8 94.7
43.3 39.8 53.6 48.5 48.3 47.5 65.9 42.8 43.7 53.2 52.9 97.0
44.9 42.3 55.0 50.7 48.6 48.9 68.8 43.3 46.2 54.2 54.6 100.5
47.5 42.8 57.7 52.6 51.0 50.9 71.4 45.0 47.3 56.8 56.7 102.7

---

作成したコード

a1=[40.0, 37.7, 51.8, 45.9, 45.2, 45.4, 62.0, 40.8, 42.1, 49.1, 51.5, 93.4]
a2=[42.6, 39.5, 53.4, 47.6, 47.0, 47.5, 63.5, 42.2, 42.9, 51.1, 51.8, 94.7]
a3=[43.3, 39.8, 53.6, 48.5, 48.3, 47.5, 65.9, 42.8, 43.7, 53.2, 52.9, 97.0]
a4=[44.9, 42.3, 55.0, 50.7, 48.6, 48.9, 68.8, 43.3, 46.2, 54.2, 54.6, 100.5]
a5=[47.5, 42.8, 57.7, 52.6, 51.0, 50.9, 71.4, 45.0, 47.3, 56.8, 56.7, 102.7]
a =[a1,a2,a3,a4,a5]

b1=[]
b2=[]
b3=[]
b4=[]
b5=[]
b=[b1,b2,b3,b4,b5]

X=0
S=[]
Max_list=[]

for num in range(5):
      X=X+(sum(a[num]))
X=X/60

for y in range(5):
  for m in range(12):
    M=(sum(x[m] for x in a))/5
    S.append(M/X)
    b[y].append(round(a[y][m]/S[m],1))
  print("b",y+1,"=",b[y])
  Max_list.append(max(b[y]))
  y=y+1

Max_num=(max(Max_list))
print("最大値は",max(Max_list))
Max_y=Max_list.index(Max_num)
print(Max_list.index(Max_num)+1,"年目の",b[Max_y].index(Max_num)+1,"ヶ月目に最大値が出ました。")

Answer 1

動いてるので無理に変えなくてもいい気はしますが、気になったところをいくつか修正してコメントつけました。ご参考まで。

a1 = [40.0, 37.7, 51.8, 45.9, 45.2, 45.4, 62.0, 40.8, 42.1, 49.1, 51.5, 93.4]
a2 = [42.6, 39.5, 53.4, 47.6, 47.0, 47.5, 63.5, 42.2, 42.9, 51.1, 51.8, 94.7]
a3 = [43.3, 39.8, 53.6, 48.5, 48.3, 47.5, 65.9, 42.8, 43.7, 53.2, 52.9, 97.0]
a4 = [44.9, 42.3, 55.0, 50.7, 48.6, 48.9, 68.8, 43.3, 46.2, 54.2, 54.6, 100.5]
a5 = [47.5, 42.8, 57.7, 52.6, 51.0, 50.9, 71.4, 45.0, 47.3, 56.8, 56.7, 102.7]
a  = [a1,a2,a3,a4,a5]

X = 0
S = []
Max_list = []

for ai in a:  # 簡略化
      X +=  sum(ai)
X = X/60

b  =  []
for y in range(5):
    b.append([])  # b1,b2,..の代わり
    for m in range(12):
        M = (sum(x[m] for x in a))/5
        S.append(M/X)
        b[y].append(round(a[y][m]/S[m],1))
    print("b",y+1,"=",b[y])
    Max_list.append(max(b[y]))
    #y=y+1  for文では不要！

Max_num=(max(Max_list))
print("最大値は",max(Max_list))
Max_y=Max_list.index(Max_num)
print(Max_list.index(Max_num)+1,"年目の",b[Max_y].index(Max_num)+1,"ヶ月目に最大値が出ました。")
# python3.6以降なら、f文字列を使うと見やすいです。例↓
# print(f"{Max_list.index(Max_num)+1}年目の{b[Max_y].index(Max_num)+1}ヶ月目に最大値が出ました。")

Answer 2

他にもありますが、とりあえずは以下の二つを直すと良いと思います。あと、途中でround関数を用いているのはなぜでしょうか。

###(1) 不要な変数(a1,a2,a3,a4,a5やb1,b2,b3,b4,b5)は定義しない方が良いです。

以下のようにして書いた方が良いと思います。また、bについてはリスト内包表記と呼ばれるものなので、調べてみてください。

a = [[40.0, 37.7, 51.8, 45.9, 45.2, 45.4, 62.0, 40.8, 42.1, 49.1, 51.5, 93.4],
     [42.6, 39.5, 53.4, 47.6, 47.0, 47.5, 63.5, 42.2, 42.9, 51.1, 51.8, 94.7],
     [43.3, 39.8, 53.6, 48.5, 48.3, 47.5, 65.9, 42.8, 43.7, 53.2, 52.9, 97.0],
     [44.9, 42.3, 55.0, 50.7, 48.6, 48.9, 68.8, 43.3, 46.2, 54.2, 54.6, 100.5],
     [47.5, 42.8, 57.7, 52.6, 51.0, 50.9, 71.4, 45.0, 47.3, 56.8, 56.7, 102.7]]

b = [[] for _ in range(5)]

###(2) それぞれの変数が何を指すかわかりにくいです。

変数の定義の際にコメントアウトで示しておくと良いです。aであれば以下のようです。

# i年目、j月のデータを aij とする
a = [[40.0, 37.7, 51.8, 45.9, 45.2, 45.4, 62.0, 40.8, 42.1, 49.1, 51.5, 93.4],
     [42.6, 39.5, 53.4, 47.6, 47.0, 47.5, 63.5, 42.2, 42.9, 51.1, 51.8, 94.7],
     [43.3, 39.8, 53.6, 48.5, 48.3, 47.5, 65.9, 42.8, 43.7, 53.2, 52.9, 97.0],
     [44.9, 42.3, 55.0, 50.7, 48.6, 48.9, 68.8, 43.3, 46.2, 54.2, 54.6, 100.5],
     [47.5, 42.8, 57.7, 52.6, 51.0, 50.9, 71.4, 45.0, 47.3, 56.8, 56.7, 102.7]]