Transcripts with beta close to zero are considered differentially expressed in a likelihood-ratio test in sleuth
1
0
Entering edit mode
4.2 years ago
elsoja ▴ 140

I'm comparing the results that I obtain when doing a DE analysis with the Wald test and the likelihood-ratio test. One the thing that I've noticed is that there are many genes with 'beta' close to zero that are considered differentially expressed between the conditions.

I know that the likelihood-ratio test does't use the beta values to calculate the p-values, but I find it strange that transcripts with similar expression values between the conditions are being considered differentially expressed.

Volcano plot (Wald test q-values):

enter image description here

Volcano plot (likelihood-ratio test q-values):

enter image description here

sleuth RNA-Seq differential-expression • 1.8k views
ADD COMMENT
0
Entering edit mode

What sort of design do you have?

ADD REPLY
0
Entering edit mode

It's a simple design with just one condition: Tissue A (6 samples) vs. Tissue B (3 samples)

ADD REPLY
0
Entering edit mode

I wonder if the plots would look similar if you changed the y-axis of the Wald-test results.

ADD REPLY
0
Entering edit mode

They don't

enter image description here

ADD REPLY
0
Entering edit mode

I suspect you'll need to contact the authors.

ADD REPLY
0
Entering edit mode

I've done this already. Waiting for an answer. Thank you!

https://groups.google.com/forum/#!topic/kallisto-sleuth-users/bRb4nRhoZZw

ADD REPLY
1
Entering edit mode
4.1 years ago
elsoja ▴ 140

This question was answered here

ADD COMMENT

Login before adding your answer.

Traffic: 2739 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