I’m a new in Scanpy and impressed by its speed and user-friendly. Since I used to be a big fan of Seurat, the most popular R package for snRNA-seq analysis, I don’t know how to do some operations I often do in Seurat with Scanpy. The next is an example.
In Seurat, I could get the average gene expression of each cluster easily by the code showed in the picture.
I think Scanpy can do the same thing as well, but I don’t know how to do right now. I would be grateful if anyone could offer me some help.