R语言笔记自定义峰度和偏度函数-创新互联

> mtcars<- function(x,na.omit=FALSE){

+ if(na.omit)

+ x<-x[!is.na(x)]

+ m<-mean(x)

+ n<-length(x)

+ s<-sd(x)

+ skew<-sum((x-m)^3/s^3)/n

+ kurt<-sum((x-m)^4/s^4)/n-3

+ return(c(n=n,mean=m,stdev=s,skew=skew,kurtosis=kurt))

+ }

vars<- c("mpg","hp","wt")

> sapply(mtcars[vars],mtstat)

        mpg      hp      wt

n     32.000000  32.0000000 32.00000000

mean   20.090625 146.6875000  3.21725000

stdev   6.026948  68.5628685  0.97845744

skew    0.610655  0.7260237  0.42314646

kurtosis -0.372766  -0.1355511 -0.02271075

>

另外有需要云服务器可以了解下创新互联cdcxhl.cn，海内外云服务器15元起步，三天无理由+7*72小时售后在线，公司持有idc许可证，提供“云服务器、裸金属服务器、高防服务器、香港服务器、美国服务器、虚拟主机、免备案服务器”等云主机租用服务以及企业上云的综合解决方案，具有“安全稳定、简单易用、服务可用性高、性价比高”等特点与优势，专为企业上云打造定制，能够满足用户丰富、多元化的应用场景需求。