比如要畫(huà)某個(gè)marker基因的featureplot
DefaultAssay(sample.integrated) <- 'RNA'
FeaturePlot(sample.integrated,features = 'LMNA',reduction = 'tsne')
如果我想要分組畫(huà)
FeaturePlot(sample.integrated,features = 'LMNA',reduction = 'tsne',split.by = 'group')
這時(shí)候會(huì)發(fā)現(xiàn)圖例不見(jiàn)了
添加圖例:
FeaturePlot(sample.integrated,features = 'LMNA',reduction = 'tsne',split.by = 'group')& theme(legend.position = "right")