描述性统计分析

用sapply 函数计算所选择的任意描述性统计量。

sapply（x,FUN,optins）, 其中x是数据框或矩阵，FUN为一个任意的函数。如果指定了option，他们将被传递给FUN。

例：mystats<-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,kurt=kurt))

}

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

> options(digits=2)
> sapply(mtcars[myvars],mystats)
       mpg     hp    wt
n    32.00  32.00 32.00
mean 20.09 146.69  3.22
stdev    6.03  68.56  0.98
skew  0.61   0.73  0.42
kurt  2.90   3.16  3.29