• 设为首页
  • 点击收藏
  • 手机版
    手机扫一扫访问
    迪恩网络手机版
  • 关注官方公众号
    微信扫一扫关注
    公众号

R语言笔记005——计算描述性统计量

原作者: [db:作者] 来自: [db:来源] 收藏 邀请

数据的分布特征:

  • 分布的集中趋势,反应各数据向其中心值靠拢或聚集的程度(平均数,中位数,四分位数,众数)
  • 分布的离散程度,反应各数据远离其中心值的趋势(极差,四分位差,方差,标准差,离散系数)
  • 分布的形状,反应数据分布的偏斜程度和峰度(偏态系数,峰度系数)

#######################

平均数(均值):一组数据相加后除以数据的个数而得到结果,称为平均数(mean)

中位数一组数据排序后处于中间位置上的变量值,称为中位数(median)

四分位数一组数据排序后处于25%(下四分位数)和75%(上四分位数)位置上的值,称为四分位数

  • 先计算位置,然后计算四分位数的值。50%处即为中位数

众数:一组数据中出现频数最多的数值(mode)

#######################

极差(全距):一组数据的最大值与最小值之差,称为极差(range)

四分位差:上四分位数与下四分位数之差,称为四分位差

平均差:各变量值与其平均数离差绝对值的平均数,称为平均差

方差:各变量值与其平均数离差平方的平均数,称为方差(variance)

标准差:方差的平方根称为标准差()

离散系数

########################

偏态:数据分布的不对称性,称为偏态

峰态:数据分布的平峰或尖峰程度,称为峰态

#########################

描述性统计量:主要包括分布的集中程度,分布的离散程度和分布的偏斜程度。

方法一:summary()函数——最大值,最小值,四分位数(上,下),均值

summary(mtcars) #结果有以下几条

1 vars<-c(\'mpg\',\'hp\',\'wt\')
2 head(mtcars[vars])

summary(mtcars[vars])

方法二:sapply()函数 #

 1 mystats<-function(x,na.omit=FALSE){
 2   if(na.omit)
 3       x<-x[!is.na(x)]
 4   m<-mean(x)
 5   n<-length(x)
 6   s<-sd(x)
 7   skew<-sum((x-m)^3/s^3)/n
 8   kurt<-sum((x-m)^4/s^4)/n
 9   return(c(n=n,mean=m,stdev=s,skew=skew,kurtosis=kurt))
10 }
11 sapply(mtcars[vars], mystats)

方法三:通过Hmisc包中的describe()函数——变量和观测的数量,缺失值和唯一值的数目、平均数,分位数,以及5个最大的值,和5个最小的值

1 library(Hmisc)
2 describe(mtcars[vars])

 

方法四:通过pastecs包中的stat.desc()函数计算描述性统计量

stat.desc(x,basic=TRUE,desc=TRUE,norm=FALSE,P=0.95)

x是一个数据框或时间序列。

basic=TRUE(默认值),则计算其中所有值、空值、缺失值的数量,以及最小值、最大值、值域,还有总和。

desc=TRUE(默认值),则计算中位数、平均数、平均数的标准误、平均数置信度为95%的置信区间、方差、标准差以及变异系数。

norm=TRUE(非默认),则返回正态分布统计量,包括偏度和峰度(以及它们的统计显著程度)和Shapiro–Wilk正态检验结果。

                                  这里使用了p值来计算平均数的置信区间(默认置信度为0.95)。

1 library(pastecs)
2 stat.desc(mtcars[vars])

 

方法五:通过psych包中的describe()函数——非缺失值的数量、平均数、标准差、中位数、截尾均值、绝对中位差、最小值、最大值、值域、偏度、峰度和平均值的标准误

1 library(psych)
2 describe(mtcars[vars])

以上方法都是为整体的数据计算描述性统计量


鲜花

握手

雷人

路过

鸡蛋
该文章已有0人参与评论

请发表评论

全部评论

专题导读
上一篇:
R语言与数据分析发布时间:2022-07-18
下一篇:
模型验证方法——R语言发布时间:2022-07-18
热门推荐
热门话题
阅读排行榜

扫描微信二维码

查看手机版网站

随时了解更新最新资讯

139-2527-9053

在线客服(服务时间 9:00~18:00)

在线QQ客服
地址:深圳市南山区西丽大学城创智工业园
电邮:jeky_zhao#qq.com
移动电话:139-2527-9053

Powered by 互联科技 X3.4© 2001-2213 极客世界.|Sitemap