Question: DESeq2:Why is there an oblique line in my MA picture, which happened before, what caused this?
0
gravatar for xuyao15927402897
5 weeks ago by
xuyao159274028970 wrote:

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

cancer deseq2 • 82 views
ADD COMMENTlink modified 5 weeks ago by Michael Love24k • written 5 weeks ago by xuyao159274028970
Answer: DESeq2:Why is there an oblique line in my MA picture, which happened before, wha
1
gravatar for Michael Love
5 weeks ago by
Michael Love24k
United States
Michael Love24k wrote:

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 COMMENTlink written 5 weeks ago by Michael Love24k

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

ADD REPLYlink modified 4 weeks ago • written 4 weeks ago by xuyao159274028970

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 REPLYlink written 4 weeks ago by Michael Love24k

Ok, thank you very much!

ADD REPLYlink written 4 weeks ago by xuyao159274028970
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 243 users visited in the last hour