计算偏度与峰度(R)

微笑、不失礼 提交于 2020-03-08 22:23:21
# 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()

 

标签
易学教程内所有资源均来自网络或用户发布的内容,如有违反法律规定的内容欢迎反馈
该文章没有解决你所遇到的问题?点击提问,说说你的问题,让更多的人一起探讨吧!