python中偏度與峰度的實作

import numpy as np
from scipy import stats #scipy中的stats可以做統計推斷
np.set_printoptions(suppress=True) #将數字表示成我們習慣的方式
x = np.random.randn(100) #randn可随機生成标準正态分布的樣本，100表示樣本數量
mu = np.mean(x, axis=0) #axis=0，那麼輸出矩陣是1行，求每一列的平均；axis=1，輸出矩陣是1列，求每一行的平均
sigma = np.std(x, axis=0) #求标準差，這裡除的是N
skew = stats.skew(x) #求偏度
kurtosis = stats.kurtosis(x) #求峰度 

print(x,'\n')
print(mu,sigma,skew,kurtosis)