Show a gene set on FeaturePlot

[schelp72]

Hope everyone is well.

I would like to display the average expression level of a set of genes (at least 5) in a FeaturePlot but the function only takes 2 genes max in its argument.

My request is similar to what is asked here:
#528
https://www.biostars.org/p/334617/

I am a beginner to R and do not know how to update this code to the newest Seurat version.

Thanks!

[Eomesodermin]

You can try generating a module score for the gene set:

Seurat.obj <- AddModuleScore(
object = Seurat.obj,
features = list(geneset.vector),
name = 'Module_score_name'
)

then plot the module score in FeaturePlot

[schelp72]

Thank you this is exactly what I needed!