Body: <p>I am doing a differential gene expression analysis to find out genes that are differentially expressed in one cell type from the lung in two conditions, control and idipathic pulmonary fibrosis (IPF). In total I have 152538 cells. Can anyone help me provide a hint to understand why only NAs get returned for the pvalues and adjusted pvalues?</p>
  <pre><code>&gt; table(seurat$group)
  
  Control     IPF 
    64652   87886 
  
  &gt; degs &lt;- FindMarkers(seurat, 
                              ident.1 = &quot;IPF&quot;,
                              ident.2 = &quot;Control&quot;)
  &gt; head(degs) 
           p_val  avg_logFC pct.1 pct.2 p_val_adj
  SPP1        NA  1.2721607 0.390 0.125        NA
  FN1         NA  1.1768771 0.741 0.409        NA
  CRIP1       NA  1.0484934 0.938 0.832        NA
  CCL18       NA  0.9830185 0.734 0.469        NA
  S100A10     NA  0.5799861 0.965 0.902        NA
  VIM         NA  0.5666287 0.992 0.964        NA
  </code></pre>
  
Comments.ContentLicense: CC BY-SA 4.0
Comments.CreationDate: 2021-01-21T10:16:46.857Z
Comments.Text: are you showing the full function called? it should ```FindMarkers(seurat,group.by="group",ident.1="IPF",ident.2="Control")``` unless you set the identity to be ```group```
Comments.UserDisplayName: StupidWolf
Comments.ContentLicense: CC BY-SA 4.0
Comments.CreationDate: 2021-01-22T11:54:12.513Z
Comments.Text: @StupidWolf Yes, I did set the identity to group. The function works for other cell types in the same data set.
Comments.UserDisplayName: Charles
Tags: <scrnaseq><seurat>
Title: Why could I be having only NAs for p_val and p_val_adj in Seurat DGE analysis?
ThreadText: # Why could I be having only NAs for p_val and p_val_adj in Seurat DGE analysis?
  
  (Asked by: username_0 on 2021-01-20)
  
  I am doing a differential gene expression analysis to find out genes that are differentially expressed in one cell type from the lung in two conditions, control and idipathic pulmonary fibrosis (IPF). In total I have 152538 cells. Can anyone help me provide a hint to understand why only NAs get returned for the pvalues and adjusted pvalues?
  
  ```
  > table(seurat$group)
  
  Control     IPF 
    64652   87886 
  
  > degs <- FindMarkers(seurat, 
                              ident.1 = "IPF",
                              ident.2 = "Control")
  > head(degs) 
           p_val  avg_logFC pct.1 pct.2 p_val_adj
  SPP1        NA  1.2721607 0.390 0.125        NA
  FN1         NA  1.1768771 0.741 0.409        NA
  CRIP1       NA  1.0484934 0.938 0.832        NA
  CCL18       NA  0.9830185 0.734 0.469        NA
  S100A10     NA  0.5799861 0.965 0.902        NA
  VIM         NA  0.5666287 0.992 0.964        NA
  
  ```
  
  **Tags:** scrnaseq, seurat
  
  > **Comments:**
  > **username_1**: are you showing the full function called? it should ```FindMarkers(seurat,group.by="group",ident.1="IPF",ident.2="Control")``` unless you set the identity to be ```group```  
  > **username_0**: @username_1 Yes, I did set the identity to group. The function works for other cell types in the same data set.  
  
  ---
  
  ## 0 answers
  
  