# 使用DBSCAN算法检测离群点,,离群点通常位于低密度区域。
umap_data <- FetchData(signac_obj_integ_flt,vars = c('umap_1','umap_2','cytpe'))
lapply(unique(umap_data $cytpe), function(x){
data <- umap_data %>% filter(cytpe== x)
### 关键参数: eps[检测范围];minPts[当前范围内最少应包含的细胞数量]
dbscanResult <- dbscan::dbscan(as.matrix(data[,1:2]), eps= 0.3, minPts= 60)$cluster
return({
rownames(data)[dbscanResult > 0]
})
}) %>% unlist -> inlier_Cells
options(repr.plot.width = 12, repr.plot.height = 9)
umap_data[inlier_Cells,] %>% ggplot(aes(x = umap_1, y = umap_2, color = cytpe)) + geom_point(size = 0.1)
利用 dbscan 检测UMAP图中的Outliers
©著作权归作者所有,转载或内容合作请联系作者
禁止转载,如需转载请通过简信或评论联系作者。
- 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
- 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
推荐阅读更多精彩内容
- 这8种学生永远拿不到高分!早看早受益! 下面是一位资深班主任总结了8种成绩提不上去的原因,分别对应8类孩子,如果你...
- 这8种学生永远拿不到高分!早看早受益! 下面是一位资深班主任总结了8种成绩提不上去的原因,分别对应8类孩子,如果你...
- 译者按: 用于OCR领域的经典网络,CTPN是目前流传最广、影响最大的开源文本检测模型,可以检测水平或微斜的文本行...
- 一、背景介绍 为什么前端工程师要用PS切图? 作为前端,我们需要负责把UI设计的图片,实现成HTML网页,在实现过...
- 1.问题描述 要使地图的中心点位置像上偏移,让地图上的轨迹完全显示出来,不被下方的自定义View遮挡。 2.解决方...