我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
人们可能会注意到,两种方法的偏度和峰度都是不同的。
我的问题是:为什么他们给出不同的结果?哪个结果是正确的?