绘制venn图的时候,我们数据的ID中有缺失值,也就是不同的ID行数不一样,需要绘制韦恩图,那么空值就读成了NA,需要去除才行;
数据如下,不同比较组的得到的差异基因ID,需要用R绘制韦恩图:
代码及注意事项如下:
library(Vennerable)
#由于行数不一样,设置fill参数;na.strings 设置空值为NA:
All_DEG_venn<-read.table("All_DEG_venn.genes",header=TRUE,comment.char="",sep = "\t",check.names=FALSE,fill=T,na.strings="")
colnames(All_DEG_venn)<-c("DAF2 vs GDAF2","DAF5 vs GDAF5","DAF11 vs GDAF11","DAF16 vs GDAF16")
#去除数据中的NA,以免多计数
All_DEP_ID=apply(All_DEG_venn,2,na.omit)
datav<-Venn(All_DEG_venn)
plot(datav,doWeight=F,type="ellipses")
如果觉得我的文章对您有用,请随意打赏。你的支持将鼓励我继续创作!