I am doing my singlecell analysis in Seurat. I want to create Pearson correlation matrix of a set of genes per cluster basis. I want to see how these set of genes are co-expressed in all clusters.
Something like this(Referenced from Deseq2 manual)
I know about Seurat has FeatureScatter but that's only one gene against another. How to do for multiple genes? Also which data slot is appropriate for this kind of analysis?