[BioC] extract genes from vennDiagram LIMMA ..How?,

Saurin Jani saurin_jani at yahoo.com
Wed Mar 23 19:48:20 CET 2005


Hi BioC,

For 3 experimental groups in LIMMA, when vennDiagram
is generated based on : 

mask <- c(1,1,1,2,2,3,3,3);
design <- model.matrix(~ -1+factor(mask));
colnames(design) <- c("Group1","Group2","Group3");
ConX<-as.matrix(c("Group1-Group2","Group1-Group3","Group2-Group3"));
fit2 <- contrasts.fit(fit,contrast.matrix);
fit2 <- eBayes(fit2);
ConRow <- nrow(ConX);

#now I am extracting all DE genes from each contrast
and collect it in array, so, later I can draw heatmap
but I have lots of probesets co in that array
for(i in 1:ConRow)
           {
            DEtable <- topTable(fit2,coef = i,number =
numberDE,adjust = "fdr",sort.by = "P");
            DEList[i] <- list(DEtable);
            AllDEGenes <-
c(AllDEGenes,as.character(DEList[[i]][,"ID"]));
            AllDE_FDR <-
c(AllDE_FDR,as.character(DEList[[i]][,"P.Value"]));
           }

For 3 experimental groups, how can someone extract
differentially expressed genes which is common in all
three groups?  

Is it possible to generate vennDiagram() for more than
3 groups..in LIMMA? 


Thank you in advance,
Saurin



More information about the Bioconductor mailing list