# 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()