1. ホーム
  2. python

[解決済み】Numpyのstd計算。TypeError: cannot perform reduce with flexible type.

2022-01-25 18:57:28

質問

7行目から始まる数値の行を読み、データがなくなるまで数値をリストにまとめ、このリストに対して標準偏差と%rmsを計算しようとしています。簡単そうに見えますが、エラーが出続けています。

Traceback (most recent call last):
  File "rmscalc.py", line 21, in <module>
    std = np.std(values)
  File "/System/Library/Frameworks/Python.framework/Versions/2.7/Extras/lib/python/numpy/core/fromnumeric.py", line 2817, in std
    keepdims=keepdims)
  File "/System/Library/Frameworks/Python.framework/Versions/2.7/Extras/lib/python/numpy/core/_methods.py", line 116, in _std
    keepdims=keepdims)
  File "/System/Library/Frameworks/Python.framework/Versions/2.7/Extras/lib/python/numpy/core/_methods.py", line 86, in _var
    arrmean = um.add.reduce(arr, axis=axis, dtype=dtype, keepdims=True)
TypeError: cannot perform reduce with flexible type

以下は私のコードです。

import numpy as np
import glob
import os

values = []
line_number = 6

road = '/Users/allisondavis/Documents/HCl'
for pbpfile in glob.glob(os.path.join(road, 'pbpfile*')): 
    lines = open(pbpfile, 'r').readlines()
    while line_number < 400 :
        if lines[line_number] == '\n':
            break
        else: 
            variables = lines[line_number].split()
            values.append(variables)
            line_number = line_number + 3
            print values

a = np.asarray(values).astype(np.float32)
std = np.std(a)
rms = std * 100
print rms

編集:rmsを生成しますが(これは間違っています - 理由はまだわかりません)、次のエラーメッセージは紛らわしいです。カウントを高くする必要があります(どんなに大きなファイルでも全部取得できるように400を選びました)。

 Traceback (most recent call last):
  File "rmscalc.py", line 13, in <module>
    if lines[line_number] == '\n':
IndexError: list index out of range

解決方法は?

values は文字列配列であり a . 変換する a を使用して数値型に変換します。 astype . 例えば

a = np.asarray(values).astype(np.float32)
std = np.std(a)