DESeq2:Why is there an oblique line in my MA picture, which happened before, what caused this?
1
0
Entering edit mode
@xuyao15927402897-20826
Last seen 17 days ago
Germany

My sample is RNA-Seq data from human Esophageal tissue, 9 control groups(normal), 2 experimental groups(Esophageal cancer),The result is Figure 1.Figure 1

I use featurecounts to count genes, this is my counts array

deseq2 cancer • 345 views
ADD COMMENT
1
Entering edit mode
@mikelove
Last seen 1 day ago
United States

Those lines are from genes with very low counts of genes, e.g. 2 counts in the treated group, 1 count in the control group. If you want to remove them you can perform some minimal filtering, e.g.:

keep <- rowSums(counts(dds) >= 10) >= x
dds <- dds[keep,]

Where you choose x to be a small number of samples.

ADD COMMENT
0
Entering edit mode

Will removing these low counts have any bad effect, will I ignore the low-expression genes?

ADD REPLY
0
Entering edit mode

It's pretty hard to achieve power for DE when there are say, only 2 samples with counts above 10. So setting x <- 3 will be fairly conservative in terms of FNR.

ADD REPLY
0
Entering edit mode

Ok, thank you very much!

ADD REPLY

Login before adding your answer.

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