热点新闻
跟着Nature Communication学作图:R语言ggplot2话点线图展示基因表达量的范围
2023-11-27 06:04  浏览:2269  搜索引擎搜索“手机低淘网”
温馨提示:信息一旦丢失不一定找得到,请务必收藏信息以备急用!本站所有信息均是注册会员发布如遇到侵权请联系文章中的联系方式或客服删除!
联系我时,请说明是在手机低淘网看到的信息,谢谢。
展会发布 发布信息 广告合作 软文发布

论文

Microbiomes in the Challenger Deep slope and bottom-axis sediments

https://www.nature.com/articles/s41467-022-29144-4#code-availability

对应代码链接

https://github.com/ucassee/Challenger-Deep-Microbes

论文里提供了大部分图的数据和代码,很好的学习材料,感兴趣的同学可以找来参考,今天的推文重复一下论文中的Figure3b

示例数据集部分截图




image.png

读取数据

dat01<-read.delim("data/20220602/NC_Figure3b.txt", header=TRUE, check.names = FALSE, sep="\t", stringsAsFactors = TRUE) head(dat01) table(dat01$Group) dat01$gene<-factor(dat01$gene, levels = dat01$gene)

作图代码

library(ggplot2) library(scales) ggplot(data=dat01,aes(x=gene,y=mean_value,color=Group))+ geom_pointrange(aes(ymin=min_value,ymax=max_value), show.legend = FALSE)+ scale_y_continuous(trans = log10_trans() ,breaks = trans_breaks("log10", function(x) 10^x), labels = trans_format("log10", math_format(10^.x)))+ scale_color_manual(values=c("#BB3D00","#6F00D2","#D9006C", "#005AB5","#00DB00","#3D7878","#D9B300"))+ theme_bw() + theme(panel.background=element_rect(fill="white",colour="black",size=0.25), axis.line=element_line(colour="black",size=0.25), axis.title=element_text(size=13,face="plain",color="black"), axis.text = element_text(size=10,face="plain",color="black"), axis.text.x = element_text(angle = 60, hjust = 1, vjust =1), #legend.position="none" )+ labs(x=NULL,y="TMP of genes")


image.png

示例数据和代码可以在公众号后台留言20220610获取

欢迎大家关注我的公众号

小明的数据分析笔记本

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

发布人:c248****    IP:117.173.23.***     举报/删稿
展会推荐61905
让朕来说2句
评论
收藏
点赞
转发