我正在使用R
和库moments
来生成一个小数据集并计算数据的四个初始时刻:
- 意味 着
- 变化
- 偏度
- 峰度
代码如下所示。我为我的 PRNG 设置了一个随机种子,并使用正态分布生成 1000 个数据点。
然后,我以两种方式打印四个时刻。首先,我打印,然后单独打印。然后,我使用方法all.moments打印它们。
library(moments)
set.seed(123)
x = rnorm(1000, sd = 0.02)
print(mean(x));
print(var(x));
print(skewness(x))
print(kurtosis(x))
print(moments::all.moments(x, order.max = 4))
输出如下所示。
print(mean(x));
0.0003225573
print(var(x));
0.0003933836
print(skewness(x));
0.06529391
print(kurtosis(x));
2.925747
print(moments::all.moments(x, order.max = 4));
1.000000e+00 3.225573e-04 3.930942e-04 8.889998e-07 4.527577e-07
人们可能会注意到,两种方法的偏度和峰度都是不同的。
我的问题是:为什么他们给出不同的结果?哪个结果是正确的?
请注意,第三和第四个时刻不是偏度和峰度。这些应该在之后计算