天天看點

描述性統計分析

用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