两组间差异分析之T检验在R中实现
做统计学分析的其中一个重要目的就是寻找组间差异,在研究中,我们最常关注的问题莫过于处理组与对照组是否存在了显著不同。就两组间比较而言,t检验是常见的分析方法之一。
本文简介怎样在R中进行t检验,以实现两组间差异分析。
本文使用的作图数据的网盘链接(提取码 6268):
https://pan.baidu.com/s/108oBkDTkuJ1HGBPX0h5Osg
文件“alpha.txt”为某16S细菌群落测序所获得的部分alpha多样性指数数据,包含3列信息:sample,样本名称;observed_species和shannon分别为两种类型的alpha多样性指数。文件“group.txt”为各样本分组信息,第一列(sample)为各样本名称;第二列(group)为各样本的分组信息。
接下来,我们期望在R中运行t检验,以查看不同分组间(两两分组之间)的各alpha多样性指数是否存在显著不同。
数据预处理及正态性假设检验
数据预处理
首先将上述两个数据表读入 R 中,并合并在一起。library(reshape2)
#读入文件,合并分组信息,数据重排
alpha group alpha
#我们期望查看 group1 和 group2 的 observed_species 指数是否存在显著差异
#选择要比较的分组
richness_12 richness_12$group head(richness_1