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

2016-09-27 10:47 阅读(364)    评论(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

分享到:

相关推荐

  • memorable-moments-1896-600x400

    百年现代奥运史,千年兴奋剂滥用史

    自古以来,人类靠着一个又一个的科技突破不断挑战新的极限。 从莱特兄弟的飞行者1号到美国宇航局研制时速上万公里的X-43A。 从古老的摩尔斯码电报机到全球互联的移动智能手机。 我们享受着这样的进程,同时也积极努力地推进技术的发展。 然而,有一个领域足以代表人类科技的最尖端水...

  • 活动行封面1080

    广州首届国际创客节开幕,大众直呼“过瘾”!

    5月20日,广州首届国际创客节——2017天英汇•IEF国际创客节在广州开幕。IEF国际创客节是每年一度针对国际创客圈展开的大型盛会,本届创客节是由SME、INNOTalk、中大创新谷、云珠沙龙和G-SICE校园直通车联合主办,吸引了来自国内外的科技爱好者、国际创客、重点科技企业代表等数千人集聚广州,共同见证了广州...

  • 0000

    “爱因斯坦二世”8天编出一篇论文,顶级期刊畅通无阻,诺奖得主全力支持

    这世上什么都能骗,骗财、骗色、骗名声。 出门在外、行骗为本的人都要讲究个唬得住人。 好比遇到要招魂抓鬼的老人,先给他表演两手赤手探油锅、招魂大法。 样子上过得去,再谈条件也就不难。 这些常见的骗子都是骗知识水平不高的人,聪明点的就不容易被唬住。 但曾经声名鹊起的扬·...

评论 抢沙发

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

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