前面说过了小提琴图的绘制(单细胞基因可视化之小提琴图)。这里继续说单细胞基因可视化的方法---气泡图,这是最受欢迎的一种展示基因表达的图了,但凡是单细胞的文章都会出现。这里我们提供的方法不仅仅局限于单细胞数据,其他的数据也可以用来做气泡热图。

Seurat包自带做气泡图的函数DotPlot,可以通过ggplot修饰,对图形的外形,颜色等修改,具体操作如下。

marker <- c("ACKR1","RAMP2","LUM","COL3A1","KRT14","KRT5","CD69","CD52")
DotPlot(scedata, features = marker)+coord_flip()+theme_bw()+theme(panel.grid = element_blank(), axis.text.x=element_text(hjust = 1,vjust=0.5))+labs(x=NULL,y=NULL)+guides(size=guide_legend(order=3))+scale_color_gradientn(values = seq(0,1,0.2),colours = c('#330066','#336699','#66CC66','#FFCC33'))

以上是最普通的作图,其他文章中还出现过更加高级的做法,例如这样的:

这样的:

还有这样的:

这里我们提供三种方法绘制单细胞基因表达气泡图,每种图形都可以经过尝试变化,形成新的图形,也可以展示不一样的东西!

更多精彩请关注我的公众号---《KS科研分享与服务》

单细胞基因个性化作图之气泡图相关推荐

  1. 单细胞基因可视化之热图改造修饰1

    热图不再过多介绍了,参考之前的内容(热图系列大全).单细胞基因可视化中热图也是比较受欢迎的,在分析完每群的marker基因之后,可以挑选显著的gene用seurat自带函数DoHeatmap可视化.当 ...

  2. 单细胞基因可视化之UMAP图修饰

    更多精彩请至<KS科研分享与服务公众号> 除了之前说过的三种常见单细胞基因可视化方法外,还有一种最常用的就是直接在UMAP或者TSNE降维图上显示表达某基因的细胞,这种方式更加直观,但是只 ...

  3. 评分addmodule 绝对值评分 8种方法可视化你的单细胞基因集打分gsea 缺氧评分

    基因集合打分不考虑基因高低表达 基因集合评分不考虑高低表达 绝对值 基因集合里面的正负表达影响addmodulescore的结果 构建的基因集进行评分时 基因的高低正负表达基因集合评分不考虑高低表达_ ...

  4. 一个R包完成单细胞基因集富集分析 (全代码)

    singleseqgset | 单细胞RNA-Seq基因集富集分析 NGS系列文章包括NGS基础.转录组分析 (Nature重磅综述|关于RNA-seq你想知道的全在这).ChIP-seq分析 (Ch ...

  5. 单细胞基因可视化之小提琴图

    完成单细胞定群后,就可以做很多东西了,首先可视化marker基因,这是很多单细胞文章中不可或缺的一部分,当然,marker基因的可视化可以拓展到对任意自己感兴趣基因的可视化. 小提琴图是其中的一种.S ...

  6. 【AUCell】单细胞基因集打分

    单细胞层次给基因集打分除了和GSVA等分析一样可以做富集以外,还可以通过cell marker gene set做亚群注释(需要注意细胞比例的问题,数量少的群体不行,可作辅助).前面写过的scMeta ...

  7. 跟着Cell学单细胞转录组分析(十一):单细胞基因评分|AUCell评分

    更多精彩请至我的公众号<KS科研分享与服务> 第一次接触基因评分是在一篇文章中,也不知道这样的叫法对不对,作者选定了几个炎症基因,利用seurat包的一个打分函数AddModuleScor ...

  8. 跟着Cell学单细胞转录组分析(八):单细胞转录组差异基因分析及多组结果可视化

    接着单细胞下游分析: 从Cell学单细胞转录组分析(一):开端!!! 跟着Cell学单细胞转录组分析(二):单细胞转录组测序文件的读入及Seurat对象构建 跟着Cell学单细胞转录组分析(三):单细 ...

  9. R统计和作图文章汇总

    R统计和作图 这个为生信学习和生信作图打造的开源R教程真香!!! 文章用图的修改和排版 (1) 文章用图的修改和排版 (2) ggplot2高效实用指南 (可视化脚本.工具.套路.配色) 在B站学习大 ...

最新文章

  1. axios不发起请求_重复的ajax请求让人很受伤
  2. python next用法
  3. 开源-一步步实现cnblogs博客采集工具-对项目需求的分析
  4. fft 重叠加法_FFT重叠相加法 参数如何定 h M?
  5. boost::process::ipstream相关的测试程序
  6. 将Sublime Text3添加到右键菜单中
  7. 演练 实现等腰三角形
  8. python调用aws接口添加安全组策略
  9. 转: android编译过程(流程图)
  10. 2020.06.25 端午节快乐
  11. IP库GeoLite2-City.mmdb的使用
  12. 新知实验室基于腾讯云实时音视频WEB端实验
  13. [转载]关于太阳(卫星)天顶角,太阳高度角,太阳方位角的整理_akala啦_新浪博客...
  14. RHEL7虚拟机实验快照
  15. 设为首页 加入收藏 html,JS设为首页和加入收藏的代码
  16. Opencv实现颜色检测
  17. 一文通俗搞懂线性无关特征向量个数≤特征值重数
  18. 【爬虫】牛客网简易爬虫(二)用Selenium看看你投过简历的公司这个月哪些有笔试+笔试时间
  19. PUSH消息是什么?
  20. java中文句号转换英文句号_java实现中文或其他语言及标点符号等转换成unicode字符串,或unicode的16进制码转换回文字或符号等...

热门文章

  1. 软件测试(概念Ⅱ) · 开发模型 · 软件的生命周期 · 瀑布模型 · 螺旋模型 · 增量模型 · 迭代模型 · 敏捷模型 · scrum · 软件测试模型之 V 模型 W 模型
  2. SQL中的空值、NULL以及类型转换那些坑
  3. centos7不显示ens33网卡的解决办法
  4. Mycat 分片表全局自增主键实现及测试
  5. 有序无序的变换 序号变图片 序号的各种形式变换
  6. 雅虎新CEO梅耶尔教你如何创新
  7. 通达OA工作流主要表的数据结构
  8. Unity --- 虚拟轴的使用
  9. [1688关键词卡首屏]1688关键词卡首屏原理,1688任意商品卡首屏接口,1688关键词卡首屏生成器,1688搜索卡首屏接口,1688搜索排名靠前,1688关键词推送商品API,1688卡首屏淘口
  10. 2023年04月个人工作生活总结