# 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) # 峰度检验 data <- read.csv(file = "C:\\users\\mike1\\desktop\\gaussanData.csv") # 读取数据 head(data) # 查看数据 skewness(data) agostino.test(x) kurtosis(data) anscombe.test(data) # 清空内存 ls() rm(list = ls()) gc() ls()
来源:https://www.cnblogs.com/zijidefengge/p/12445145.html