Fold Change for rnaseq DTU
1
0
Entering edit mode
rbenel ▴ 20
@rbenel-13642
Last seen 3 days ago
Israel

Hi!

I have been using the rnaseqDTU workflow and the function testForDEU internally uses the DESeq2 function nbinormLRT. I was wondering what the ramifications are if I wanted to calculate fold change for example for specific timepoints.. Hypothetically could I use the wald test instead of LRT in this case to get FCs?

Thanks!

dxd <- estimateSizeFactors(dxd) 
  dxd <- estimateDispersions(dxd, quiet=FALSE)
  dxd <- testForDEU(dxd, reducedModel=~sample + exon + replicate:exon)
DEXSeq rnaseqDTU • 213 views
ADD COMMENT
1
Entering edit mode
@mikelove
Last seen 12 hours ago
United States

I’m not sure DEXSeq offers access to those individual LFCs but I’ll defer to the package authors.

ADD COMMENT
0
Entering edit mode

OK. Thanks Michael! Should I add an additional tag to this post?

ADD REPLY
0
Entering edit mode

Any update on this post authors of DEXSeq?

ADD REPLY
1
Entering edit mode

I pinged the authors.

ADD REPLY
0
Entering edit mode

Thank You!

ADD REPLY
0
Entering edit mode

This is a good point. It should be possible to use the Wald test vs LRT to get the LFC of the interaction exon:condition. But I've not compared the different approaches myself. If any (differences) I guess they would be similar to those observed for gene level tests.

ADD REPLY
0
Entering edit mode

OK. Interesting. Thank you for your input!

ADD REPLY

Login before adding your answer.

Traffic: 206 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6