用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