
Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
how to do test on RNAseq data like this  shuixiangyuer  RNA Sequencing  1  01172013 10:27 AM 
RNASeq: The bench scientist's guide to statistical analysis of RNASeq data.  Newsbot!  Literature Watch  0  09182012 03:00 AM 
RNASeq: A Statistical Framework for eQTL Mapping Using RNAseq Data.  Newsbot!  Literature Watch  0  08162011 03:00 AM 
RNASeq: A survey of statistical software for analysing RNAseq data.  Newsbot!  Literature Watch  12  12202010 10:10 PM 
PubMed: Statistical Design and Analysis of RNASeq Data.  Newsbot!  Literature Watch  0  05092010 08:00 PM 

Thread Tools 
01242013, 08:14 AM  #1 
Junior Member
Location: Odense Join Date: Aug 2012
Posts: 5

Which statistical test: GSEA on duplicated RNASeq data
Hi everybody
I'm struggling a bit trying to do GSEA of RNASeq data. I've ended up settling with a package known as GAGE (Generally Applicable Geneset Enrichment). The main reason is that this algorithm is the only one I've been able to find that does not require 10+ biological replicates. The algorithms employed by GAGE is targetted toward microarray data and as such there are some adjustments that are necessary prior to analysis. Basically, I need to do a transformation to make the data homoscedastic and I need to do length bias correction.I'm then asking GAGE to do a paired comparison between treatment and control and for each pair this will give me some enrichment score. What I'm struggling with is what sort of method I should ask GAGE to use for statistical testing, as I have only two replicates. The options are: I think the correct test to use is the rankbased two sample ttest, but it would be nice if someone with more statistical knowledge could comment on my workflow. Last edited by JesperGrud; 01242013 at 08:16 AM. Reason: Layout of post 
05022013, 10:29 PM  #2 
Senior Member
Location: MDC, Berlin, Germany Join Date: Oct 2009
Posts: 317

Hi, If you haven't solved this out, I'd like to suggest you try our newly developed R package SeqGSEA, which is available at http://bioconductor.org/packages/rel...l/SeqGSEA.html.
It can integrate differential expression and differential splicing together for GSEA. By using negative binomial to model read counts, SeqGSEA can correctly capture technical and biological variance in RNASeq data. It doesn't require a large number of biological replicates, but I'd like to know how many you have. Cheers Xi
__________________
Xi Wang 
05032013, 01:19 AM  #3 
Junior Member
Location: Odense Join Date: Aug 2012
Posts: 5

Hi.
In the end i used a rankbased Wilcoxontest. This might yield some false negative in the sense that pvalue scoring is too conservative, if the data are in fact normally distributed. However, since im only interested in the say top 5 pathways from a biological point of view, it does not matter that much. Our current setup is to have just 2 replicates from cell culture experiments. Our tests show that 2 replicates gives us sufficient data to detect most differential expression using DESeq. I havnt read the details concerning your package, but the questions that spring to mind is if it requires paired end reads? and if the differential splicing analysis can be omitted? 
05042013, 02:40 AM  #4  
Senior Member
Location: MDC, Berlin, Germany Join Date: Oct 2009
Posts: 317

Quote:
I am afraid 2 replicates are too few for SeqGSEA, as its based on sample label permutation for statistical significance. Regarding your questions, SeqGSEA doesn't require PE reads. It only takes readcount data. SeqGSEA can work with DEonly GSEA. Hope this info would help with your future data analysis and other researchers. Cheers.
__________________
Xi Wang 

Thread Tools  

