Pythonで統計-モーメント

モーメント

１次モーメント １次モーメントは、平均を表します。 import numpy as np vals = np.random.normal(0, 0.5, 10000) #１次モーメント　平均 moment1 = np.mean(vals) ２次モーメント

２次モーメントは、分散を表します。

import numpy as np

vals = np.random.normal(0, 0.5, 10000)
#２次モーメント　分散
moment2 = np.var(vals)

３次モーメント

３次モーメントは、歪度を表します。
分布がどれだけ偏っているかを示します。
scipy.statsをインポートします。

import numpy as np
import scipy.stats as sp

vals = np.random.normal(0, 0.5, 10000)

#３次モーメント　歪度　分布がどれだけ偏っているか

moment3 = sp.skew(vals)

４次モーメント

４次モーメントは、尖度を表します。
分布がどれだけ尖っているかを示します。
3次モーメント同様にscipy.statsをインポートします。

import numpy as np
import scipy.stats as sp

vals = np.random.normal(0, 0.5, 10000)

#４次モーメント　尖度　どれだけ尖っているか

moment4 = sp.kurtosis(vals)

まとめ

１から４までのモーメントを出力します。
ヒストグラムをおまけに表示します。

import numpy as np
import matplotlib.pyplot as plt
import scipy.stats as sp

vals = np.random.normal(0, 0.5, 10000)

plt.hist(vals, 50)
plt.show()

moment = []
#１次モーメント　平均
moment.append(np.mean(vals))

#２次モーメント　分散
moment.append(np.var(vals))

#３次モーメント　歪度　分布がどれだけ偏っているか

moment.append(sp.skew(vals))

#４次モーメント　尖度　どれだけ尖っているか

moment.append(sp.kurtosis(vals))

for i in range(4):
   print("{0}モーメント：{1}".format(i+1,moment[i]))

参考