Seurat繪製單細胞千里江山圖（山巒圖）

山巒圖，英文名叫做ridge plot，其本質上是密度圖（density plot）,其實中文並沒有統一的稱呼，但是由於其外貌高低不平，形似山川，一般稱呼其爲“山巒圖”。山巒圖的用法與小提琴圖、箱形圖相似，主要用於可視化數據的分佈情況，而外形上卻具有中國山水畫的美感。

library(Seurat)
library(SeuratData)
AvailableData()
library(ggthemes)
library(ggimage)

library(pbmc3k.SeuratData) 
pbmc3k.final
An object of class Seurat 
13714 features across 2638 samples within 1 assay 
Active assay: RNA (13714 features, 2000 variable features)
 2 dimensional reductions calculated: pca, umap

先選一個作爲好看的模板。

imcell  <- 'cell.png'
classmk <-c("IL7R", "CCR7",     "IL7R", "S100A4","CD14", "LYZ","MS4A1","CD8A","FCGR3A", "MS4A7",    "GNLY",
            "NKG7","FCER1A", "CST3","PPBP")
RidgePlot(pbmc3k.final,features =classmk) 

看起來，這個幾個基因不錯：

qljs<- c("S100A4","LYZ","CST3","IL7R")

但是最好的還是LYZ，LaiYunZhou的首字母，就它了。

RidgePlot(pbmc3k.final,features ="LYZ") +  # 畫山
  geom_jitter(width = .5,height =.4,alpha=0.5,colour="green",shape=21) +  # 畫雨
  geom_image(image = imcell, size = .025,alpha=0.4)+ # 畫水
  theme_solarized(light = F ) + ggtitle("單細胞·千里江山圖\n天街小雨潤如酥\n草色遙看近卻無") +
  NoLegend()