Search
Question: Can gage not use large gene sets for pathway analysis?
0
gravatar for Guest User
4.5 years ago by
Guest User12k
Guest User12k wrote:
Hi Weijun Luo, Using your two packages of gage and pathview, I have finished pathway analysis of signal transduction and identified two disturbed pathways from 40 pathways. But I met problem with metabolism pathway analysis, because I could not find any disturbed pathway with the metabolism gene sets (82 pathways) from Kegg. More weird, If I used whole gene sets for honeybees from kegg including the met/sig gene sets 122 pathways, the two disturbed signal pathways could not be find too. Some guys told me that this resulted from more gene sets, which increase p. values of each gene set. Is it true, and how can I solve this problem? Thanks a lot!! Richard -- output of sessionInfo(): XML Parsing Error: no element found Location: file:///home/wenfu/CAseqanalysis/ameNA.xml Line Number 2, Column 1: -- Sent via the guest posting facility at bioconductor.org.
ADD COMMENTlink modified 4.5 years ago by Luo Weijun1.4k • written 4.5 years ago by Guest User12k
0
gravatar for Luo Weijun
4.5 years ago by
Luo Weijun1.4k
United States
Luo Weijun1.4k wrote:
Hi Richard, It is normal that you don?t get any significant calls in a pathway analysis (with multiple pathways/tests) because none of the p-values (or q-values) is small enough. Very likely, there is not enough testing power with your data given its sample size, noise level and experiment quality. The adjusted p-values (or q-values) would be different when the total number of tests/pathways changes. With your current dataset, you may do 2 things: -Loosen the selection criteria (q-value cutoff), the option cutoff = 0.1 in sigGeneSet function and q.cutoff = 0.1in gagePipe function can be set to a bigger value, say 0.2 etc. -change the gene set size filter to include more pathways that are actually tested, the argument in gage function is set.size = c(10, 500). You can set it to be set.size = c(10, 2000) or even set.size = c(10, Inf). HTH -------------------------------------------- On Wed, 12/4/13, Richard <richard@> wrote: Subject: Can gage not use large gene sets for pathway analysis? To: bioconductor at r-project.org, mao at life.uiuc.edu Date: Wednesday, December 4, 2013, 10:28 AM Hi Weijun Luo, Using your two packages of gage and pathview, I have finished pathway analysis of signal transduction and identified two disturbed pathways from 40 pathways. But I met problem with metabolism pathway analysis, because I could not find any disturbed pathway with the metabolism gene sets (82 pathways) from Kegg. More weird, If I used whole gene sets for honeybees from kegg including the met/sig gene sets 122 pathways, the two disturbed signal pathways could not be find too. Some guys told me that this resulted from more gene sets, which increase p. values of each gene set. Is it true, and how can I solve this problem? Thanks a lot!! Richard -- output of sessionInfo(): XML Parsing Error: no element found Location: file:///home/wenfu/CAseqanalysis/ameNA.xml Line Number 2, Column 1: -- Sent via the guest posting facility at bioconductor.org.
ADD COMMENTlink written 4.5 years ago by Luo Weijun1.4k
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.2.0
Traffic: 403 users visited in the last hour