[BioC] Gene names in limma Venn Diagrams
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