Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
DEXSeq splicing plot from CuffDiff output jasmin9x Bioinformatics 1 06-25-2012 02:08 AM
DEXSeq output not convincing r_sitaram RNA Sequencing 10 05-21-2012 05:15 AM
DEXSeq - able to analyze paired (matched) data? spabinger RNA Sequencing 13 02-08-2012 11:40 AM
DEXSeq vs htseq-count/DESeq counting model jdsv Bioinformatics 2 11-20-2011 08:48 PM
How to output Normalised count data from DESeq? fabrice Bioinformatics 1 08-17-2011 12:02 AM

Thread Tools
Old 09-19-2012, 08:09 AM   #1
Junior Member
Location: San Francisco, CA

Join Date: Sep 2009
Posts: 8
Default Can DEXSeq output variance stabilized data like DESeq?

I am interested in how exposure to a certain drug affects gene expression levels and splicing in cell lines. We have cell lines derived from over 100 different (human) individuals and have exposed each to control buffer and drug in vitro, so pairing our RNA-seq data when making any differential expression/splicing comparisons is vital. (There is so much variability in expression levels between cell lines that it swamps the drug-induced expression changes, even if every cell line changes in the same direction, if the data are not paired for the comparison.)

I found that outputting variance stabilized data from DESeq and then conducting a simple paired t-test was much more powerful than any other differential gene expression approach I have tried thus far for this dataset (including using cell line identity as a covariate in the generalized linear model.) I would like to try a similar approach using exon-level data to investigate splicing, and was wondering if there is a way to export variance-stabilized exon-level data from DEXSeq, and if so, how? (I am a relative R novice, so example code would be much appreciated.)

If there isn't a way to do this in DEXSeq, and one isn't in the pipeline for being implemented soon, I am wondering if there are any other suggestions for conducting differential splicing/exon analysis without sacrificing the inherent power in our paired dataset?
elizzybethy is offline   Reply With Quote

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 02:46 AM.

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