完成單細胞定群后，就可以做很多東西了，首先可視化marker基因，這是很多單細胞文章中不可或缺的一部分，當(dāng)然，marker基因的可視化可以拓展到對任意自己感興趣基因的可視化。

小提琴圖是其中的一種。Seurat自帶函數(shù)VlnPlot可以進行基因可視化，其實在前面QC結(jié)果的戰(zhàn)事中已經(jīng)使用過一次了。

VlnPlot(scedata, 
        features = c("RAMP2","VWF","PECAM1","LUM"),
        pt.size = 0,
        ncol = 2) 

image

通過添加split.by參數(shù)，VlnPlot可以顯示基因在各個樣本中的表達情況。

VlnPlot(scedata, 
        features = c("RAMP2","VWF","PECAM1","LUM"),
        pt.size = 0,
        ncol = 2,
        split.by = "orig.ident")

image

但是，VlnPlot有一個問題就是，只能展示單個圖，不能批量展示基因的表達。不能實現(xiàn)很多文章中出現(xiàn)堆疊的效果。這就需要額外的修飾改造達到這種效果。這里我們提供兩種方法可以近似獲得堆疊小提琴圖效果（圖A，B），但細節(jié)問題還需進一步修飾。

image

image

更多精彩請關(guān)注---《KS科研分享與服務(wù)》