Question: removing outliers before running DESeq2
0
Entering edit mode

Dear all,

After running DESeq2 program successfully I found out that I had to detect and remove outliers from my htseq count datasets before running DESeq2. However, I am a little confused about how to perform this task.

Can u advise me the most straightforward way to do this?

Thank u in advance

Nazanin Hosseinkhan

ADD COMMENTlinkeditmoderate 18 months ago nazaninhoseinkhan • 360 • updated 18 months ago arup ♦ 1.3k
Entering edit mode
1

Why do you think you need to remove outliers? I tend to remove the genes that doesnt have more than 5 counts on average across all samples but nothing more.

ADD REPLYlinkeditmoderate 18 months ago
firatuyulur
• 280
Entering edit mode
0

Thank u so much. Yes, I've already removed genes with lower than 10 reads

ADD REPLYlinkeditmoderate 18 months ago
nazaninhoseinkhan
• 360
Entering edit mode
0

Why 10? Why not any other?

ADD REPLYlinkeditmoderate 18 months ago
Arindam Ghosh
• 160
2
Entering edit mode

There is no need to remove out layers but it's better to remove the technical artifacts(noise incorporated by the sequencer). Genes/transcripts having read count within 0-10 rage is generally considered as artifacts(if the expression pattern is inconstant throughout all the samples). For more details check the following thread.

https://support.bioconductor.org/p/95755/

ADD COMMENTlinkeditmoderate 18 months ago arup ♦ 1.3k
Entering edit mode
0

Thanks a lot for your advice.

ADD REPLYlinkeditmoderate 18 months ago
nazaninhoseinkhan
• 360

Login before adding your answer.

Powered by the version 2.0