【转】偏度与峰度

?

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23

# test the skewness 以及峰度   library(moments) x <- rnorm(100,0,1) # 产生100均值为0， 标准差为1 的正太分布随机数 x <- rnorm(100,5,3)   skewness(x) # 计算偏度 agostino.test(x) # 偏度的检验 kurtosis(x) # 计算峰度 由于没有减3 ，所以正太应该为3 anscombe.test(x) # 峰度检验   head(data) # 查看数据 skewness(data) agostino.test(x) kurtosis(data) anscombe.test(data)     # 清空内存 ls() rm(list = ls()) gc() ls()