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

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

分享到:

相关推荐

  • 微信截图_20170721164152

    生育专家骗腹生子12年造75人,开创“共享云爹”斩获搞笑诺奖

    人类对传宗接代的欲望大都高涨,因为生育闹出的荒唐事也着实不少。 1926年,加拿大一位富豪查尔斯·万斯·米勒与世长辞,但他仍为加拿大人了举办了一场盛大的竞赛。 他让律师保管他留下的900万加元,在10年后将钱交给多伦多市内这10年生孩子最多的家庭。 协约书 刚巧的是,这10年间经济大萧...

  • ad1ea50f4bfbfbeddb2b859f7ef0f736aec31f15

    千年法老不懂诅咒,掘墓者却仍相继离奇死去?

    “谁要是干扰了法老的安宁,死亡就会降临到他的头上。” 1923年4月23日,法老图坦卡蒙陵墓挖掘工作的资助人卡纳冯伯爵在旅馆中神秘去世。 英国的各大报刊相继报道了这起离奇的事件。 报道中称,卡纳冯伯爵在法老陵墓被发掘不到半年内,左脸被奇怪的蚊子叮咬。 他在刮胡子的时候不小心刮...

  • PaulKarason

    无毒无害的贵金属,做成保健品却把人吃成了“蓝血贵族”

    作为贵金属界的扛把子,金和银两者通常都是以装饰、货币、等价物的形式示人的。 但随着人类社会的发展,银和金的命运却变得有些不太一样。 若抛开价值属性只谈实用性,白银比黄金要强上不少。 人类社会由石器时代进入金属时代后,就有人发现银质器皿中的水要比其他容器中的要清洁。 ...

评论 抢沙发

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

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