跟着 Cell Reports 学画图: Ribo-seq frame 与 reads length 分布

远方

1引言

Ribo-seq 的质量展示图包括 不同 reads 长度的 frame 分布, reads length 长度分布 等,一般都是分开展示的,如下:

偶然间看到 2019 年 Cell Reports 的一篇文章换一种方式来呈现,感觉还挺不错的,level 一下就提升了:

文章标题:

The ribosome collision sensor Hel2 functions as preventive quality control in the secretory pathway

今天来学习一下如何绘制这篇图。首先展示一下最终效果:

付费的小伙伴们记得私聊我,我把 文献,测试数据,代码分享给你。

2load data

library(ggplot2)
library(cowplot)

# load data
len <- read.csv('len.csv',header = T)
len.frame.sum <- read.csv('len.frame.sum.csv',header = T)

# 查看数据
head(len,3)
  pos   Freq
1  24 169937
2  25 192774
3  26 244397

head(len.frame.sum,3)
  aln.qwidth   frame count         per
1         24 Frame 0 81107 0.004560747
2         24 Frame 1 33477 0.001882453
3         24 Frame 2 55353 0.003112568

观察文章图片,发现图例标签的样式比较特别:

3绘图