如何在R中进行双样本T检验

587 阅读2分钟

双样本t检验是用来检验两个种群的平均值是否相等的。

你可以使用下面的基本语法在R中进行双样本t检验:

t.test(group1, group2, var.equal=TRUE) 

注意:通过指定var.equal=TRUE,我们告诉R假设两个样本之间的变异是相等的。

如果你不想做这个假设,只需省略这个参数,R将转而执行韦尔奇的t检验,它不假设样本之间的变异相等。

下面的例子展示了如何在R中实际执行两个样本的t检验。

例子:R语言中的双样本T检验

假设我们想知道两种不同种类的植物是否具有相同的平均高度。

为了测试这一点,我们从每个物种中收集了一个简单的随机样本,即12株植物。

下面的代码显示了如何在R中进行双样本T检验,以确定两个物种之间的平均高度是否相等:

#create vectors to hold plant heights from each sample
group1 <- c(8, 8, 9, 9, 9, 11, 12, 13, 13, 14, 15, 19)
group2 <- c(11, 12, 13, 13, 14, 14, 14, 15, 16, 18, 18, 19)


下面是如何解释检验的结果。

数据:包含样本数据的向量的名称。

t:t检验统计量,在本例中,它是-2.5505

df:自由度,计算为n1 + n2 - 2 = 12 + 12 - 2 =22

p值:与t检验统计量为-2.5505和df=22相对应的p值。结果p值为0.01823。我们可以通过使用T分数到P值计算器来确认这个值。

95%的置信区间。两组之间的真实差异的95%置信区间。结果是**[-5.59, -.576]**。

样本估计值:每组的样本平均数。在这种情况下,第一组的样本平均数是11.667,第二组的样本平均数是14.75

这个特殊的双样本t检验的无效假设和备选假设如下:

H0: µ1 = µ2 (两个群体的平均值相等)

HA :µ1≠µ2(两个群体的平均值相等)

因为我们检验的P值 (0.01823) 小于0.05,所以我们拒绝无效假设。

这意味着我们有足够的证据得出结论,两个物种之间的平均植物高度不相等。

技术说明

R中的t.test()函数使用以下语法:

t.test(x, y, alternative="two.sided", mu=0, paired=FALSE, var.equal=FALSE, conf.level=0.95)

其中:

  • x,y:包含数据的两个向量的名称。
  • 备选假设:备选假设。选项包括 "双面"、"更小 "或 "更大"。
  • mu:假设为平均值的真实差异的值。
  • 成对的:是否使用成对的t检验。
  • var.equal:两组之间的变异量是否相等。
  • conf.level:用于测试的置信度。

当你进行自己的t检验时,可以随意改变这些参数,这取决于你要进行的特定检验。