テキストファイルを読み込み分散を行う方法

Question

テキストファイルを読み込んで分散を行う方法が分かりません。
def文を使い、二乗したものの合計をsum2(a,b)でs2に代入し、s2をデータ数で割ったものをa2で代入し、
分散をvと定義してv=a2-a*aで表す方法を教えてください。

import numpy as np
def sum(a,b):
    y=0
    for i in range(b):
        y+=a[i]
    return y
x=np.loadtxt('data.txt')
s=sum(x,x.size)
a=s/x.size

s2= sum2(a,b)

v=a2-a*a
print('number of data=',x.size,'sum=',s,'average=',a'variance=',v)

Comments

[jeanbiego]

・data.txtの内容はどういった構造ですか。
・記載しているコードを動かすとどうなりますか。
・sumという関数は組み込み関数として既に存在しているので、違う名前にしたほうが面倒がありません。

[shu-..]

data.txtはそのファイルにいくつかの数字が入っているものです。
この状態で起動すると
print('number of data=',x.size,'sum=',s,'average=',a'variance=',v)
^
SyntaxError: invalid syntax
と表示されます

[jeanbiego]

・ここに返答するだけでなく、質問欄を編集してください
・print(x)してみると、どう出ますか
・aと'variance='の間にカンマが抜けてませんか

Answer 1

numpy使っているのだから、numpyの機能を使ってみました。
pythonはfor文を使うと実行時間が長くなるので、可能な限り使わないほうが賢明です。

print文で変数を表示したい場合はformat文を使います。

import numpy as np
x = np.loadtxt('data.txt')
s = x.sum()
a = x.mean()
v = x.var()
print('number of data={0} sum={1} average={2} variance={3}'.format(x.size, s, a, v))

data.txtとしては以下のような形式を仮定しました。

1
5
6

実行結果です。

number of data=3 sum=12.0 average=4.0 variance=4.666666666666667

for文を使った書き方は推奨できませんが、念の為。

import numpy as np
def sum_array(a):
    y=0
    for i in a:
        y+=i
    return y
def sum2_array(a):
    y=0
    for i in a:
        y+=i*i
    return y

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

# 総和
s = sum_array(x)
# 平均
a = s/x.size

s2 = sum2_array(x)
# 二乗平均
a2 = s2/x.size
# 分散
v = a2-a*a
print('number of data={0} sum={1} average={2} variance={3}'.format(x.size, s, a, v))