当前位置:网站首页>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)
边栏推荐
- Remote connection to MySQL under windows and Linux system
- PMP personal sprint preparation experience
- Mmsegmentation series training and reasoning their own data set (3)
- Is bone conduction earphone better than traditional earphones? The sound production principle of bone conduction earphones is popular science
- Xiaomi, a young engineer, was just going to make soy sauce
- Redis cluster
- 2022-2028 global encryption software industry research and trend analysis report
- Golang configure export goprivate to pull private library code
- 3124. Word list
- 只需简单几步 - 开始玩耍微信小程序
猜你喜欢
Verilog 状态机
Formatting logic of SAP ui5 currency amount display
The capacity is upgraded again, and the new 256gb large capacity specification of Lexar rexa 2000x memory card is added
ZABBIX API creates hosts in batches according to the host information in Excel files
创业了...
Golang configure export goprivate to pull private library code
Pychart creates new projects & loads faster & fonts larger & changes appearance
Common means of modeling: aggregation
数据传输中的成帧
Design details of SAP e-commerce cloud footernavigationcomponent
随机推荐
Pychart creates new projects & loads faster & fonts larger & changes appearance
MVVM and MVC
The video number will not be allowed to be put on the shelves of "0 yuan goods" in the live broadcasting room?
Use usedeferredvalue for asynchronous rendering
[JSON] gson use and step on the pit
C reflection practice
3124. Word list
Mongodb base de données non relationnelle
[staff] the direction of the symbol stem and the connecting line (the symbol stem faces | the symbol stem below the third line faces upward | the symbol stem above the third line faces downward | the
QT implementation interface jump
批量检测url是否存在cdn—高准确率
STM32__ 05 - PWM controlled DC motor
QT uses sqllite
2022-2028 global nano abrasive industry research and trend analysis report
Mongodb non relational database
实现一个自定义布局的扫码功能
How does proxy IP participate in the direct battle between web crawlers and anti crawlers
2022 hoisting machinery command examination paper and summary of hoisting machinery command examination
Share the basic knowledge of a common Hongmeng application
2022 safety officer-c certificate examination questions and mock examination