python求解峰度kurtosis 和matlab求得的值差很远 解决办法
正在做一个把matlab程序转python的工作,遇到如题问题
实际的代码中输入数据有点大,下面简单的以 x = [1,2,3,4,5]
来举例说明
matlab
>> x = [1,2,3,4,5];
>> kurtosis(x)
ans =
1.7000
python
from scipy import stats
x = [1,2,3,4,5]
print(stats.kurtosis(x))
输出
-1.3
解决办法
调用 kurtosis 函数时指定不使用Fisher的定义,使用Pearson的定义,详细参数见 scipy.stats.kurtosis
from scipy import stats
x = [1,2,3,4,5]
print(stats.kurtosis(x, fisher=False))