#常用统计分析
#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(mtcars
m
p
g
,
m
t
c
a
r
s
mpg,mtcars
mpg,mtcarshp)) #也不能出正确结果
sapply(mtcars[vars],f1) #sapply的语法sapply(x,fun,option),x是数据框,fun是任意函数,option是描述统计量(mean,sd,var,min,max,median,length,range,quantile)
mean(mtcars
m
p
g
)
;
s
d
(
m
t
c
a
r
s
mpg);sd(mtcars
mpg);sd(mtcarsmpg);var(mtcars
m
p
g
)
;
m
i
n
(
m
t
c
a
r
s
mpg);min(mtcars
mpg);min(mtcarsmpg);max(mtcars
m
p
g
)
;
m
e
d
i
a
n
(
m
t
c
a
r
s
mpg);median(mtcars
mpg);median(mtcarsmpg);length(mtcars
m
p
g
)
;
r
a
n
g
e
(
m
t
c
a
r
s
mpg);range(mtcars
mpg);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() #查找已加载的程序包
R常用描述性统计
最新推荐文章于 2024-08-04 01:00:00 发布