Question: DESeq result "LFC > 0 (up) : 0, 0% ,LFC > 0 (down) : 0, 0%"
gravatar for Hughie
2.8 years ago by
Hughie80 wrote:

Hello! All.
I'm using DESeq to check differential gene expression , but I got in doubt recent days and meet strange result which different from DESeq munual's demo , below I post my code and wish your kindly help:

workDir <- "C:/Users/Administrator/Desktop/rawcounts"
sampleFiles <- grep(".reads",list.files(directory),value=TRUE)
stage <- factor(c("B","CD4","CD8","CLP","CMP","EryA","EryB","GMP","Granulocyte","HSC","LT_HSC",
sampleTable<-data.frame(sampleName=sampleFiles, fileName=sampleFiles,stage = stage)
dds<-DESeqDataSetFromHTSeqCount(sampleTable=sampleTable, directory=directory,design =~ stage)
dds <- dds[rowSums(counts(dds)) > 1,]
resOrdered <- res[order(res$padj),] 

And I got the result:


Note:My datas have no replicates, so I wonder if this is the problem.
Thank you again for your attention!

rna-seq • 1.4k views
ADD COMMENTlink modified 2.8 years ago • written 2.8 years ago by Hughie80
gravatar for Santosh Anand
2.8 years ago by
Santosh Anand5.2k
Santosh Anand5.2k wrote:

IMHO, you need some replicates to have a meaningful comparison. And yes, the FDR/p-values depend very much on replicates, so this might explain the things. You may see this post from Simon Anders for further details regarding running DESeq2 w/o replicates

Apart from that, you also need to choose which comparisons are you trying to make. In the default setting, the results() returns the comparison of the last level of the last variable in the design formula over the first level of this variable. For example, for a simple two-group comparison, this would return the log2 fold changes of the second group over the first group (the reference level).

see ?results for details

Also, see this for multifactor design

ADD COMMENTlink modified 2.8 years ago • written 2.8 years ago by Santosh Anand5.2k

Thank you! Santosh. These two pages are very useful and I will try more. BTW,except edgeR, is there any packages can be recommended for DEG analysis. Thank you again for your reply.

ADD REPLYlink written 2.8 years ago by Hughie80

Newer and faster alternative pipeline is to use transcript abundance quantification methods such as Salmon (Patro et al. 2017), Sailfish (Patro, Mount, and Kingsford 2014), kallisto (Bray et al. 2016), or RSEM (B. Li and Dewey 2011), to estimate abundances without aligning reads.

Also see,

ADD REPLYlink written 2.8 years ago by Santosh Anand5.2k

Thank you again! Santosh.

I will try to figure out these packages and the link is very specific.

ADD REPLYlink written 2.8 years ago by Hughie80
Please log in to add an answer.


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