我有几个基因特征载体,其中包含了它们所在物种的名称,我制作了一个UpSetR图,显示了基因间共同的物种数量.现在我想做相反的事情:绘制物种间共同基因的数量,但我不知道该怎么做.
我的例子:
gene1 <- c("Panda", "Dog", "Chicken")
gene2 <- c("Human", "Panda", "Dog")
gene3 <- c("Human", "Panda", "Chicken")
...#About 20+ genes with 100+ species each
我希望得到的结果示例:
Panda <- c("gene1", "gene2", "gene3")
Dog <- c("gene1", "gene2")
Human <- c("gene2", "gene3")
Chicken <- c("gene1", "gene3")
...
我知道它在概念上很容易,但后勤更复杂.任何人都可以给我一个线索吗？
谢谢!