mat <- matrix(rnorm(1e6), ncol = 20)
design <- model.matrix(~factor(1:4, each=5))
colnames(design) <- LETTERS[1:4]
contrast <- matrix(c(1,-1,0,0,1,0,-1,0,1,0,0,-1,0,1,-1,0,0,1,0,-1),
ncol = 5)
colnames(contrast) <- paste(LETTERS[c(1,1,1,2,2)],
LETTERS[c(2,3,4,3,4)], sep = " vs ")
fit <- lmFit(mat, design)
fit2 <- contrasts.fit(fit, contrast)
fit2 <- eBayes(fit2)
## two Venn diagrams - a 3-way Venn with the first three contrasts
## and a 2-way Venn with the last two contrasts
collist <- list(1:3,4:5)
venn <- makeVenn(fit2, contrast, design, eset, collist = collist)
vennreport <- vennPage(venn, "index.html", "Venn diagrams")
indexPage <- HTMLReport("index", "My results", reportDirectory =
".", baseUrl = ".")
publish(vennreport)
finish(indexPage)
Run the code above in your browser using DataLab