This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Extracting Statistically Significant Differentially Expressed Exons From Dexseq Results

Hi

I am trying to extract some data from my DEXseq results. I was able to get the final HTML report for my analysis. However, I was wondering if there's a way I could extract specific exons that are significantly differentially expressed between two conditions. I am interested in getting the gene id (ensembl gene id ok), exon or exon bin id (if applicable), their exon boundaries, fold change & adjusted p value. I am a R/bioconductor newbie. Any help is appreciated.

Thanks

rna-seq bioconductor r splicing
• 5,306 views • link
updated by Fahmi • Fahmi • written by pm2013 • pm2013

2 answers

In part, you can get that from res <- DEUresultTable(ecs), where ecs is an ExonCountSet. The coordinates are in featureData(ecs). You can then get the log2 foldchange with estimatelog2FoldChanges(ecs, ...)

• 74 views • link
written by Devon Ryan Devon Ryan

Thank you for your comment. I will try your suggestions. Looks like there's no straightforward way to this.

• 0 views • link
written by pm2013 • pm2013

Hello, how could I get the statistically significant for DEXSeq result. The output which I got from html only the sorted result based on Ensembl name without the p-value or q-value. How do I get that result? Thank you

• 0 views • link
written by Fahmi • Fahmi

Get the non-HTML results. If you didn't do the analysis yourself then ask whomever did do it for a csv file or something like that.

• 0 views • link
written by Devon Ryan Devon Ryan

Is it from the dxr result? Do you mind to show me the code for this. I tried before didn't get it.

• 0 views • link
updated by Ram Ram • written by Fahmi • Fahmi

Assuming "dxr" is the output of DEXSeqResults() then yes, that's a table with everything in it. Just write it to a file.

• 0 views • link
written by Devon Ryan Devon Ryan

Log in to answer this question.

More posts like this