[BioC] Gene names in limma Venn Diagrams

James Wettenhall wettenhall at wehi.edu.au
Tue Sep 23 11:48:35 MEST 2003


On Tue, 23 Sep 2003 c.chaffer at pgrad.unimelb.edu.au wrote:
> Hi James,
> How's it going?  I have a question for you!!
> I've been playing with my data in Limma and am wondering if there is a way 
> to get a list of the genes that are represented in the intersection of the 
> venn diagram??

Use classifyTests to get a classification object "clas" which 
you can use in vennCounts.

Then you can use :
clasWithGenes <- cbind(genelist,clas$classification)
to combine the genelist information with the classification. 
Make sure that nrow(genelist)==nrow(clas$classification).  
If not, you might need to run unwrapdups on genelist.

Then to see only the rows of clasWithGenes with genes 
upregulated in both comparisons, use :

The number of rows of this result should be the same as the 
corresponding count in 


More information about the Bioconductor mailing list