【发布时间】:2022-01-09 22:01:32
【问题描述】:
我有一份测序基因列表。这个基因列表被注释为一个 GEN-ID。下面的链接显示了一个列表的示例。
在此列表中,数字表示在样本中找到了 gen。 NA 表示找不到gen。
我想以两个这样的列表为例,用它们做一个维恩图。我为此使用的代码如下:
listA <- read.csv("young.csv", header = FALSE)
A <- listA
A
listB <- read.csv("old.csv", header = FALSE)
B <- listB
B
length(A$V2)
length(B$V2)
A[is.na(A)] <- ""
B[is.na(B)] <- ""
library(VennDiagram)
xx.1 <- venn.diagram(list("young_control" =A$V2, "old_control" = B$V2), fill = c("yellow","cyan"), cex
=4.0, filename = "venn_excersice.png")
通过这个,我生成了一个维恩图,显示了仅在 listA 中、仅在 listB 中或同时在 listA 和 listB 中的基因。
我现在的问题是: 如果我想从生成的 Venn 的每个区域获取 Gen-ID 列表,我该怎么做? 我尝试了 attr、output、intersect 和其他函数,但不知何故它不起作用,我不知道为什么会这样。
此外,我想生成一个交互式维恩图。您可以在其中单击该区域并立即获得该区域中 Gen-ID 的可视化表示。
如果有人对我可以尝试达到这个目标有什么建议,我会非常高兴!
列表格式:
列表 A
ID.........Y
G-1........1
G-2.......不适用
G-3........3
G-4........4
列表 B
ID.........O
G-1........1
G-2....... 2
G-3........3
G-4........不适用
【问题讨论】:
-
您能与我们分享您的
listA和listB的数据吗?您可以使用dput(listA)和dput(listB) -
我不确定如何在堆栈溢出中使用 dput。我希望我对列表的表示以及您可以在超链接上找到的图片对您有所帮助。
-
您在 R 中使用
dput(),然后将输出复制并粘贴到您的问题中 -
请提供足够的代码,以便其他人更好地理解或重现问题。
标签: r list bioinformatics venn-diagram genetics