1

R与库一起使用moments来生成一个小数据集并计算我的数据的四个初始时刻:

  1. 意思是
  2. 变化
  3. 偏度
  4. 峰态

代码如下所示。我为我的 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

人们可能会注意到,两种方法的偏度和峰度都是不同的。

我的问题是:为什么他们给出不同的结果?哪个结果是正确的?

4

0 回答 0