library (VennDiagram)
install.packages("VennDiagram")
lmff= read.csv("lfmm.csv",header = T,stringsAsFactors = F)
dim(lmff)
objects(lmff)
Fst= read.csv("Fst.csv",header = T,stringsAsFactors = F)
dim(Fst)
objects(Fst)
XPEHH= read.csv("XPEHH.csv",header = T,stringsAsFactors = F)
dim(XPEHH)
objects(XPEHH)
SAM= read.csv("SAM.csv",header = T,stringsAsFactors = F)
dim(SAM)
X=list(lmff=lmff$lfmm,
Fst=Fst$Fst,
XPEHH=XPEHH$XPEHH,
SAM=SAM$SAM)
venn.diagram(X,filename = "ven1.tiff",fill=c("red","blue","green","yellow"))
##提取信息
inter <- get.venn.partitions(X)
for (i in 1:nrow(inter)) inter[i,'values'] <- paste(inter[[i,'..values..']], collapse = ', ')
write.table(inter[-c(5, 6)], 'venn4_inter.txt', row.names = FALSE, sep = '\t', quote = FALSE)
lfmm1 = read.csv("lfmm-snp.csv",header = T,stringsAsFactors = F)
dim(lfmm1)
objects(lfmm1)
Fst1= read.csv("Fst-snp.csv",header = T,stringsAsFactors = F)
dim(Fst1)
objects(Fst1)
XPEHH1= read.csv("XPEHH-snp.csv",header = T,stringsAsFactors = F)
dim(XPEHH1)
objects(XPEHH1)
SAM1= read.csv("SAM-snp.csv",header = T,stringsAsFactors = F)
dim(SAM1)
objects(SAM1)
X1=list(lmff1=lfmm1$lmff1,
Fst1=Fst1$Fst1,
XPEHH1=XPEHH1$XPEHH1,
SAM1=SAM1$SAM1)
venn.diagram(X1, filename = "ven2.tiff",fill=c("red","blue","green","yellow"))
inter <- get.venn.partitions(X1)
for (i in 1:nrow(inter)) inter[i,'values'] <- paste(inter[[i,'..values..']], collapse = ', ')
write.table(inter[-c(5, 6)], 'venn5_inter.txt', row.names = FALSE, sep = '\t', quote = FALSE)
维恩图绘制
最后编辑于 ：
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着，像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上，一...
- 文/苍兰香墨 我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
推荐阅读更多精彩内容
- 生信分析miRNA靶点预测R语言韦恩图 制作的代码 install.packages("VennDiagram")...
- 笔者最近几年家里催着结婚，虽然这几年各大亲网站红娘、实体店、家人亲戚朋友也介绍了不少妹纸，各种类型都不少，但由于笔...
