Seurat FindMarkers() output interpretation

Question

I am using FindMarkers() between 2 groups of cells, my results are listed but i'm having hard time in choosing the right markers. Do I choose according to both the p-values or just one of them? If one of them is good enough, which one should I prefer? I am interested in the marker-genes that are differentiating the groups, so what are the parameters i should look for?
             p_val avg_logFC pct.1 pct.2  p_val_adj
UBD   3.080608e-06 1.7755753     1 0.000 0.07683959
CFB   1.262305e-05 1.3989233     1 0.067 0.31485675
RGS13 1.548593e-05 0.9009480     1 0.200 0.38626551
PYDC1 1.896309e-05 0.9622537     1 0.133 0.47299636
BIRC3 4.099414e-05 1.0129472     1 1.000 1.00000000
ICAM1 4.379895e-05 0.8219610     1 0.200 1.00000000

as you can see, p-value seems significant, however the adjusted p-value is not. and when i performed the test i got this warning In wilcox.test.default(x = c(BC03LN_05 = 0.249819542916203,  ... : cannot compute exact p-value with ties
I am completely new to this field, and more importantly to mathematics. Please help me understand in an easy way.

geek_y · Answer

The p-values are not very very significant, so the adj. p-value. You need to plot the gene counts and see why it is the case. It could be because they are captured/expressed only in very very few cells. VlnPlot or FeaturePlot functions should help. You haven't shown the TSNE/UMAP plots of the two clusters, so its hard to comment more.

Answered by geek_y on October 3, 2021

Seurat FindMarkers() output interpretation

One Answer

Add your own answers!

Related Questions

maftools tcgaCompare plot: x-axis labels getting clipped

scRNA-Seq: Account for sequencing depth and gene length?

How to get a list of genes corresponding to the list of SNPs (rs ids)?

Merge / Reconciliate several de novo transcriptome assemblies with different kmers

How to extract Genbank sequence and identify ORFs from the sequence in R?

How are the values of prop.part() and prop.clades() calculated?

Principle of TMT Tags in Multiplex Proteomics

How does picard’s MarkDuplicate handle unmapped reads?

Can Maven (mass spectrometry software) read mzML files?

Java error when trying to run UBCG jar (java.lang.ClassNotFoundException)

How can i get data set of illness and not illness people?

How I change a function not to need an input

DESeqDataSetFromTximport all(lengths > 0) is not TRUE

Download multiple fastq files using fastq-dump

How do I get cromwell to create soft-links?

VCF spec: is it possible to have other alleles in addition to the MISSING value (‘.’) in ALTS?

How can I restrict the search space when querying GEO?

Adding treatment groups via metadata to Seurat object?

Human Herpes Virus, fusion motif in protein vs. SARS-CoV-2?

Identification of starting point for the synthesis of a protein

Ask a Question