在單細胞的研究中醋粟,經(jīng)常需要做的就是聚類靡菇。在聚類的時候,對于聚出多少個cluster,一個相關(guān)的參數(shù)是resolution,當(dāng)數(shù)值小的時候袄简,cluster少混槐,數(shù)值大的時候,cluster多立倍。當(dāng)我們想看到隨著resolution變化,cluster變化的過程朝卒,可以考慮用clustree這個R包可視化博烂。
(為單細胞分析聚類可視化只是其中的一個功能香椎,其tutorial如下:https://cran.r-project.org/web/packages/clustree/vignettes/clustree.html)
# === 以上省略Seurat包進行單細胞分析的過程====
library(clustree)
test <- FindClusters(FrontalCortex_counts, resolution = seq(0.1, 2, 0.5))
clus.tree.out <- clustree(test)
pdf(file = "test.pdf", width = 15, height = 10)
print(clus.tree.out)
dev.off()
可以得到:
可以看到哪些群會隨著resolution的增大而繼續(xù)分裂增多。