注释气泡图函数(更新)

之前我们写过一个原创可视化函数Dotplot_anno.R,nature级别图表:一个注释气泡热图函数(适用于单细胞及普通数据)。主要解决的问题是1) 单细胞基因可视化分组注释。2) Bulk RNA差异基因热图、气泡图。3) 富集分析结果气泡图展示。这个函数还是受到很多小伙伴的喜爱的。但是在实际应用过程中发现了一些问题,所以这里我们对这个函数进行一下更新,解决一些问题,之前购买过此函数的小伙伴可微信联系作者获取更新版。微信VIP群已更新!

主要更新的内容是一些之前遗留的疏忽错误,例如气泡过大的问题等等。此外小伙伴反应之后,我们也发现,之前的函数在做单细胞气泡图的时候函数没有提供修改因子顺序,因为order=T设置之后导致气泡是按照从大到小排列的。所以这里我们修改一下,让顺序排列更加灵活。

图片

为了演示,我们换一个数据,其他的内容不变,只演示单细胞的内容。我们先按照默认的做一下,order=T。

library(Seurat)library(ggplot2)library(dittoSeq)source('./Dotplot_anno.R')human_data <- readRDS("D:/KS项目/公众号文章/human_data.rds")DefaultAssay(human_data) <- "RNA"#marker基因的排序还是以细胞为一组,主要是为了后续好分组markers <- c("SPARCL1","GNG11","IFITM3","AQP1",#macrophage             "CD3D","CD3E","CD3G","CD2",#T cell             "AIF1","MS4A6A","MNDA","LYZ",#DC             "S100A8","S100A9","ANXA1","CALML3",#Neutrophil             "CPA3","SLC18A2","RGS13","IL1RL1")#Mast#默认,order = TDotplot_anno(human_data,              features = markers,              celltype_color = dittoColors(),             group = c(rep('MAC',4), rep('Tcell',4),rep('DC',4),                       rep('NEU',4), rep('MAST',4)),             color = colorRampPalette(c("navy","white","firebrick3"))(100),             order = T)

图片

可以看到,顶部分组排序是按照首字母排序的,数据不是从到大到小从左到右排列。我们可以将分组排序,让其调整顺序。需要引入的一个参数是groups。​​​​​​​

#可以看到,顶部分组排序是按照首字母排序的#数据不是从到大到小从左到右排列。我们可以将分组排序,让其调整顺序#需要引入的一个参数是groupsDotplot_anno(human_data,              features = markers,              celltype_color = dittoColors(),             group = c(rep('MAC',4), rep('Tcell',4),rep('DC',4),                       rep('NEU',4), rep('MAST',4)),             color = colorRampPalette(c("navy","white","firebrick3"))(100),             order = T,             groups = c('MAC','MAST','DC','NEU','Tcell'))

图片

假设我们对于左侧celltype的排序也要按照自己顺序调整。那么order参数选择F,用level设置自己需要的顺序,groups也是可以设置的。

#假设我们对于左侧celltype的排序也要按照自己顺寻调整#那么order参数选择F,用level设置自己需要的顺序,groups也是可以设置的Dotplot_anno(human_data,              features = markers,              celltype_color = dittoColors(),             group = c(rep('MAC',4), rep('Tcell',4),rep('DC',4),                       rep('NEU',4), rep('MAST',4)),             color = colorRampPalette(c("navy","white","firebrick3"))(100),             order = F,             groups = c('MAC','MAST','DC','NEU','Tcell'),             level = c("T cell","Neutrophil","mDC","Mast","Macrophage"))

图片

其他的功能使用没有变化,函数的使用说明见B站视频,链接:

https://www.bilibili.com/video/BV16h411g7C1/?spm_id_from=333.999.0.0&vd_source=05b5479545ba945a8f5d7b2e7160ea34

这就是更新内容了,觉得分享有用的,点个赞再走呗!

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/3415.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

ElementUI —— Upload 文件上传

前言&#xff1a; 实现单文件上传&#xff01;此文章仅记录代码如何实现&#xff0c;便于后期项目学习使用&#xff01;&#xff01; 代码实现&#xff1a; <!--********** View ********** --> <el-uploadclass"upload-demo":http-request"UploadFi…

通过concurrently 包,在一个终端中通过npm script运行两个项目

npm script可以通过 && 串行&#xff0c;& 并行执行两个任务&#xff0c;但是不能在同个打开的终端中运行两个项目&#xff0c;可以通过concurrently包解决 安装 npm i concurrently -Dscript命令 "dev":"concurrently \"npm run start\&quo…

Numpy—ndarray对象

NumPy 最重要的一个特点是其N维数组对象 ndarray&#xff0c;它是一系列同类型数据的集合&#xff0c;以 0 下标为开始进行集合中元素的索引。 ndarray 对象是用于存放同类型元素的多维数组。ndarray 中的每个元素在内存中都有相同存储大小的区域。ndarray 内部由以下内容组成…

【分布式训练】基于Pytorch的分布式数据并行训练

基于Pytorch的分布式数据并行训练 动机为什么要并行分布数据&#xff1f;现有资料的不足 Outline整体框架图带解释的最小demo示例没有multiprocessing开启multiprocessing 分布式训练启动方式混合精度训练&#xff08;采用apex&#xff09; 参考资料 简介&#xff1a; 在PyTorc…

【ROS】ROS1人机界面开发:第一个最简ROS+QtGui程序(按钮启动发布者)

【ROS】郭老二博文之:ROS目录 1、创建工程 1)新建工程:Other Project --> ROS Workspace 2)设置工程名称、路径 3)可以通过点击“Browse”来创建目录 注意:使用自带ros插件的qtcreator-ros,无法创建目录、也不能选择目录,这是个bug,因此需要在终端手动创建目录…

密码学学习笔记(十六 ):Schnorr签名算法

交互式零知识证明 零知识证明(ZKP)就是不会将证据泄露给验证者的知识证明。Schnorr身份认证识别协议是一个交互式ZKP&#xff0c;它满足了完备性、可靠性、零知识性。所谓的交互式ZKP方案通常包含3个步骤&#xff08;承诺、挑战和证明&#xff09;&#xff0c;在文献中通常被称…

elastic elasticsearch 源码解析之选主选举过程

选主 选举算法 角色定义 这里的选主为什么提角色? 是因为不同角色在选主中起到不同的作用.master的非voting_only节点不但参与投票同时还可以参与竞选, master 的voting_only角色仅投票不参与竞选,其余角色不参与. 支持的角色 master data data_content data_hot data_wa…

django框架向DRF框架演变过程详解

一、Django框架实现项目查询接口 主要知识点&#xff1a; Django框架视图函数 1、在 Django 项目中创建一个应用&#xff08;如果还没有创建&#xff09;&#xff1a; python manage.py startapp projects 2、在项目的 models.py 文件中定义项目模型 from django.db impor…

React面试题汇总

1、面试官&#xff1a;说说对 React 的理解&#xff1f;有哪些特性&#xff1f; React遵循组件设计模式&#xff0c;使用虚拟 DOM 来有效地操作 DOM&#xff0c;遵循从高阶组件到低阶组件的单向数据流。React 特性有很多&#xff0c;如&#xff1a;JSX 语法、单向数据绑定、虚…

【C++】-stack和queue的具体使用以及模拟实现(dqeue的介绍+容器适配器的介绍)

&#x1f496;作者&#xff1a;小树苗渴望变成参天大树&#x1f388; &#x1f389;作者宣言&#xff1a;认真写好每一篇博客&#x1f4a4; &#x1f38a;作者gitee:gitee✨ &#x1f49e;作者专栏&#xff1a;C语言,数据结构初阶,Linux,C 动态规划算法&#x1f384; 如 果 你 …

Java八股文总结(二)

Java八股文总结&#xff08;续&#xff09; 接上篇笔记&#xff1a;Jhttps://blog.csdn.net/weixin_44780078/article/details/130192373 文章目录 Java八股文总结&#xff08;续&#xff09;六、MySql 相关1. InnoDB 与 MyISAM 的区别&#xff1f;2. 为什么 InnoDB 存储引擎表…

kubesphere部署kafka单点

zookeeper 存储: #zk-cm nfsRWO /bitnami/zookeeper 配置: 秘钥: 服务: 有状态服务 名称: zookeeper 镜像: bitnami/zookeeper:3.8.0 副本数量: 环境变量: - TZAsia/Shanghai - ZOO_SERVER_ID1 - ZOO_PORT_NUMBER2181 - ZOO_HEAP_SIZE128 - ALLOW_ANONYMOUS_LOGINyes 启动…

JavaFx 用户界面控件2——ListView

1.列表显示ListView 下面是一个JavaFX ListView的示例代码和使用方法&#xff1a; public class ListViewExample extends Application {Overridepublic void start(Stage primaryStage) {// 创建一个可观察的列表&#xff0c;用于存储ListView中的数据ObservableList<Str…

【MongoDB】SpringBoot整合MongoDB

【MongoDB】SpringBoot整合MongoDB 文章目录 【MongoDB】SpringBoot整合MongoDB0. 准备工作1. 集合操作1.1 创建集合1.2 删除集合 2. 相关注解3. 文档操作3.1 添加文档3.2 批量添加文档3.3 查询文档3.3.1 查询所有文档3.3.2 根据id查询3.3.3 等值查询3.3.4 范围查询3.3.5 and查…

UML与SYSML的关系

UML与SysML的联系 UML&#xff08;统一建模语言&#xff09;和SysML&#xff08;系统建模语言&#xff09;是两种与建模相关的语言&#xff0c;它们之间存在联系和区别。 SysML的图分类如下图所示。 联系 SysML是基于UML的&#xff0c;它重用了UML 2的子集&#xff0c;并提…

wordpress后台的路径都是空白,进不了后台怎么办?

主题或插件冲突&#xff1a;某些主题或插件可能与其他已安装的主题或插件不兼容&#xff0c;导致登录页面显示为空白。您可以通过禁用所有插件和将主题更改为默认主题来查看是否解决了问题。如果解决了问题&#xff0c;则可以逐个重新启用插件以找出是哪个插件造成冲突。 PHP错…

CGT Asia嘉年华|3D细胞培养与类器官研发峰会10月广州召开

类器官指利用成体干细胞或多能干细胞进行体外三维&#xff08;3D&#xff09;培养而形成的具有一定空间结构的组织类似物&#xff0c;是近10年来干细胞领域发展最快的研究热点之一。2022年&#xff0c;FDA 通过现代化法案 2.0&#xff0c;批准全球首个完全基于“类器官芯片”研…

Win11中的Swapfile.sys

除了 pagefile.sys 和 hiberfil.sys 文件外&#xff0c;在系统根目录会多出一个 swapfile.sys 虚拟内存文件。Windows 10/8 系统为什么会同时使用 SWAP 交换文件和 Page 页面文件呢&#xff1f; 其实 swapfile.sys 文件目前只被用来交换 Universal App (其实就是Metro App)的个…

线程、并发、多线程并发

单线程并发 单线程并发&#xff0c;并非真正意义上的单线程。而是只有单一的用户线程。还包括数据库&#xff0c;socket等系统多线程。 单个用户线程&#xff1a;对于十万个用户同时访问服务器&#xff0c;有两种方式处理并发。 为每个用户开一个新的用户线程&#xff0c;每个…

12.postgresql--plsql

[ <<label>> ] [ DECLAREdeclarations ] BEGINstatements END [ label ];do $$ beginraise notice 无Label标签与声明部分; end $$; 或者 do $$ <<test>> beginraise notice 使用label标签; end test $$; 或者 do $$ <<test>> declarei in…