当前位置:网站首页>Intersection of Venn graph
Intersection of Venn graph
2022-07-02 03:06:00 【youngleeyoung】
# intersection
my_disease_markers_for_AM3=openxlsx::read.xlsx("G:/silicosis/sicosis/Supplementary Table.xlsx",
sheet = "Sheet1")
head(my_disease_markers_for_AM3)
library(dplyr)
disease_markers_for_AM3=list()
for (eachcolumn in colnames(my_disease_markers_for_AM3)) {
disease_markers_for_AM3[[paste0(eachcolumn)]]=my_disease_markers_for_AM3[,eachcolumn] %>%
na.omit() %>% unique()
}
head(disease_markers_for_AM3)
library(VennDiagram)#https://www.jianshu.com/p/b5a4c40c3a33
venn_list <- list(AM3_Silicosis = disease_markers_for_AM3[[1]],
MoAM_GSE184854 = disease_markers_for_AM3[[2]])
venn.diagram(venn_list, filename = 'venn2_for_humanPF-Silicosi-Bleomycin.png', imagetype = 'png',
fill = c('red', 'blue'), alpha = 0.50, cat.col = rep('black', 2),
col = 'black', cex = 1.3, fontfamily = 'serif',
cat.cex = 1.3, cat.fontfamily = 'serif')
getwd()
# Continue with the above 4 Groups as an example , Intersecting elements between groups get
inter <- get.venn.partitions(venn_list)
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)
for (i in 1:nrow(inter)) inter[i,'values'] <- paste(inter[[i,'..values..']], collapse = ', ')
openxlsx::write.xlsx(inter[-c(5, 6)], 'venn4_inter.xlsx', row.names = FALSE, sep = ',', quote = FALSE)
边栏推荐
- [JSON] gson use and step on the pit
- 【JVM】创建对象的流程详解
- [JS reverse series] analysis of a customs publicity platform
- 2022-2028 global deep sea generator controller industry research and trend analysis report
- GB/T-2423.xx 环境试验文件,整理包括了最新的文件里面
- Actual battle of financial risk control - under Feature Engineering
- Calculation of page table size of level 2, level 3 and level 4 in protection mode (4k=4*2^10)
- JS <2>
- Set status bar color
- Find duplicates [Abstract binary / fast and slow pointer / binary enumeration]
猜你喜欢
随机推荐
What are the common proxy servers and what are the differences?
表单自定义校验规则
Verilog 时序控制
Framing in data transmission
V-model of custom component
verilog 并行块实现
Leetcode question brushing (10) - sequential question brushing 46 to 50
What is the principle of bone conduction earphones and who is suitable for bone conduction earphones
Just a few simple steps - start playing wechat applet
在QML中加载不同字体
Mongodb non relational database
2022-2028 global military computer industry research and trend analysis report
tarjan2
Golang configure export goprivate to pull private library code
Gradle 笔记
Competition and adventure burr
Baohong industry | four basic knowledge necessary for personal finance
Batch detect whether there is CDN in URL - high accuracy
QT uses sqllite
What kind of good and cost-effective Bluetooth sports headset to buy