嗨,我是R的新手,我必须用它来制作维恩图。我已经谷歌搜索了一段时间,我能找到的所有例子都处理二进制变量。但是我有2个列表(实际上是2个csv文件)。列表中的项目只是字符串,如PSF113_xxxx。我必须比较它们,看看每个共享的内容和共享内容。我如何在R中制作一个维恩图?
此外,文件中的文件数量不同,一个文件的数量略多于另一个文件,这意味着cbind函数会返回错误。
到目前为止,我已经想出了这个,但这只是给我一个带有一个名为group 1的圆圈的图像,其中1为内部,0为外部。
matLis <- list(matrix(A), matrix(B))
n <- max(sapply(matLis, nrow))
do.call(cbind, lapply(matLis, function (x)
rbind(x, matrix(, n-nrow(x), ncol(x)))))
x = vennCounts(n)
vennDiagram(x)
这是数据的一个例子
2 PSF113_0018
3 PSF113_0079
4 PSF113_0079a
5 PSF113_0079b
左边的编号不是我做的任何事情,它补充说当我从excel将文件导入R时
head(A)
> head(A)
V1
1 PSF113_0016a
2 PSF113_0018
3 PSF113_0079
4 PSF113_0079a
5 PSF113_0079b
6 PSF113_0079c
> head(b,10)
V1
1 PSF113_0016a
2 PSF113_0021
3 PSF113_0048
4 PSF113_0079
5 PSF113_0079a
6 PSF113_0079b
7 PSF113_0079c
8 PSF113_0295
9 PSF113_0324a
10 PSF113_0324b
答案 0 :(得分:2)
您的代码仍然不具备可重现性,因为您还没有定义A或B.这里是venneuler
包中维恩图的指南,因为我发现它更灵活。
List1 <- c("apple", "apple", "orange", "kiwi", "cherry", "peach")
List2 <- c("apple", "orange", "cherry", "tomatoe", "pear", "plum", "plum")
Lists <- list(List1, List2) #put the word vectors into a list to supply lapply
items <- sort(unique(unlist(Lists))) #put in alphabetical order
MAT <- matrix(rep(0, length(items)*length(Lists)), ncol=2) #make a matrix of 0's
colnames(MAT) <- paste0("List", 1:2)
rownames(MAT) <- items
lapply(seq_along(Lists), function(i) { #fill the matrix
MAT[items %in% Lists[[i]], i] <<- table(Lists[[i]])
})
MAT #look at the results
library(venneuler)
v <- venneuler(MAT)
plot(v)
编辑:头脑非常有帮助,因为它为我们提供了一些工作。试试这种方法:
#For reproducibility (skip this and read in the csv files)
A <- structure(list(V1 = structure(1:6, .Label = c("PSF113_0016a",
"PSF113_0018", "PSF113_0079", "PSF113_0079a", "PSF113_0079b",
"PSF113_0079c"), class = "factor")), .Names = "V1",
class = "data.frame", row.names = c("1",
"2", "3", "4", "5", "6"))
B <- structure(list(V1 = structure(1:10, .Label = c("PSF113_0016a",
"PSF113_0021", "PSF113_0048", "PSF113_0079", "PSF113_0079a",
"PSF113_0079b", "PSF113_0079c", "PSF113_0295", "PSF113_0324a",
"PSF113_0324b"), class = "factor")), .Names = "V1",
class = "data.frame", row.names = c("1",
"2", "3", "4", "5", "6", "7", "8", "9", "10"))
#after reading in the csv files start here
Lists <- list(A, B) #put the word vectors into a list to supply lapply
Lists <- lapply(Lists, function(x) as.character(unlist(x)))
items <- sort(unique(unlist(Lists))) #put in alphabetical order
MAT <- matrix(rep(0, length(items)*length(Lists)), ncol=2) #make a matrix of 0's
colnames(MAT) <- paste0("List", 1:2)
rownames(MAT) <- items
lapply(seq_along(Lists), function(i) { #fill the matrix
MAT[items %in% Lists[[i]], i] <<- table(Lists[[i]])
})
MAT #look at the results
library(venneuler)
v <- venneuler(MAT)
plot(v)
这种方法的不同之处在于我没有列出两个数据帧(如果它们是数据帧),然后将它们转换为字符向量。我认为这应该有用。