#安装UpSetR
install.packages("UpSetR")
library(UpSetR)
#install.packages("UpSetR")
library(UpSetR)
library(Cairo)
# 从CSV文件中读取数据
setwd("C:/Users/fordata/Desktop/研究生/第二个想法(16s肠型+宏基因组功能)/第二篇病毒组/result/otherDB")
data <- read.csv("3462+11.csv")
data <- data[, colSums(is.na(data)) == 0]
# 创建一个新的PDF文件
CairoPDF("your_file.pdf", width = 7, height = 5)
upset(data,
sets = c("URPC", "RVD", "IMG_VR4",
"GVD", "IGVD", "NEXT_PILOT"),
number.angles = 0,
point.size = 2,
line.size = 1,
mainbar.y.label = "vOTU Counts",
sets.x.label = "vOTUs Per Database",
text.scale = c(1.1, 1.1, 1, 1, 1, 1),
mb.ratio = c(0.7, 0.30),
nintersects = NA,
sets.bar.color = "#456990",
main.bar.color = "#EA8379"
)
dev.off()
R 数据可视化 —— 集合可视化 UpSetR_r语言如何提取upset图中各个子集的对象?-CSDN博客
超详细的教程可见这个!!