【转】偏度与峰度

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