R语言中的t检验丨数析学院

2016-09-27 10:47 阅读(471)    评论(0)   

问题

在R中,我们如何检验从总体中抽样得到的两组样本是否有不同的均值,或是通过总体中的某一组样本检验总体均值与某一理论均值间的差异。

指南

示例数据 我们选择内置的sleep数据集作为示例数据。

接下来,我们将sleep数据集处理为宽数据;在之前发布的教程中,我们曾谈到过数据集长宽转换的方法(R语言行动指南-36:长数据与宽数据的转换)。

两组数据的对比:独立双样本t检验

首先假设我们的两组数据完全由独立抽样得来;为了达成这个目的,我们暂且忽略掉ID这个变量的存在。

t.test函数可以直接用于检验像sleep这样的长数据,在下面的例子中,列extra中记录着用于检验的数值,而列group则为样本的分组依据;如果数据没有以数据框形式存储,我们同样可以通过指定两个独立的向量来完成这一操作。

默认地,t.test不会假设样本具有相同的方差,因此该函数默认地调用Welch t检验方法而不是student t检验。可以看到在上文的Welch t检验中,自由度df=17.776,这是由于方法内对于可能存在的非均等方差的调整。如果我们要调用student t检验方法,那么我们需要设置参数var.equal=TRUE。

配对样本t检验

有时需要我们检验的数据可能是同一个实验对象在实验前后的某一关键指标数值,或是相互匹配的两组对象在进行不同实验之后的数值反应,也就是说,我们得到的两组样本数据间具有某种配对关系。此时,我们可以通过配对样本t检验的方式对其进行检验。

同样的,在上文中我们提及了t-test函数可以兼用于一个带有分组变量的数据框或是两个独立的向量。在配对样本中,样本的配对关系取决于其对应的位置。如果我们的数据集为包含分组变量的数据框,那么程序将默认group=1的数据行中的第一行与group=2的数据行中的第一行相互匹配。所以我们需要特别注意数据的排列顺序并确保其中没有缺失值,否则样本间的配对就不得不被打破。在下面的例子中,我们运用group和ID两个变量来确保数据排序的正确。

配对t检验的实质等同于检验每组相互配对的样本数据的差值的总体均值是否为0。(详情可见下文中单样本t检验的内容)

样本与外生总体均值的对比:单样本t检验

假设现在我们的目的为检验列extra数据对应的总体均值是否为0,在下面的例子中,我们暂时忽略了变量group与变量ID。


复制链接,通过电脑学习效果最佳: http://datacademy.io/lesson/64

分享到:

相关推荐

  • 微信图片_20170827200359

    你狠心冲掉的便便,正成为许多肠胃病人的救星

    自从1928年抗生素被发现以来,人类就好比在对抗细菌上获得了超级武器。 遗憾的是,超级武器并没有带来全面的胜利。 人类与细菌旷日持久的军备竞赛催生出了一批更加难缠的超级细菌。 实际上,如果你在医院注射抗生素的话,可能也会被它们盯上。 这是一种叫艰难梭菌*(Clostridium diff...

  • 1

    神奇女侠背后鲜为人知的故事,”神奇博士“居然还发明了测谎仪

    神奇女侠是历史上第一位超级女英雄,与超人、蝙蝠侠并驾齐驱,被称为的“DC漫画三巨头”。 超级英雄都有属于自己的超级武器,但神奇女侠的武器却比较奇特,名叫“真言套索”(Lasso of Truth)。 顾名思义,除了能在战斗中作绳鞭武器外,它最奇妙的设定莫过于只要被其套住,无论是人是神都得吐露真言...

  • NewsImage_Hugh-Herr

    攀岩天才17岁失去双腿,重拾学业自制假肢,依靠科技再攀高峰

    一场意外,有的时候会夺走人的生命,有的时候夺走的是人的自由。 他们失去了双脚,被剥夺了自由行走、奔跑的能力,没得选择地戴上了假肢、坐上了轮椅。 对他们而言,假肢虽然赋予他们重新行走的能力,却也成为他们的一个雷区。 但有一种发明,在生理作用上足以媲美真腿。 这种发明能...

评论 抢沙发

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

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