SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Problems running pBWA swNGS Bioinformatics 3 04-29-2014 05:03 AM
Problems running Apollo in Windows7 gwilymh Bioinformatics 0 07-29-2013 10:24 AM
Problems running ExomeCNV lkadalayil Bioinformatics 2 07-23-2013 12:06 AM
Gage-b pari_89 Bioinformatics 1 07-06-2013 03:43 AM
Problems with running Quake vivienne_lovely Bioinformatics 1 01-18-2012 08:32 AM

Reply
 
Thread Tools
Old 10-20-2013, 06:34 AM   #1
sindrle
Senior Member
 
Location: Norway

Join Date: Aug 2013
Posts: 266
Default Problems running GAGE

Hi!
Here is my code, my input is from DESeq2 results.

deseq2.fc=timepointsResults$log2FoldChange;
names(deseq2.fc)=rownames(timepointsResults);
exp.fc=deseq2.fc;
out.suffix="timepointsResultsDESeq2"
require(gage);
data(kegg.gs);
dir.create("/Volumes/timemachine/HTseq_DEseq2/1_Results/GAGE_upregulated_timepointsResults");
setwd("/Volumes/timemachine/HTseq_DEseq2/1_Results/GAGE_upregulated_timepointsResults");
fc.kegg.p <- gage(exp.fc, gsets = kegg.gs, ref = NULL, samp = NULL, same.dir = TRUE) ;
sel <- fc.kegg.p$greater[, "q.val"] > 0.1 & !is.na(fc.kegg.p$greater[, "q.val"]) ;
path.ids <- rownames(fc.kegg.p$greater)[sel] ;
path.ids2 <- substr(path.ids, 1, 8) ;
require(pathview) ;
pv.out.list <- sapply(path.ids2[1:3], function(pid) pathview(gene.data = exp.fc, pathway.id = pid, species = "hsa", out.suffix=out.suffix))

I just want to look for the top 3 upregulated pathways, but I only get an "404 not found" error.

Anyone understand this?
sindrle is offline   Reply With Quote
Old 10-20-2013, 07:28 AM   #2
sindrle
Senior Member
 
Location: Norway

Join Date: Aug 2013
Posts: 266
Default

This is just wired..
I updated from R 3.0.1 to newest and installed a new package "httr" for GAGE.

And it works.....

Now I just have to understand how a "significantly upregulated pathway" is calculated as significant...

Last edited by sindrle; 10-20-2013 at 12:33 PM.
sindrle is offline   Reply With Quote
Old 10-20-2013, 06:42 PM   #3
bigmw
Senior Member
 
Location: US

Join Date: Aug 2013
Posts: 123
Default

It is always important to update R/Bioconductor to use the latest features of gage and pathview pacakges. In addition, if you installed the packages under current release versions of R(3.0.2)/ Bioconductor(2.13), you don’t really need to worry about any dependencies like httr.
I guess you are using the joint workflow described in the vignette “RNA-Seq Data Pathway and Gene-set Analysis Workflows”. Notice that the native GAGE/Pathview workflow may give you more sensitive results. You may check the main vignette “Generally Applicable Gene-set/Pathway Analysis” for more details of gage package:
http://bioconductor.org/packages/rel...html/gage.html
To understand how GAGE tests for significant pathways, check the paper at: http://www.biomedcentral.com/1471-2105/10/161. HTH.
bigmw is offline   Reply With Quote
Old 10-20-2013, 11:22 PM   #4
sindrle
Senior Member
 
Location: Norway

Join Date: Aug 2013
Posts: 266
Default

Thanks!
Ill check it out.
sindrle is offline   Reply With Quote
Old 10-21-2013, 12:43 AM   #5
sindrle
Senior Member
 
Location: Norway

Join Date: Aug 2013
Posts: 266
Default

"Notice that the native GAGE/Pathview workflow may give you more sensitive results"

How come?
sindrle is offline   Reply With Quote
Old 10-21-2013, 08:28 AM   #6
bigmw
Senior Member
 
Location: US

Join Date: Aug 2013
Posts: 123
Default

The joint workflows are there for users’ convenience. You can do the differential expression analysis (at individual gene level) with well established tools like DESeq2/DESeq, edgeR, Limma and Cufflinks and input the results into GAGE/Pathview workflow for pathway analysis and visualization.

The reason the native GAGE/Pathview workflow is more sensitive is due to the special design of GAGE method. For two-state comparison, it does pair-wise comparison between experiments vs controls (disease vs normal etc), i.e. conducts gene set or pathway test between each sample pairs and then summarizes the results. The more samples (hence independent experiment-control sample pairs) you have, the more testing power you get this way. In the joint workflow, we can’t really decompose the two groups into experiment-control pairs, and we can only take the log2 fold changes from upstream analysis, which is a summary or average fold change without sample size or experiment design info. Therefore, gene set test on this is single column data is not as powerful as the GAGE default pair-wise comparison scheme. In addition, you can analyze and visualize pathway changes in every single experiment or sample (Figure 2 in vignette).
To see the difference between the two types of workflows, you can try the native workflow out and compare the results to joint workflows. You may also read the GAGE paper for more details, at http://www.biomedcentral.com/1471-2105/10/161.

Last edited by bigmw; 10-21-2013 at 08:31 AM.
bigmw is offline   Reply With Quote
Old 10-21-2013, 08:30 AM   #7
sindrle
Senior Member
 
Location: Norway

Join Date: Aug 2013
Posts: 266
Default

Ah, well thanks alot! I will redo my analysis after I have read some more about it.
sindrle is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 07:12 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO