How can I obtain the percentage gene expression per identity class in Seurat as further processible numbers (e.g. matrix)?
This can be solved like this:

library(Seurat)

my_genes &lt;- c("gene1", "gene2", "gene3")

exp &lt;- FetchData(object, my_genes)

matrix &lt;- as.matrix(colMeans(exp  &gt; 0))*100