我有一组微阵列数据,我想根据附图(用PowerPoint完成)表示它。我尝试了各种R软件包(VennDiagram、venneuler、limma)。但是,我无法按比例绘制数据或显示值。
有人知道怎么做吗?
代码玩可以在下面找到!
谢谢大家!
# you might need these:
source("http://www.bioconductor.org/biocLite.R")
biocLite("limma")
library(VennDiagram)
library(limma)
library(venneuler)
# Trial 1: Kind of proportional but difficult to customize for labels and other stuff
su <- venneuler(c(A=162, B=104, C=86, "A&B"=206, "A&C"=112, "B&C"=90 ,"A&B&C"=2433))
plot(su)
# Trial 2: looks nice but not proportional
hw <- c(F,F,F,F,T,T,T,T)
hm <- c(F,F,T,T,F,F,T,T)
hr <- c(F,T,F,T,F,T,F,T)
c4 <- cbind(hw,hm,hr)
e <- vennCounts(c4)
e[1:8,4] <- c(3193,86,104,90,162,112,206,2433)
vennDiagram(e)
2条答案
按热度按时间erhoui1w1#
我推荐VennDiagram软件包:http://cran.r-project.org/web/packages/VennDiagram/VennDiagram.pdf
在pake19上,你会发现10个非常好的例子(包括高级和简化的),到目前为止,我还没有发现任何我需要它做的事情。
作为对其他人的一般帮助谁结束在这里关于这个问题(不直接相关的这个问题,但相关),我会建议寻找您的答案使用此搜索
oogrdqng2#
在我的
nVennR
包中,可以通过两种方式来完成,如vignette所示。包可以表示更大数量的集合。