This is done to emulate the plots from flow cytometry.

but using data from single cell rna-seq (seurat object …)

library(Seurat)
library(sctree)

Relies on having done the standard seurat workflow before FindAllMarkers

plot_flowstyle(Seurat::pbmc_small, c("HLA-DPB1", "LYZ", "TYMP"))