Pan-African genome demonstrates how population-specific genome graphs improve high-throughput sequencing data analysis
https://www.nature.com/articles/s41467-022-31724-3
本地pdf s41467-022-31724-3.pdf
论文中公布了大部分图的数据，但是没有公布对应的作图代码，没有关系，我们可以自己写代码试着模仿，今天的推文重复一下论文中的Figure 2D 的下半部分
部分示例数据截图
library(readxl)
dat01<-read_excel("data/20220806/fig2d.xlsx",
sheet = "Sheet1")
library(tidyverse)
library(stringr)
dat01 %>%
pivot_longer(-var) %>%
mutate(value=replace_na(value,"B")) -> new.dat01
library(ggplot2)
ggplot(data=new.dat01,aes(x=name,y=var))+
geom_point(size=15,aes(color=value),
show.legend = FALSE)+
geom_segment(data=data.frame(x=7:11,
xend=7:11,
y=6,
yend=5:1),
aes(x=x,y=y,xend=xend,yend=yend),
color="#68c3a6",size=3)+
scale_color_manual(values = c("A"="#68c3a6",
"B"="#d3d3d3"))+
geom_rect(data = data.frame(xmin=0.5,
xmax=11.5,
ymin=c(0.5,2.5,4.5),
ymax=c(1.5,3.5,5.5)),
aes(xmin=xmin,xmax=xmax,ymin=ymin,ymax=ymax),
fill="#f6f6f6",
alpha=0.5,
inherit.aes = FALSE)+
geom_point(data=new.dat01,aes(x=name,y=var,color=value),
size=15,
show.legend = FALSE)+
geom_segment(data=data.frame(x=7:11,
xend=7:11,
y=6,
yend=5:1),
aes(x=x,y=y,xend=xend,yend=yend),
color="#68c3a6",size=3)+
theme_void() -> p1
p1
这里遇到了一个问题是 想把条状的灰色背景放到最底层，正常是把作图函数写到最前面就可以了，但是因为添加条状背景是连续数据，气泡图是离散数据，写到前面会报错，写到后面就没有这个问题，暂时不知道什么原因。这里的解决办法是有叠加了一层气泡图和线段的函数
部分示例数据集截图
dat02<-read_excel("data/20220806/fig2d.xlsx",
sheet = "Sheet2")
dat02
ggplot(data=dat02,aes(x=var,y=value))+
geom_hline(yintercept = 10,lty="dashed",
color="#d3d3d3")+
geom_col(fill="#68c3a6")+
geom_text(aes(label=paste0(value,"M")),hjust=1)+
coord_flip()+
scale_y_reverse(expand=expansion(mult = c(0.4,0)),
breaks=c(0,10),
labels=c(0,"10M"))+
scale_x_discrete(position = "top")+
theme_classic()+
theme(axis.line.y = element_blank(),
axis.ticks.y = element_blank())+
labs(x=NULL,y="Total size") -> p2
p2
library(patchwork)
p2 + p1 +
plot_layout(widths = c(1,3))
背景的浅色灰色条不知道如何跨越两个图，这里是出图后编辑图片实现的
示例数据和代码可以 给公众号推文打赏一元获取
欢迎大家关注我的公众号
小明的数据分析笔记本
小明的数据分析笔记本 公众号 主要分享：1、R语言和python做数据分析和数据可视化的简单小例子；2、园艺植物相关转录组学、基因组学、群体遗传学文献阅读笔记；3、生物信息学入门学习资料及自己的学习笔记！
我正在寻找执行以下操作的正确语法(在Perl、Shell或Ruby中):#variabletoaccessthedatalinesappendedasafileEND_OF_SCRIPT_MARKERrawdatastartshereanditcontinues. 最佳答案 Perl用__DATA__做这个:#!/usr/bin/perlusestrict;usewarnings;while(){print;}__DATA__Texttoprintgoeshere 关于ruby-如何将脚
几个月前，我读了一篇关于rubygem的博客文章，它可以通过阅读代码本身来确定编程语言。对于我的生活，我不记得博客或gem的名称。谷歌搜索“ruby编程语言猜测”及其变体也无济于事。有人碰巧知道相关gem的名称吗？ 最佳答案 是这个吗:http://github.com/chrislo/sourceclassifier/tree/master 关于ruby-寻找通过阅读代码确定编程语言的rubygem？，我们在StackOverflow上找到一个类似的问题：
?博客主页：https://xiaoy.blog.csdn.net?本文由呆呆敲代码的小Y原创，首发于CSDN??学习专栏推荐：Unity系统学习专栏?游戏制作专栏推荐：游戏制作?Unity实战100例专栏推荐：Unity实战100例教程?欢迎点赞?收藏⭐留言?如有错误敬请指正！?未来很长，值得我们全力奔赴更美好的生活✨------------------❤️分割线❤️-------------------------
嗨~大家好，这里是可莉！今天给大家带来的是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.创建临时变量来
关闭。这个问题是off-topic.它目前不接受答案。想改进这个问题吗？Updatethequestion所以它是on-topic用于堆栈溢出。关闭11年前。Improvethisquestion我不经常使用ruby-通常它加起来相当于每两个月或更长时间编写一次脚本。我的大部分编程都是使用C++进行的，这与ruby有很大不同。由于我与ruby之间的差距如此之大，我总是忘记语言的基本方面(比如解析文本文件和其他简单的东西)。我想每天练习一些基本的东西，我想知道是否有一些我可以订阅的网站，并且会向我发送当天的Ruby问题或类似的东西。有人知道这样的站点/Internet服务吗？
如果特定语言环境中缺少翻译，如何配置i18n以使用en语言环境翻译？当前已插入翻译缺失消息。我正在使用RoR3.1。 最佳答案 找到相似的question这里是答案:#application.rb#railswillfallbacktoconfig.i18n.default_localetranslationconfig.i18n.fallbacks=true#railswillfallbacktoen,nomatterwhatissetasconfig.i18n.default_localeconfig.i18n.fallback
在我的双语Rails4应用程序中，我有一个像这样的LocalesController:classLocalesController用户可以通过此表单更改其语言环境:deflocale_switcherform_tagurl_for(:controller=>'locales',:action=>'change_locale'),:method=>'get',:id=>'locale_switcher'doselect_tag'set_locale',options_for_select(LANGUAGES,I18n.locale.to_s)end这有效。但是，目前用户无法通过URL更改
我使用Ruby编程已经有一段时间了，现在只使用Ruby的标准MRI实现，但我一直对我经常听到的其他实现感到好奇。前几天我在读有关Rubinius的文章，这是一个用Ruby编写的Ruby解释器。我试着在不同的地方查找它，但我很难弄清楚这样的东西到底是如何工作的。我在编译器或语言编写方面从来没有太多经验，但我真的很想弄明白。一门语言究竟如何才能被自己解释？编译中是否有一个我不明白这有意义的基本步骤？有人可以像我是个白痴一样向我解释这个吗(因为无论如何这都不会太离谱) 最佳答案 它比你想象的要简单。Rubinius并非100%用Ruby编
Ruby是完全面向对象的语言。在ruby中，一切都是对象，因此属于某个类。例如5属于Objectclass1.9.3p194:001>5.class=>Fixnum1.9.3p194:002>5.class.superclass=>Integer1.9.3p194:003>5.class.superclass.superclass=>Numeric1.9.3p194:005>5.class.superclass.superclass.superclass=>Object1.9.3p194:006>5.class.superclass.superclass.superclass.su
我有一个Highstock图表(带有标记和阴影的线条)，并且想以编程方式显示一个highstock工具提示，例如，当我选择某个表上的一行(包含图表数据)我想显示相应的highstock工具提示。这可能吗？ 最佳答案 股票图表thissolution不起作用:在thisexample你必须更换这个:chart.tooltip.refresh(chart.series[0].data[i]);为此:chart.tooltip.refresh([chart.series[0].points[i]]);解决方案可用here.