Question: Discrepancy in gene set size between gage() and geneData()
gravatar for serpalma.v
4.1 years ago by
serpalma.v30 wrote:

Dear community,

I have run GAGE and Pathview on my microarray data. Everything worked fine.

However, the number of genes assigned to a kegg pathway after using thegage() function is different from the number of genes produced by thegeneData() function, which is the function to use in order to extract the gene IDs assigned to a significantly enriched pathway.

For example, one pathway has181 genes in the variable set.size after using gage(), but when that pathway is used in geneData() with the whole expression set, it returns 196. I tried two more pathways and similar results appeared 34 vs 37 and 15 vs 17.

How can I know exactly which are the genes that are contained in the variable set.size after using gage()?

Many thanks!

ADD COMMENTlink modified 18 months ago by RamRS26k • written 4.1 years ago by serpalma.v30
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 797 users visited in the last hour