ggtree绘制进化树

setwd("C:\\Users\\Administrator\\Desktop\\tree")
library("ggtree")
library("ggplot2")
tree <- read.tree("species.tree")
ggtree(tree,layout = "circular")

 

ggtree(tree, layout="fan", ladderize = FALSE, branch.length = "none")+
  geom_tiplab2(size=3)+ 
  theme(legend.position = "right")

 

ggtree(tree) + geom_text(aes(label=node), hjust=-.3) + geom_tiplab(size=3)

 

MRCA(tree, tip=c("Pelobacter_carbinolicus", "uncultured_nitrogen-fixing_bacterium"))
ggtree(tree) + 
  geom_cladelabel(node=40, label="Some random clade", color="red")

 

ggtree(tree) + 
  geom_tiplab() + 
  geom_cladelabel(node=5, label="Some random clade", 
                  color="red2", offset=-.3, align=TRUE) + 
  geom_cladelabel(node=30, label="A different clade", 
                  color="blue", offset=-.3, align=TRUE) + 
  theme_tree2() + 
#  xlim(0, 1) + 
  theme_tree()

 

ggtree(tree) + 
  geom_tiplab() + 
  geom_hilight(node=38, fill="gold") + 
  geom_hilight(node=17, fill="purple")

 

ggtree(tree) + 
  geom_tiplab() + 
  geom_taxalink("Pelobacter_carbinolicus", "uncultured_nitrogen-fixing_bacterium", color="blue3") +
  geom_taxalink("Dysgonomonas_gadei", "uncultured_low_G+C_Gram-positive_bacterium", color="orange2", curvature=-.9)

 

library("ggstance")
library("vegan")
abun<-read.table("otu_table_group.s.relative.mat.draw",header = T)
colnames(abun)

library("reshape2")

mydata<-melt(abun,id=colnames(abun)[1],measure=colnames(abun)[2:length(colnames(abun))],
             variable.name = "group",value.name = "value")

p <- ggtree(tree)+geom_tiplab(size=3) #+geom_text(aes(label=node), hjust=-.3,size=3) + xlim(0,2) 

p2 <- facet_plot(p, panel="Stacked Barplot", data=mydata, geom = geom_barh,offset = .15,scales="free",
                 mapping = aes(x = decostand(value ,"log"), fill = as.factor(group)), 
                 stat='identity' )+theme_tree2(legend.position="right") #+ theme(legend.position="right") #+theme(legend.position=c(.1, .8))#
p2

 

p <- ggtree(tree,layout = "fan",ladderize = TRUE)+geom_tiplab2(size=5, align = TRUE, offset = 0.5)
otu_table = read.delim("otu_table_group.s.relative.mat.draw", row.names= 1,  header=T, sep="\t")
otu_table<-log(otu_table+0.0000001)
p3<-gheatmap(p, otu_table, offset = -.010, width=0.5, font.size=3,low = "green",
             high = "red", color = "white", colnames_level=c("S","R","SN","H","K","KN","HN"), 
             hjust = 0)#, colnames_angle=-45, hjust=-.1)
p3

 

p <- ggtree(tree)+geom_tiplab(size=5, align = TRUE)
p
otu_table = read.delim("otu_table_group.s.relative.mat.draw", row.names= 1,  header=T, sep="\t")
otu_table<-log(otu_table+0.0000001)
p3<-gheatmap(p, otu_table, offset = 0.6, width=0.5, font.size=3,low = "green",
             high = "red", color = "white", colnames_level=c("S","R","SN","H","K","KN","HN"), 
             hjust = 0)#, colnames_angle=-45, hjust=-.1)
p3



作者:周运来就是我
链接:https://www.jianshu.com/p/6bfbae082cf5
来源:简书
著作权归作者所有。商业转载请联系作者获得授权,非商业转载请注明出处。

  • 7
    点赞
  • 40
    收藏
    觉得还不错? 一键收藏
  • 2
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值