SEQanswers

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 01:08 AM
DEXSeq output not convincing r_sitaram RNA Sequencing 10 05-21-2012 04:15 AM
DEXSeq - able to analyze paired (matched) data? spabinger RNA Sequencing 13 02-08-2012 10:40 AM
DEXSeq vs htseq-count/DESeq counting model jdsv Bioinformatics 2 11-20-2011 07:48 PM
How to output Normalised count data from DESeq? fabrice Bioinformatics 1 08-16-2011 11:02 PM

Reply
 
Thread Tools
Old 09-19-2012, 07:09 AM   #1
elizzybethy
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
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 01:51 PM.


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