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
>
名称栏目:R语言笔记自定义峰度和偏度函数
分享路径:http://myzitong.com/article/gojgch.html