单细胞UMAP降维图的修饰(一些偏僻问题)

2024-03-26 20:50

本文主要是介绍单细胞UMAP降维图的修饰(一些偏僻问题),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

本来呢,单细胞UMAP或者TSNE降维图正常展示就好了,最多修改颜色坐标就可以,但是奈何总是有文章挑战别人的眼睛,非要秀一波非常规操作,那我们也只能把他们打下神坛。
事情的起因是看了这篇文章中的UMAP图,有啥特点呢?每个点都是有边框的,这不就是ggplot中选择shape=21的形状吗?很显然,Dimplot是实现不了的。UMAP降维图本质也是散点图,只需要将作图数据导出,ggplot2就可以实现任何你想要的修饰了。

首先我们设置下颜色,并将作图的数据导出,导出的数据包含UMAP两个坐标的数据。


allcolour=c("#20B2AA","#FFA500","#9370DB","#98FB98","#1E90FF","#7CFC00","#FFFF00","#808000","#FF00FF","#FA8072","#7B68EE","#9400D3","#800080","#A0522D","#D2B48C","#D2691E","#87CEEB","#40E0D0","#5F9EA0","#FF1493","#0000CD","#008B8B","#FFE4B5","#8A2BE2","#228B22","#E9967A","#4682B4","#32CD32","#F0E68C","#FFFFE0","#EE82EE","#FF6347","#6A5ACD","#9932CC","#8B008B","#8B4513","#DEB887")
umap = NCsce@reductions$umap@cell.embeddings %>% as.data.frame() %>%cbind(cell_type = NCsce@meta.data$celltype)

然后就可以绘图了,很简单:


ggplot(umap, aes(x= UMAP_1 , y = UMAP_2 ,fill = cell_type))+geom_point(size = 3.5,colour="grey35",shape=21) +scale_fill_manual(values = allcolour)+ theme_classic()

可以像有些文章中一样,干脆不显示坐标轴,就是单纯的图。

ggplot(umap, aes(x= UMAP_1 , y = UMAP_2 ,fill = cell_type))+geom_point(size = 3,colour="grey35",shape=21) +scale_fill_manual(values = allcolour)+ theme_bw()+theme(panel.grid =element_blank(),panel.border = element_blank(),axis.text.y = element_blank(),axis.ticks.y = element_blank(),axis.ticks.x = element_blank(),axis.text.x = element_blank())+labs(x='', y='')

好了,这里也就是仅仅说明有这么一种情况存在,想要继续修饰可以发挥你们的聪明才智,争取有一个惊艳众人的可视化结果。

这篇关于单细胞UMAP降维图的修饰(一些偏僻问题)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/849808

相关文章

解决pandas无法读取csv文件数据的问题

《解决pandas无法读取csv文件数据的问题》本文讲述作者用Pandas读取CSV文件时因参数设置不当导致数据错位,通过调整delimiter和on_bad_lines参数最终解决问题,并强调正确参... 目录一、前言二、问题复现1. 问题2. 通过 on_bad_lines=‘warn’ 跳过异常数据3

解决RocketMQ的幂等性问题

《解决RocketMQ的幂等性问题》重复消费因调用链路长、消息发送超时或消费者故障导致,通过生产者消息查询、Redis缓存及消费者唯一主键可以确保幂等性,避免重复处理,本文主要介绍了解决RocketM... 目录造成重复消费的原因解决方法生产者端消费者端代码实现造成重复消费的原因当系统的调用链路比较长的时

深度解析Nginx日志分析与499状态码问题解决

《深度解析Nginx日志分析与499状态码问题解决》在Web服务器运维和性能优化过程中,Nginx日志是排查问题的重要依据,本文将围绕Nginx日志分析、499状态码的成因、排查方法及解决方案展开讨论... 目录前言1. Nginx日志基础1.1 Nginx日志存放位置1.2 Nginx日志格式2. 499

kkFileView启动报错:报错2003端口占用的问题及解决

《kkFileView启动报错:报错2003端口占用的问题及解决》kkFileView启动报错因office组件2003端口未关闭,解决:查杀占用端口的进程,终止Java进程,使用shutdown.s... 目录原因解决总结kkFileViewjavascript启动报错启动office组件失败,请检查of

SpringBoot 异常处理/自定义格式校验的问题实例详解

《SpringBoot异常处理/自定义格式校验的问题实例详解》文章探讨SpringBoot中自定义注解校验问题,区分参数级与类级约束触发的异常类型,建议通过@RestControllerAdvice... 目录1. 问题简要描述2. 异常触发1) 参数级别约束2) 类级别约束3. 异常处理1) 字段级别约束

Python错误AttributeError: 'NoneType' object has no attribute问题的彻底解决方法

《Python错误AttributeError:NoneTypeobjecthasnoattribute问题的彻底解决方法》在Python项目开发和调试过程中,经常会碰到这样一个异常信息... 目录问题背景与概述错误解读:AttributeError: 'NoneType' object has no at

Spring的RedisTemplate的json反序列泛型丢失问题解决

《Spring的RedisTemplate的json反序列泛型丢失问题解决》本文主要介绍了SpringRedisTemplate中使用JSON序列化时泛型信息丢失的问题及其提出三种解决方案,可以根据性... 目录背景解决方案方案一方案二方案三总结背景在使用RedisTemplate操作redis时我们针对

Kotlin Map映射转换问题小结

《KotlinMap映射转换问题小结》文章介绍了Kotlin集合转换的多种方法,包括map(一对一转换)、mapIndexed(带索引)、mapNotNull(过滤null)、mapKeys/map... 目录Kotlin 集合转换:map、mapIndexed、mapNotNull、mapKeys、map

nginx中端口无权限的问题解决

《nginx中端口无权限的问题解决》当Nginx日志报错bind()to80failed(13:Permissiondenied)时,这通常是由于权限不足导致Nginx无法绑定到80端口,下面就来... 目录一、问题原因分析二、解决方案1. 以 root 权限运行 Nginx(不推荐)2. 为 Nginx

解决1093 - You can‘t specify target table报错问题及原因分析

《解决1093-Youcan‘tspecifytargettable报错问题及原因分析》MySQL1093错误因UPDATE/DELETE语句的FROM子句直接引用目标表或嵌套子查询导致,... 目录报js错原因分析具体原因解决办法方法一:使用临时表方法二:使用JOIN方法三:使用EXISTS示例总结报错原