R语言中的数据汇总

2016-08-19 04:06 阅读(445)    评论(0)   

问题:

在R中,如何对分组的数据进行汇总(统计其均值、方差等)?

指南:

下文中将介绍三种可用于对根据特定变量分组后的数据分别调用汇总函数(例如均值、标准差等等)的方法。

  • ddply()函数:这是最为简便的方法,需要调用plyr包,我们在此大力推荐大家使用。
  • summarizeBy() 函数:这也是个相对简单的方法,需要调用doBy包。
  • aggregate()函数:这个函数内置于R的原封安装中,相对来说使用起来会较为复杂一点。

假设我们拥有如下数据,根据数据中样本性别和条件的不同,我们可以将其分为F-placeboF-aspirinM-placeboM-aspirin四个组别。我们的任务是,分别统计每个组别中的样本个数、变量change的均值、标准差以及其均值标准误。

1

使用ddply:

2

缺失值处理

如果数据中存在确实值NA,那我们就需要对每一个函数设置参数na.rm=TRUE。函数length()没有内置na.rm选项,所以对于这种情况我们需要使用sum(!is.na(...))语句,来确认数据集中有多少的非缺失值。

3
一个用于样本计数以及计算均值、标准差、均值标准误和置信区间的函数

有时我们不需要像上文所述方法一样手动的取指定我们需要用于标准差计算的值,构建一个函数便能帮我们搞定以下所有的事:

  • 计算均值、计算标准差、样本计数
  • 计算均值标准误(友情提示,这个值在计算过程中不会用到组内变量。 更多详情见../../Graphs/Plotting means and error bars (ggplot2) )
  • 计算95%的置信区间(或是其他你所需求的数值)
  • 对列进行重命名从而使结果数据集更便于进一步的操作

为了使用这个函数,我们需要在进行计算前先把下面的代码粘贴到自己的程序中:

45

运用实例:summarySE函数帮你一步搞定

6

用NA来填充数据集中为空的组合情况

有些时候,汇总数据框中会存在一些空的因素组合情况,也就是说,这些组合理论上是可以存在的,但是并没有实际出现在我们的原始数据集中。有时,我们需要在数据框中用NA值对将这样的组合对应的内容进行替换填充。要完成这一点,我们需要在调用ddplysummarySE的相关语句时设置参数.drop=FALSE

运用实例:6.1

使用

使用summarizeBy()对数据进行折叠汇总:7

需要注意的是,这样的均值标准误并不适用于组内变量的比较 更多详情见../../Graphs/Plotting means and error bars (ggplot2)

缺失值处理

如果数据中存在缺失值NA,那我们就需要对每一个函数设置参数na.rm=TRUE。函数length()没有内置na.rm选项,所以对于这种情况我们需要使用sum(!is.na(...))语句,来确认数据集中有多少的非缺失值。8

一个用于样本计数以及计算均值、标准差、均值标准误和置信区间的函数

和对ddply包的使用类似的,通过组合运用doBy包中的函数,我们也可以构建一个类似的函数summarySE,帮我们搞定以下所有的事:

  • 计算均值、计算标准差、样本计数
  • 计算均值标准误(友情提示,这个值在计算过程中不会用到组内变量。 更多详情见../../Graphs/Plotting means and error bars (ggplot2) )
  • 计算95%的置信区间(或是其他你所需求的数值)
  • 对列进行重命名从而使结果数据集更便于进一步的操作

为了使用这个函数,我们需要在进行计算前先把下面的代码粘贴到自己的程序中:9

10

运用实例:summarySE函数帮你一步搞定

11
 用NA来填充数据集中为空的组合情况

有些时候,汇总数据框中会存在一些空的因素组合情况,也就是说,这些组合理论上是可以存在的,但是并没有实际出现在我们的原始数据集中。有时,我们需要在数据框中用NA值对将这样的组合对应的内容进行替换填充。

下面的函数将帮助我们将0值填充入所有的空缺组合对应的内容:

12运用实例:13

使用

aggregate函数使用起来可能相对繁琐一些,但由于它是R自带的内置函数,因而我们在使用它之前不必安装其他的程序包。141516

如果你想要无视数据集中存在的NA值,请设置参数na.rm=TRUE,如下所示:

17

公司二维码

 

 

分享到:

相关推荐

  • 12

    史上最腹黑病毒,18个月感染世界1/4人口,搞死3000万人,至今仍在变异!

    第一次世界大战是欧洲历史上破坏性最强的战争之一。 这场大战中,人类创造了许多个“首次”,如首次坦克战、首次毒气战、首次重机枪屠杀战以及超级重炮的首次登场。 这些深深反映出人类在战争领域的“强大创造力”。 马克沁重机枪 历史书为了世人牢记这些血和泪筑成的“创造力”,总不吝篇幅地...

  • vangogh_color and black and white

    色盲曾是超能力?人类远祖靠它度过了恐龙时代,如今却成了最高发的遗传病

    “原来你是红绿色盲(色弱)啊?!” “那你看国旗是什么颜色的?还有宝强的帽子是什么颜色的?” 这是一个色觉认知障碍者常常要面对的情景。 天知道他在决定公开这个秘密的时候经历了多久的挣扎。 很多人对色盲或色弱有着不小的误解,虽然关于色盲症的研究早在两百多年前就已经出现了。 ...

  • 15

    你们做梦都想吃的米其林餐厅,其实全是轮胎厂的阴谋

    欧洲吃货界一直流传着一本“美食圣经”——《米其林指南》。 据说《米其林指南》收录了全欧洲最值得去的餐饮店。 这本书中最高级的便是三颗星的餐厅,而每一颗星的增减都会影响餐厅每月百万欧元的收入。 在业界人士眼里,《米其林指南》是不折不扣的“餐饮界的江湖榜”。 1900年出版第一册《...

评论 抢沙发

  • 昵称 (必填)
  • 邮箱 (必填)
  • 网址

SME 发掘你不知道的科技故事