質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

ただいまの
回答率

88.82%

Python 期待通りの計算結果になりません。正しい計算ロジックを教えてください。

解決済

回答 2

投稿

  • 評価
  • クリップ 0
  • VIEW 293

moto_i

score 20

行列演算に関する質問です。下記のリストがあります。
list = [[7, 2, 1, 3], [15, 4, 1, 5], [20, 5, 1, 10]]

各行の要素を各行の先頭の値と末尾の値を足した合計で除算した後100倍にしたいです。
[7, 2, 1, 3] であれば 先頭の値が7,末尾の値が3なので、合計が10となり
各要素の値は10で割られて
0.7, 0.2, 0.1, 0.3 となります。
これを100倍して
70, 20, 10, 30というような計算処理をしたいです。
が思ったように計算ができません。
正しいロジックを教えてください。

list = [[7, 2, 1, 3], [15, 4, 1, 5], [20, 5, 1, 10]]

calc_list = [[0]*4]*3
print(f"calc_list : {calc_list}")
for i in range(3):
    for j in range(4):
        calc_list[i][j] = ( list[i][j] / (list[i][0] + list[i][3]) )*100
print(f"calc_list : {calc_list}")


実行結果
calc_list : [[80.0, 16.666666666666664, 3.3333333333333335, 20.0], [80.0, 16.666666666666664, 3.3333333333333335, 20.0], [80.0, 16.666666666666664, 3.3333333333333335, 20.0]]
```
期待する結果
calc_list : [[70.0, 0.20, 0.10, 0.30], [75.0, 0.20, 0.05, 0.25], [66.6666, 16.6666, 3.3333, 33.3333]]

  • 気になる質問をクリップする

    クリップした質問は、後からいつでもマイページで確認できます。

    またクリップした質問に回答があった際、通知やメールを受け取ることができます。

    クリップを取り消します

  • 良い質問の評価を上げる

    以下のような質問は評価を上げましょう

    • 質問内容が明確
    • 自分も答えを知りたい
    • 質問者以外のユーザにも役立つ

    評価が高い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。

    質問の評価を上げたことを取り消します

  • 評価を下げられる数の上限に達しました

    評価を下げることができません

    • 1日5回まで評価を下げられます
    • 1日に1ユーザに対して2回まで評価を下げられます

    質問の評価を下げる

    teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。

    • プログラミングに関係のない質問
    • やってほしいことだけを記載した丸投げの質問
    • 問題・課題が含まれていない質問
    • 意図的に内容が抹消された質問
    • 過去に投稿した質問と同じ内容の質問
    • 広告と受け取られるような投稿

    評価が下がると、TOPページの「アクティブ」「注目」タブのフィードに表示されにくくなります。

    質問の評価を下げたことを取り消します

    この機能は開放されていません

    評価を下げる条件を満たしてません

    評価を下げる理由を選択してください

    詳細な説明はこちら

    上記に当てはまらず、質問内容が明確になっていない質問には「情報の追加・修正依頼」機能からコメントをしてください。

    質問の評価を下げる機能の利用条件

    この機能を利用するためには、以下の事項を行う必要があります。

回答 2

checkベストアンサー

+1

数値計算は numpy を使いましょう。

import numpy as np

a = np.array([[7, 2, 1, 3], [15, 4, 1, 5], [20, 5, 1, 10]])

s = a[:, 0] + a[:, -1]
out = a / s[:, np.newaxis] * 100
print(out)
# [[70.         20.         10.         30.        ]
#  [75.         20.          5.         25.        ]
#  [66.66666667 16.66666667  3.33333333 33.33333333]]

解説

1. 最初の列目と最後の列を足す。

s = a[:, 0] + a[:, -1]
print(s)  # [10 20 30]

2. 配列 s の形状を (3,) から (3, 1) にする。

s[:, np.newaxis]
[10 20 30]


[[10]
 [20]
 [30]]

3. 除算する際にブロードキャストが行われる。

左辺 a が (3, 4) の配列、右辺 s[:, np.newaxis] が (3, 1) の配列なので、
右辺が (3, 1) -> (3, 4) に broadcast されてから除算される。

[[10]
 [20]
 [30]]


array([[10, 10, 10, 10],
       [20, 20, 20, 20],
       [30, 30, 30, 30]])


としてから、配列 a と上の配列の要素ごとの除算

4. 除算後に各要素を 100 倍する。

a / s[:, np.newaxis] * 100

division by zero はおこならい (最初の列目+最後の列=0ではない)ことを仮定しています。

zero 除算対策版

最初の列目+最後の列=0となることがある場合は、割れないので値を0とする場合

import numpy as np

a = np.array([[7, 2, 1, 3], [15, 4, 1, 5], [20, 5, 1, 10]], dtype=float)

s = (a[:, 0] + a[:, -1])[:, np.newaxis]
out = np.divide(a, s, out=np.zeros_like(a), where=s != 0)
print(out)
# [[0.7        0.2        0.1        0.3       ]
#  [0.75       0.2        0.05       0.25      ]
#  [0.66666667 0.16666667 0.03333333 0.33333333]]

投稿

編集

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

  • 2019/04/06 17:19

    回答ありがとうございました。

    キャンセル

0

質問文に記載されている処理を, ほぼそのままコードにしてみました。

d.py

vals = [[7, 2, 1, 3], [15, 4, 1, 5], [20, 5, 1, 10]]

divs = [x[0] + x[-1] for x in vals]
result = [[e * 100 / divs[idx] for e in vals[idx]] for idx in range(len(divs))]

print(*result, sep="\n")

実行結果

イメージ説明

投稿

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

  • 2019/04/07 20:39

    回答ありがとうございました

    キャンセル

15分調べてもわからないことは、teratailで質問しよう!

  • ただいまの回答率 88.82%
  • 質問をまとめることで、思考を整理して素早く解決
  • テンプレート機能で、簡単に質問をまとめられる

関連した質問

同じタグがついた質問を見る