計算偏度與峰度（R）

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