<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>

CC BY-SA 4.0
2021-01-21T10:16:46.857Z
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```
StupidWolf
CC BY-SA 4.0
2021-01-22T11:54:12.513Z
@StupidWolf Yes, I did set the identity to group. The function works for other cell types in the same data set.
Charles
<scrnaseq><seurat>
Why could I be having only NAs for p_val and p_val_adj in Seurat DGE analysis?
# 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**: @StupidWolf Yes, I did set the identity to group. The function works for other cell types in the same data set.  

---

## 0 answers

