在今天的内容,我将向大家介绍如何在R语言里进行t检验分析,暂不涉及相关统计理论。
R语言提供t.test()函数可以进行各种各样的t检验。与其他统计包不同的是,t.test()函数默认数据是异方差的,并采用Welch方法矫正自由度,关于异方差和Welch方法我会在后续的内容中详细介绍,大家先有一个印象即可。
在这里我将利用R里内置的鸢尾花数据集(iris)向大家展示如何进行t检验,这里iris数据集是由150朵鸢尾花的花瓣长度、花瓣宽度、花萼长度、花萼宽度以及鸢尾花种类组成。
setosa <- iris[which(iris$Species=='setosa'),] #提取setosa类的鸢尾花
versicolor <- iris[which(iris$Species=='versicolor'),] #提取versicolor类的鸢尾花
mydata <- rbind(setosa,versicolor) #按行合并数据集
# 独立双样本t检验
t.test(y~x) # y是一个数值型向量,而x是一个二分类变量
t.test(mydata$Sepal.Length~mydata$Species) #检验不同鸢尾花花萼长度差异
从这里我们可以看出两种不同的鸢尾花(setosa和versicolor&#x