6

我有一组微阵列数据,我想按照附图(用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)

在此处输入图像描述

4

2 回答 2

0

我会推荐 VennDiagram 包: http ://cran.r-project.org/web/packages/VennDiagram/VennDiagram.pdf

在 pake 19 中,您会发现 10 个非常好的示例(包括高级示例和简化示例)。到目前为止,我还没有发现我需要它做的任何事情。

作为对最终来到这里的其他人的一般帮助(与此问题不直接相关,但相关),我建议您使用此搜索来寻找您的答案

于 2014-03-21T13:47:35.247 回答
0

使用我的nVennR包,这可以通过两种方式完成,如小插图所示。在这个例子中,

> library(nVennR)
> myV <- createVennObj(nSets = 3, sNames = c('A', 'B', 'C'), sSizes = c(0, 86, 104, 90, 162, 112, 206, 2433))
> myV <- plotVenn(nVennObj = myV)

三组,比例

包可以代表更大数量的集合。

于 2019-05-24T11:19:29.257 回答