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)