CNS图表复现04—单细胞聚类分群的resolution参数问题
时间:2022-07-25
本文章向大家介绍CNS图表复现04—单细胞聚类分群的resolution参数问题,主要内容包括其使用实例、应用技巧、基本知识点总结和需要注意事项,具有一定的参考价值,需要的朋友可以参考一下。
分享是一种态度
回顾
我们的CNS图表复现之旅已经开始,前面3讲是;
如果你也想加入交流群,自己去:你要的rmarkdown文献图表复现全套代码来了(单细胞)找到我们的拉群小助手哈。
眼尖的同学发现,我的代码里面,关于单细胞聚类分群的resolution参数一直在变化:
sce <- FindNeighbors(sce, dims = 1:15)
sce <- FindClusters(sce, resolution = 0.2)
table(sce@meta.data$RNA_snn_res.0.2)
sce <- FindClusters(sce, resolution = 0.8)
table(sce@meta.data$RNA_snn_res.0.8)
sce <- FindClusters(sce, resolution = 0.5)
table(sce@meta.data$RNA_snn_res.0.5)
然后提出疑问,这个resolution参数到底应该是选择多少呢?难道这个步骤没有一个绝对的标准吗?
我之前给大家举例是使用balloonplot这个可视化函数,代码如下:
library(gplots)
tab.1=table(sce@meta.data$RNA_snn_res.0.2,sce@meta.data$RNA_snn_res.0.8)
balloonplot(tab.1)
就可以很直观的看到,我们把resolution参数分别赋值为0.2和0.8的效果,如下:
0.2和0.8的效果对比
很明显,这个resolution越小呢,我们得到的分群数量就越少,所以0.2的时候是17个群,但是0.8的时候是31个群。
而且我们根据balloonplot的可视化,可以看到,在0.2的时候的17个群里面的有一些群,会随着resolution的调高,继续裂变成为多个群。
有意思的是,我恰好在你要的rmarkdown文献图表复现全套代码来了(单细胞) 看到了一个更好的可视化方法:
# Check clustering stability at given resolution
# Set different resolutions
res.used <- seq(0.1,1,by=0.2)
res.used
# Loop over and perform clustering of different resolutions
for(i in res.used){
sce <- FindClusters(object = sce, verbose = T, resolution = res.used)
}
# Make plot
library(clustree)
clus.tree.out <- clustree(sce) +
theme(legend.position = "bottom") +
scale_color_brewer(palette = "Set1") +
scale_edge_color_continuous(low = "grey80", high = "red")
clus.tree.out
就是借用clustree包,可视化如下:
可以非常清晰的看到,随着resolution的调高,具体是哪些群在不停地继续裂变成为多个群。
但是呢, 仍然是没有回答粉丝的问题,就是resolution设置多少,难道说没有一个绝对的指标吗?
我这里只能说,确实没有,不仅仅是resolution参数,生物信息学数据分析过程中,就比如这个单细胞吧,质控的时候去除多少个质量差的细胞去除多少基因,选择高变基因数量多少,PCA降维后选择多少个PC,基本上每个步骤都是可以灵活调整的。
这就是,数据分析的魅力吧。
- JavaScript 教程
- JavaScript 编辑工具
- JavaScript 与HTML
- JavaScript 与Java
- JavaScript 数据结构
- JavaScript 基本数据类型
- JavaScript 特殊数据类型
- JavaScript 运算符
- JavaScript typeof 运算符
- JavaScript 表达式
- JavaScript 类型转换
- JavaScript 基本语法
- JavaScript 注释
- Javascript 基本处理流程
- Javascript 选择结构
- Javascript if 语句
- Javascript if 语句的嵌套
- Javascript switch 语句
- Javascript 循环结构
- Javascript 循环结构实例
- Javascript 跳转语句
- Javascript 控制语句总结
- Javascript 函数介绍
- Javascript 函数的定义
- Javascript 函数调用
- Javascript 几种特殊的函数
- JavaScript 内置函数简介
- Javascript eval() 函数
- Javascript isFinite() 函数
- Javascript isNaN() 函数
- parseInt() 与 parseFloat()
- escape() 与 unescape()
- Javascript 字符串介绍
- Javascript length属性
- javascript 字符串函数
- Javascript 日期对象简介
- Javascript 日期对象用途
- Date 对象属性和方法
- Javascript 数组是什么
- Javascript 创建数组
- Javascript 数组赋值与取值
- Javascript 数组属性和方法
- 多线程应用 - 阻塞队列ArrayBlockingQueue详解
- Java虚拟机 - 超级详细的类加载说明
- ReentrantLock
- 多线程应用 - 阻塞队列LinkedBlockingDeque详解
- AtomicInteger详解
- 多线程应用 - 阻塞队列LinkedBlockingQueue详解
- 多线程应用 - 超详细的AQS详情
- 多线程应用 - 基于AQS的Condition
- HashMap源码分析 - jdk8
- HashMap源码分析 - JDK7和JDK8有什么区别
- 线程应用 - 中断interrupt详解
- ThreadLocal全面解析
- tomcat类加载机制了解一下
- tomcat-超详细的启动流程(init)
- tomcat-超详细的启动流程(start)