venn diagram怎么做_venn diagram 中文一般_venn diagram英语解释(3)

n12 = 188,

n13 = 191,

n14 = 184,

n15 = 177,

n23 = 194,

n24 = 197,

n25 = 190,

n34 = 190,

n35 = 173,

n45 = 186,

n123 = 112,

n124 = 108,

n125 = 108,

n134 = 111,

n135 = 104,

n145 = 104,

n234 = 111,

n235 = 107,

n245 = 110,

n345 = 100,

n1234 = 61,

n1235 = 60,

n1245 = 59,

n1345 = 58,

n2345 = 57,

n12345 = 31,

category = c("A", "B", "C", "D", "E"),

fill = c("dodgerblue", "goldenrod1", "darkorange1", "seagreen3", "orchid3"),

cat.col = c("dodgerblue", "goldenrod1", "darkorange1", "seagreen3", "orchid3"),

cat.cex = 2,

margin = 0.05,

cex = c(1.5, 1.5, 1.5, 1.5, 1.5, 1, 0.8, 1, 0.8, 1, 0.8, 1, 0.8, 1, 0.8,

1, 0.55, 1, 0.55, 1, 0.55, 1, 0.55, 1, 0.55, 1, 1, 1, 1, 1, 1.5),

ind = TRUE

);

library("multtest"); library("annaffy"); library("hgu95av2.db");library("ALL")

library("genefilter")

data(ALL, package = "ALL")

ALLB <- ALL[,ALL$BT %in% c("B1","B2","B3","B4")]

f1 <- function(x) (shapiro.test(x)$p.value > 0.05)

sel1 <- genefilter(exprs(ALL[,ALLB$BT=="B1"]), filterfun(f1))

sel2 <- genefilter(exprs(ALL[,ALLB$BT=="B2"]), filterfun(f1))

sel3 <- genefilter(exprs(ALL[,ALLB$BT=="B3"]), filterfun(f1))

sel4 <- genefilter(exprs(ALL[,ALLB$BT=="B4"]), filterfun(f1))

selected <- sel1 & sel2 & sel3 & sel4

library(limma)

/# Three-set Venn diagram

x <- matrix(as.integer(c(sel2,sel3,sel4)),ncol = 3,byrow=FALSE)

colnames(x) <- c("sel2","sel3","sel4")

vc <- vennCounts(x, include="both")

vennDiagram(vc)

/# Four-set Venn diagram

x <- matrix(as.integer(c(sel1,sel2,sel3,sel4)),ncol = 4,byrow=FALSE)

colnames(x) <- c("sel1","sel2","sel3","sel4")

vc <- vennCounts(x, include="both")

vennDiagram(vc)

本文来自电脑杂谈，转载请注明本文网址：
http://www.pc-fly.com/a/jisuanjixue/article-48795-3.html