【R语言】avereps函数处理相同基因表达值

前面给大家介绍了，如果同一个基因有多行表达值，该如何处理。

前面给大家介绍了如何通过R的aggregate函数来实现对相同的基因名按列取平均。今天小编给大家介绍另外一种方法，使用limma包中的avereps函数来处理。大家对limma这个包应该都不陌生。它是目前主流的三种差异表达分析包（limma,edgeR和DEseq2）中的一个，也是应用比较广泛的一个。

前面小编也给大家介绍了

☞R代码TCGA差异表达分析

☞零代码TCGA差异表达分析

☞GEO芯片数据差异表达分析

今天我们来聊聊如何使用limma包中的avereps，对相同基因的表达值取平均。

#加载limma包library(limma)#生成一个表达矩阵df=matrix(1:20,nrow=4)#设置行名，基因名字rownames(df)=paste0("gene",c(1,1,2,3))#设置列名，样本名字colnames(df)=paste0("sample",1:5)#对相同的基因的表达取平均result=avereps(df)result

原始矩阵如下

通过avereps取完平均之后的矩阵如下

我们也可以使用☞ 表达谱数据中相同基因如何处理中讲到的aggregate函数来处理这个矩阵得到相同的结果。

#使用aggregate函数#将基因名字添加到数据框中df1=data.frame(gene=rownames(df),df)#使用aggregate函数取平均result1=aggregate(.~gene,df1,mean)#设置行名，并删掉gene名字这一列rownames(result1)=result1$generesult2=result1[,-1]result2

可以看到，跟avereps方法得到的结果是一样的，殊途同归。

为了方便大家交流学习，共同进步，我特地创建了微信交流群

本站仅提供存储服务，所有内容均由用户发布，如发现有害或侵权内容，请点击举报。