R语言有两个包可以计算峰度和偏度

1、moments

实现：

library(moments);

x=rnorm(10000)

偏度

 skewness(x)    或者    moments:skewness(x)

峰度

kurtosis(x)     或     moments::kurtosis(x)

2、fBasics

library(fBasics);

x=rnorm(10000)

偏度

 skewness(x) 

峰度

kurtosis(x)     

example

> x=rnorm(1000000)
> moments::kurtosis(x)
[1] 3.002087
> kurtosis(x)
[1] 0.002080558
attr(,"method")
[1] "excess"
> skewness(x)
[1] -0.003034968
attr(,"method")
[1] "moment"
> moments::skewness(x)
[1] -0.003034973

区别是moment包的峰度没有减去3