Answers.Body: <p>I'm guessing here, but I suspect this is at least in part due to the fact that all your samples in condition 1 are not expressed at all. </p>
  
  <p>Formally you have an infinitely large fold-change between condition1 and condition2. Most tools solve this by adding a small pseudocount to unexpressed samples, but this will still result in a very high fold change. So far for the big fold-change.</p>
  
  <p>If I had to guess at the reason why this is significant I would think in terms of dispersion estimation. <a href="https://support.bioconductor.org/p/100764/" rel="nofollow noreferrer">This Bioconductor post</a> has some discussion with the authors of DESeq2 on how this works.</p>
  
  <p>Finally: calculating differential expression on these genes in these samples will never be really informative, as they really are only expressed in one sample. The way I see it you have two options:</p>
  
  <ol>
  <li>Exclude the sample from your analysis. Straightforward, but does not solve this problem for other samples.</li>
  <li>Filter your genes before calculating differential expression, i.e. remove all genes that are not expressed in a large number of samples. This approach is standard in for instance <a href="https://pachterlab.github.io/sleuth/about" rel="nofollow noreferrer">sleuth</a>. A good starting point would be to remove genes that are not expressed in 75% of your samples.</li>
  </ol>
  
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-11-17T09:44:42.133Z
Answers.LastActivityDate: 2017-11-17T09:44:42.133Z
Answers.UserDisplayName: holmrenser
Body: <p>I am using DESeq2 to perform a differential expression analysis, but I obtained very strange results for some genes. For some genes, I have very high log2FoldChange with very low p-value as displayed in the following table.</p>
  
  <p><a href="https://i.sstatic.net/V35Fp.png" rel="nofollow noreferrer"><img src="https://i.sstatic.net/V35Fp.png" alt="enter image description here"></a></p>
  
  <p>When I look the normalized expression values for these genes.  I can see that these genes are expressed in only one single sample (condition1= samples from 1 to 7, condition2 = samples 8 to 13).</p>
  
  <p><a href="https://i.sstatic.net/Ooejw.png" rel="nofollow noreferrer"><img src="https://i.sstatic.net/Ooejw.png" alt="enter image description here"></a></p>
  
  <p>Could you please tell me how to avoid DESeq2 to detect these genes as differentially expressed ?</p>
  
  <p>I am using to following commands:</p>
  
  <pre><code>dds                      &lt;- DESeqDataSetFromHTSeqCount(sampleTable = sampleTable, directory = directory, design= ~ condition)
  dds                      &lt;- DESeq(dds)
  res                      &lt;- results(dds, contrast=c("condition",condition1,condition2))
  </code></pre>
  
Comments.ContentLicense: CC BY-SA 3.0
Comments.CreationDate: 2017-11-07T11:17:13.750Z
Comments.Text: What are the size factors for these samples (`sizeFactors(dds)`)? I suspect that sample10 is just really wonky and throwing everything off, so you'll probably want to just exclude it.
Comments.UserDisplayName: Devon Ryan
Comments.ContentLicense: CC BY-SA 3.0
Comments.CreationDate: 2017-11-07T12:48:34.273Z
Comments.Text: I have check that but the factor is in the range of other samples... I have similar situations with other samples..
Comments.UserDisplayName: Bob
Comments.ContentLicense: CC BY-SA 3.0
Comments.CreationDate: 2017-11-07T12:52:52.617Z
Comments.Text: Have you made any of the diagnostic plots related to cook's cutoff distance? If so, post some. If not, make some (see the [here](https://bioconductor.org/packages/release/bioc/vignettes/DESeq2/inst/doc/DESeq2.html#count-outlier-detection) for an example).
Comments.UserDisplayName: Devon Ryan
Tags: <differential-expression><r><deseq2>
Title: Why do I obtain very strange results with DESeq2?
ThreadText: # Why do I obtain very strange results with DESeq2?
  
  (Asked by: username_0 on 2017-11-07)
  
  I am using DESeq2 to perform a differential expression analysis, but I obtained very strange results for some genes. For some genes, I have very high log2FoldChange with very low p-value as displayed in the following table.
  
  [![enter image description here](https://i.sstatic.net/V35Fp.png)](https://i.sstatic.net/V35Fp.png)
  
  When I look the normalized expression values for these genes. I can see that these genes are expressed in only one single sample (condition1= samples from 1 to 7, condition2 = samples 8 to 13).
  
  [![enter image description here](https://i.sstatic.net/Ooejw.png)](https://i.sstatic.net/Ooejw.png)
  
  Could you please tell me how to avoid DESeq2 to detect these genes as differentially expressed ?
  
  I am using to following commands:
  
  ```
  dds                      <- DESeqDataSetFromHTSeqCount(sampleTable = sampleTable, directory = directory, design= ~ condition)
  dds                      <- DESeq(dds)
  res                      <- results(dds, contrast=c("condition",condition1,condition2))
  
  ```
  
  **Tags:** differential-expression, r, deseq2
  
  > **Comments:**
  > **username_1**: What are the size factors for these samples (`sizeFactors(dds)`)? I suspect that sample10 is just really wonky and throwing everything off, so you'll probably want to just exclude it.  
  > **username_0**: I have check that but the factor is in the range of other samples... I have similar situations with other samples..  
  > **username_1**: Have you made any of the diagnostic plots related to cook's cutoff distance? If so, post some. If not, make some (see the [here](https://bioconductor.org/packages/release/bioc/vignettes/DESeq2/inst/doc/DESeq2.html#count-outlier-detection) for an example).  
  
  ---
  
  ## 1 answers
  
  ### Answer by username_2 on 2017-11-17. Score: 2
  
  I'm guessing here, but I suspect this is at least in part due to the fact that all your samples in condition 1 are not expressed at all.
  
  Formally you have an infinitely large fold-change between condition1 and condition2. Most tools solve this by adding a small pseudocount to unexpressed samples, but this will still result in a very high fold change. So far for the big fold-change.
  
  If I had to guess at the reason why this is significant I would think in terms of dispersion estimation. [This Bioconductor post](https://support.bioconductor.org/p/100764/) has some discussion with the authors of DESeq2 on how this works.
  
  Finally: calculating differential expression on these genes in these samples will never be really informative, as they really are only expressed in one sample. The way I see it you have two options:
  
  1. Exclude the sample from your analysis. Straightforward, but does not solve this problem for other samples.
  2. Filter your genes before calculating differential expression, i.e. remove all genes that are not expressed in a large number of samples. This approach is standard in for instance [sleuth](https://pachterlab.github.io/sleuth/about). A good starting point would be to remove genes that are not expressed in 75% of your samples.
  
  
  ---
  
  