percentMarkerOverlap
percentMarkerOverlap(markers1, markers2, mode = "min", marker_col = "gene")
| markers1 | Data frame of markers for one dataset as returned by Seurat::FindAllMarkers() |
|---|---|
| markers2 | Data frame of markers for second dataset as returned by Seurat::FindAllMarkers() |
| mode | "s1", "s2", "min", "max" |
| marker_col | String specifying the column in the markers data frames which specifies the cluster. By default, Seurat calls this "gene" (the default here); in objects produced by the lab's pipeline, it may be called "external_gene_name". |
Data frame with three columns: "s1_cluster" giving the cluster for dataset 1, "s2_cluster" giving the cluster for dataset 2, and "marker_overlap" which gives the computed min percent overlap of markers between the two clusters
percentMarkerOverlap(markers_pbmc, markers_pbmc)#> s1_cluster s2_cluster marker_overlap #> 1 0 0 1.00000000 #> 2 0 1 0.19178082 #> 3 0 2 0.13333333 #> 4 0 3 0.27659574 #> 5 1 0 0.19178082 #> 6 1 1 1.00000000 #> 7 1 2 0.13698630 #> 8 1 3 0.15068493 #> 9 2 0 0.13333333 #> 10 2 1 0.13698630 #> 11 2 2 1.00000000 #> 12 2 3 0.06382979 #> 13 3 0 0.27659574 #> 14 3 1 0.15068493 #> 15 3 2 0.06382979 #> 16 3 3 1.00000000