我有三个 ID 列表。
我想比较 3 个列表,并绘制维恩图。在获得的维恩图中,我将在交叉点中显示的不是数字而是 ID。
我需要在 R 中做到这一点,但我真的不知道怎么做。
你可以帮帮我吗?
那是我的代码。它有效,但只显示数字,我会将“术语”显示为交叉点
set1 <- unique(goterm1)
set2 <- unique(goterm2)
set3 <- unique(goterm3)
require(limma)
Diagram <- function(set1, set2, set3, names)
{
stopifnot( length(names) == 3)
# Form universe as union of all three sets
universe <- sort( unique( c(set1, set2, set3) ) )
Counts <- matrix(0, nrow=length(universe), ncol=3)
colnames(Counts) <- names
for (i in 1:length(universe))
{
Counts[i,1] <- universe[i] %in% set1
Counts[i,2] <- universe[i] %in% set2
Counts[i,3] <- universe[i] %in% set3
}
vennDiagram( vennCounts(Counts) )}
Diagram(set1, set2, set3, c("ORG1", "ORG2", "ORG3"))
Venn
最佳答案
您也可以使用 limma 完成这项壮举。请参阅下面的示例。
这个想法基本上与您发布的代码完全相同,但它没有被包装到一个函数中(因此可能更容易调试)。
你让它与下面的代码一起工作吗?如果没有,请发布您收到的可能的错误消息和警告。
# Load the library
library(limma)
# Generate example data
set1<-letters[1:5]
set2<-letters[4:8]
set3<-letters[5:9]
# What are the possible letters in the universe?
universe <- sort(unique(c(set1, set2, set3)))
# Generate a matrix, with the sets in columns and possible letters on rows
Counts <- matrix(0, nrow=length(universe), ncol=3)
# Populate the said matrix
for (i in 1:length(universe)) {
Counts[i,1] <- universe[i] %in% set1
Counts[i,2] <- universe[i] %in% set2
Counts[i,3] <- universe[i] %in% set3
}
# Name the columns with the sample names
colnames(Counts) <- c("set1","set2","set3")
# Specify the colors for the sets
cols<-c("Red", "Green", "Blue")
vennDiagram(vennCounts(Counts), circle.col=cols)
该代码应给出类似于以下内容的图:
关于r - 如何使用R绘制维恩图,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/19467073/