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

520 阅读2分钟

单样本t检验是用来确定一个群体的平均值是否等于某个值的。

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

t.test(data, mu=10)

下面的例子展示了如何在实践中使用这种语法。

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

假设一位植物学家想知道某一种植物的平均高度是否等于15英寸。

她收集了12种植物的简单随机样本,并记录了每一种植物的高度(英寸)。

她可以使用下面的代码在R中进行单样本t检验,以确定该种植物的平均高度是否真的等于15英寸:

#create vector to hold plant heights
my_data <- c(14, 14, 16, 13, 12, 17, 15, 14, 15, 13, 15, 14)

#perform one sample t-test
t.test(my_data, mu=15)

	One Sample t-test

data:  my_data
t = -1.6848, df = 11, p-value = 0.1201
alternative hypothesis: true mean is not equal to 15
95 percent confidence interval:
 13.46244 15.20423
sample estimates:
mean of x 
 14.33333 

下面是如何解释输出中的每个值:

data:t检验中使用的向量的名称。在这个例子中,我们使用了my_data

t:t检验统计量,计算为(x-μ)/(s√n)=(14.333-15)/(1.370689/√12)=-1.6848

df:自由度,计算为n-1=12-1=11

p值:与-1.6848和11个自由度的t检验统计量相对应的双尾p值。在这种情况下,p=0.1201

95%的置信区间:真实人口平均数的95%置信区间,计算为**[13.46244, 15.20423]**。

这个单样本t检验的无效假设和备选假设如下。

H0: µ = 15 (该种植物的平均高度为15英寸)

HA : µ ≠15 (平均高度不是 15英寸)

因为我们测试的P值 (0.1201) 大于0.05,所以我们不能拒绝测试的无效假设。

这意味着我们没有足够的证据说这个特定种类的植物的平均高度与15英寸不同。