R语言使用moments包计算偏度(Skewness)和峰度(Kurtosis)实战计算偏度(Skewness)和峰度(Kurtosis)、确定样本数据是否具有与正态分布匹配的偏度和峰度(假设检验)。
直接百度公式自己手算吧。你也是够懒的,首先m是期望,s是标准差
g1<- n/((n-1)*(n-2))*sum((x-m)^3)/s^3g2<- ((n*(n+1))/((n-1)*(n-2)*(n-3))*sum((x-m)^4)/s^4-(3*(n-1)^2)/((n-2)*(n-3)))
第一个是峰度,第二个偏度
好消息,FBasics包中有,e1071中也有。我记得有人说算出来有微小差别
skewness()
kurtosis()