VennDiagram.R

From Organic Design wiki
Revision as of 01:58, 28 May 2007 by Sven (talk | contribs)
# This illustrates a problem with missing values
library(limma)
p <- 100
n <- 3

x <- matrix(0, nr=p,nc=n)
set.seed(1)
x[,1] <- rbinom(p, 1, 1/10)
x[,2] <- rbinom(p, 1, 1/10)
x[,3] <- rbinom(p, 1, 1/10)
x

vennCounts(x)[7,]
vennDiagram(x)

x[4,3] <- NA

vennCounts(x)[7,]
vennDiagram(x)