#常用统计分析
#1.描述性统计分析
mtcars #车辆路试数据集,每加仑汽油行驶英里数(mpg),马力(hp),车重(wt),变速箱类型(am,0自动挡,1手动挡,)气缸数cyl
head(mtcars) #只看前六行数据
summary(mtcars) #对每个变量输出最小值,最大值,均值,中位数,四分位数
vars=c(“mpg”,“hp”,“wt”);summary(mtcars[vars])
summary(mtcarsKaTeX parse error: Expected 'EOF', got '#' at position 6: mpg)
#̲基础安装包并没有提供偏度和峰度…mpg) #mpg的偏度和峰度
f1(mtcars[vars]) #不能出来结果,提示错误
f1(c(mtcarsmpg,mtcarshp)) #也不能出正确结果
sapply(mtcars[vars],f1) #sapply的语法sapply(x,fun,option),x是数据框,fun是任意函数,option是描述统计量(mean,sd,var,min,max,median,length,range,quantile)
mean(mtcarsmpg);sd(mtcarsmpg);var(mtcarsmpg);min(mtcarsmpg);max(mtcarsmpg);median(mtcarsmpg);length(mtcarsmpg);range(mtcarsmpg);quantile(mtcars$mpg)
#虽然基本的包中没有包含偏度,峰度等描述统计量,但有学者扩展了描述统计,形成了其它的包
install.packages(“pastecs”)
library(pastecs)
stat.desc(mtcars[vars]) #pastecs包中的函数
#还有psych包中有describe可以解出偏度和峰度
install.packages(“psych”)
library(psych)
describe(mtcars[vars])
sessionInfo() #查找已加载的程序包