Question: Discrepancy in gene set size between gage() and geneData()
0
gravatar for serpalma.v
3.1 years ago by
serpalma.v20
Germany
serpalma.v20 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 6 months ago by RamRS20k • written 3.1 years ago by serpalma.v20
Please log in to add an answer.

Help
Access

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