jjzjj

跟着BMC Plant Biology学作图：R语言ggtree圆形树形图并添加分组背景色

小明的数据分析笔记本 2024-01-12 原文

论文

Comparative analysis of de novo genomes reveals dynamic intra‑species divergence of NLRs in pepper

数据和代码

https://github.com/sdaf11111/NLR-map-in-pepper

论文中Figure2的示例数据和代码作者放到了github主页，我们可以学习一下他的代码

示例数据是一个nwk格式的树文件 和 一个csv格式的分组文件

学到的新知识点

R包 svglite 输出图片如果保存为svg格式可能会用到这个R包

函数split()可以把数据框根据某一列分组转换成列表格式，文字表达可能有点看不明白，看一下函数的输出效果

df<-data.frame(x=c("A","A","B","B","B"),
               y=c(1,4,5,2,7))
df
split(df$y,df$x)
image.png

函数geom_abline() ggtree里的函数，可以把树末端补平

接下来是实际的代码

首先是加载需要用到的R包

library(dplyr)
library(ggtree)
library(ggplot2)
library(svglite)
library(scales)

读取数据

info <- read.csv("data/20230202/Annuum.Intact.NBARC.group.trimal92.csv")
head(info)
tree <- read.tree("data/20230202/Annuum.Intact.NBARC.tree.trimal92.nwk.txt")

给树添加分组信息

groupInfo<-split(info$ID,info$Group)
tree<- groupOTU(tree, groupInfo)

指定颜色

heatmap.colours <- c("#be9fe1","#8ac6d1","#e1ccec","#fddb3a",
                     "#C0C0C0","#c9b6e4","#d5c455","#ffb6b9","#fae3d9",
                     "#9aceff","#d7cde6","#bbded6","#ede59a","#4f98ca",
                     "#4a69bb","#f5cdaa","NA",
                     "#c3d14a","#63b637","#FF0000","#008000","#FF1493","#FF4500")
names(heatmap.colours) <- c("G1","G2","G3","G4",
                            "G5","G6","G7","G8","G9",
                            "G10","G11","G12","G13","GT",
                            "GR","G14","NG",
                            "CHIL","Know","CANN","CECW","CZUN","CASF")

作图代码

p <- ggtree(tree, layout='circular', size=0.2) %<+% info +
  geom_aline(linetype="solid", 
             size=0.5, aes(color=group),
             alpha=0.5) +   
  scale_colour_manual(values=heatmap.colours,
                      breaks=c("G1","G2","G3","G4","G5",
                               "G6","G7","G8","G9","G10",
                               "G11","G12","G13","GT","GR",
                               "G14","CHIL","Know","CANN",
                               "CECW","CZUN","CASF"), 
                      name="Group") +
  theme(legend.position="right")+
  geom_tippoint(aes(color=Species), size=0.2)+
  guides(color=guide_legend(ncol=6))
q <- flip(p, 3001, 4090) %>% rotate(3507)

ggsave(file="NLR_tree.92.pdf", plot=q,
       width = 9.4,height = 4)
image.png

这里有一个问题是geom_abline()函数的效果如果是在Rstudio的图片显示界面是看不到的，如果保存为pdf就可以看到效果。暂时不明白是什么原因

论文中提供的代码还有一部分是计算MRCA，这部分我暂时没有想明白，想明白了再来介绍吧

示例数据和代码可以到论文中提到的链接处下载，或者给推文点赞，点击在看，最后留言获取

欢迎大家关注我的公众号

小明的数据分析笔记本

小明的数据分析笔记本 公众号 主要分享：1、R语言和python做数据分析和数据可视化的简单小例子；2、园艺植物相关转录组学、基因组学、群体遗传学文献阅读笔记；3、生物信息学入门学习资料及自己的学习笔记！

作图树形sectioncodeblockquote

有关跟着BMC Plant Biology学作图：R语言ggtree圆形树形图并添加分组背景色的更多相关文章

  1. ruby-on-rails - 按天对 Mongoid 对象进行分组 - 2

    在控制台中反复尝试之后，我想到了这种方法，可以按发生日期对类似activerecord的(Mongoid)对象进行分组。我不确定这是完成此任务的最佳方法，但它确实有效。有没有人有更好的建议，或者这是一个很好的方法？#eventsisanarrayofactiverecord-likeobjectsthatincludeatimeattributeevents.map{|event|#converteventsarrayintoanarrayofhasheswiththedayofthemonthandtheevent{:number=>event.time.day,:event=>ev

  2. ruby - 如何将脚本文件的末尾读取为数据文件(Perl 或任何其他语言) - 2

    我正在寻找执行以下操作的正确语法(在Perl、Shell或Ruby中):#variabletoaccessthedatalinesappendedasafileEND_OF_SCRIPT_MARKERrawdatastartshereanditcontinues. 最佳答案 Perl用__DATA__做这个:#!/usr/bin/perlusestrict;usewarnings;while(){print;}__DATA__Texttoprintgoeshere 关于ruby-如何将脚

  3. ruby - 寻找通过阅读代码确定编程语言的ruby gem？ - 2

    几个月前，我读了一篇关于ruby​​gem的博客文章，它可以通过阅读代码本身来确定编程语言。对于我的生活，我不记得博客或gem的名称。谷歌搜索“ruby编程语言猜测”及其变体也无济于事。有人碰巧知道相关gem的名称吗？ 最佳答案 是这个吗:http://github.com/chrislo/sourceclassifier/tree/master 关于ruby-寻找通过阅读代码确定编程语言的rubygem？，我们在StackOverflow上找到一个类似的问题：

  4. Unity 热更新技术 | （三） Lua语言基本介绍及下载安装 - 2

    ?博客主页：https://xiaoy.blog.csdn.net?本文由呆呆敲代码的小Y原创，首发于CSDN??学习专栏推荐：Unity系统学习专栏?游戏制作专栏推荐：游戏制作?Unity实战100例专栏推荐：Unity实战100例教程?欢迎点赞?收藏⭐留言?如有错误敬请指正！?未来很长，值得我们全力奔赴更美好的生活✨------------------❤️分割线❤️-------------------------

  5. 7个大一C语言必学的程序 / C语言经典代码大全 - 2

    嗨~大家好，这里是可莉！今天给大家带来的是7个C语言的经典基础代码~那一起往下看下去把【程序一】打印100到200之间的素数#includeintmain(){ inti; for(i=100;i 【程序二】输出乘法口诀表#includeintmain(){inti;for(i=1;i 【程序三】判断1000年---2000年之间的闰年#includeintmain(){intyear;for(year=1000;year 【程序四】给定两个整形变量的值，将两个值的内容进行交换。这里提供两种方法来进行交换，第一种为创建临时变量来进行交换，第二种是不创建临时变量而直接进行交换。1.创建临时变量来

  6. ruby - 在 Ruby 中创建按公共(public)键值分组的新哈希 - 2

    假设我有一个在Ruby中看起来像这样的哈希:{:ie0=>"Hi",:ex0=>"Hey",:eg0=>"Howdy",:ie1=>"Hello",:ex1=>"Greetings",:eg1=>"Goodday"}有什么好的方法可以将它变成如下内容:{"0"=>{"ie"=>"Hi","ex"=>"Hey","eg"=>"Howdy"},"1"=>{"ie"=>"Hello","ex"=>"Greetings","eg"=>"Goodday"}} 最佳答案 您要求一个好的方法来做到这一点，所以答案是:一种您或同事可以在六个月后理解

  7. arrays - 如何在下面的示例中将两个值数组分组为 n 个值数组？ - 2

    我已经有很多两个值数组，例如下面的例子ary=[[1,2],[2,3],[1,3],[4,5],[5,6],[4,7],[7,8],[4,8]]我想把它们分组到[1,2,3],[4,5],[5,6],[4,7,8]因为意思是1和2有关系，2和3有关系，1和3有关系，所以1,2,3都有关系我如何通过ruby​​库或任何算法来做到这一点？ 最佳答案 这是基本Bron–Kerboschalgorithm的Ruby实现:classGraphdefinitialize(edges)@edges=edgesenddeffind_maximum_

  8. ruby - 如何将相同的相邻数字分组 - 2

    如果至少有两个相邻的数字相同，格式为，我需要打包.这是我的输入:[2,2,2,3,4,3,3,2,4,4,5]以及预期的输出:"2:3,3,4,3:2,2,4:2,5"到目前为止我试过:a=[1,1,1,2,2,3,2,3,4,4,5]a.each_cons(2).any?do|s,t|ifs==t如果相等，也许可以尝试计数器，但那是行不通的。 最佳答案 您可以使用Enumerable#chunk_while(如果你使用的是Ruby>=2.3):a.chunk_while{|a,b|a==b}.flat_map{|chunk|chu

  9. ruby - 如何保持我不常用的编程语言技能 - 2

    关闭。这个问题是off-topic.它目前不接受答案。想改进这个问题吗？Updatethequestion所以它是on-topic用于堆栈溢出。关闭11年前。Improvethisquestion我不经常使用ruby​​-通常它加起来相当于每两个月或更长时间编写一次脚本。我的大部分编程都是使用C++进行的，这与ruby​​有很大不同。由于我与ruby​​之间的差距如此之大，我总是忘记语言的基本方面(比如解析文本文件和其他简单的东西)。我想每天练习一些基本的东西，我想知道是否有一些我可以订阅的网站，并且会向我发送当天的Ruby问题或类似的东西。有人知道这样的站点/Internet服务吗？

  10. ruby-on-rails - 如果特定语言环境中缺少翻译，如何配置 i18n 以使用 en 语言环境？ - 2

    如果特定语言环境中缺少翻译，如何配置i18n以使用en语言环境翻译？当前已插入翻译缺失消息。我正在使用RoR3.1。 最佳答案 找到相似的question这里是答案:#application.rb#railswillfallbacktoconfig.i18n.default_localetranslationconfig.i18n.fallbacks=true#railswillfallbacktoen,nomatterwhatissetasconfig.i18n.default_localeconfig.i18n.fallback

随机推荐