Question: DESeq result "LFC > 0 (up) : 0, 0% ,LFC > 0 (down) : 0, 0%"
0
gravatar for Hughie
14 months ago by
Hughie70
guangzhou
Hughie70 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:

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

And I got the result:

image

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

rna-seq • 709 views
ADD COMMENTlink modified 14 months ago • written 14 months ago by Hughie70
4
gravatar for Santosh Anand
14 months ago by
Santosh Anand4.7k
Santosh Anand4.7k 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

http://seqanswers.com/forums/showpost.php?p=107433&postcount=2

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 https://support.bioconductor.org/p/67600/#67612

ADD COMMENTlink modified 14 months ago • written 14 months ago by Santosh Anand4.7k

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 14 months ago by Hughie70

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, https://liorpachter.wordpress.com/2017/08/02/how-not-to-perform-a-differential-expression-analysis-or-science/

ADD REPLYlink written 14 months ago by Santosh Anand4.7k

Thank you again! Santosh.

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

ADD REPLYlink written 14 months ago by Hughie70
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 2.3.0
Traffic: 1937 users visited in the last hour