Question: Seurat "FindMarkers" and "FindallMarkers" v.s. Scanpy "rank_genes_groups"
gravatar for krorange
5 months ago by
krorange0 wrote:

I am processing the same dataset with both Seurat and Scanpy. In Seurat, I got 3 clusters and cluster 2 seems like the target cell type; I got 2 clusters in Scanpy and cluster 1 seems like the target. I am trying to get the marker genes that shows up in both target clusters. But I have two questions...

  1. My understanding is that if I didn't specify ident.2 in Seurat's "FindMarkers", Seurat will find the marker genes comparing the ident.1=0 to the rest, am I right? If so, what's the difference between FindMarkers(object = dataa, ident.1 = 2, min.pct = 0.25), FindMarkers(object = dataa, ident.1 = 2, ident.2 = c(0, 1), min.pct = 0.25) and FindAllMarker(object=dataa, min.pct=0.25)?

  2. Since I'm comparing Seurat result with Scanpy's "", which processing method in question 1 should I compare with?

I'm really confused, it would be helpful if someone can explain these to me.

Thank you so much!

scanpy seurat scrna R single-cell • 884 views
ADD COMMENTlink modified 5 months ago by TriS4.3k • written 5 months ago by krorange0
gravatar for TriS
5 months ago by
United States, Buffalo
TriS4.3k wrote:

part of the answer can be found here for Seurat: FindConservedMarkers vs FindMarkers vs FindAllMarkers Seurat

for Scanpy's code (

    The default method is `'t-test'`,
    `'t-test_overestim_var'` overestimates variance of each group,
    `'wilcoxon'` uses Wilcoxon rank-sum,
    `'logreg'` uses logistic regression. See [Ntranos18]_,
    `here <>`__ and `here
    for why this is meaningful.

looks like the approaches are not that different, and Scanpy's rank_genes_group is similar to Seurat FindMarkers

ADD COMMENTlink modified 5 months ago • written 5 months ago by TriS4.3k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2336 users visited in the last hour