允许a = c(1:9,1) library(timeDate)（也在library(fBasics)中导入），skewness(a) = 0.10901和kurtosis(a) = -1.634，。
在library(moments)，skewness(a) = 0.1277和kurtosis(a) - 3 = -1.313，与Matlab的结果相符。
这是其中一个包中的错误吗？
答案 0 :(得分：9)
有很多方法可以估算人口时刻。您可以使用e1071包中的偏度和峰度函数来比较这些不同的方法
> library(e1071)
> skewness(a,,1)
[1] 0.1277249
> skewness(a,,2)
[1] 0.1514631
> skewness(a,,3)
[1] 0.1090534
> kurtosis(a,,1)
[1] -1.313042
> kurtosis(a,,2)
[1] -1.356985
> kurtosis(a,,3)
[1] -1.633564
所有这些措施都是有效的，无偏见的（至少在正常情况下）。根据D. N. Joanes和C. A. Gill（1998）（比较样本偏度和峰度的测量，统计学家，47,183-189。），类型3在非正常人群中表现更好，而2型在正常人群中表现最佳。类型1是较旧的定义。
在大样本中，这些估计量之间的差异可以忽略不计。看？所用公式的偏度。