今天要学做的是一篇发表在《nature genetic》上的图片,说是学做,其实我们只是利用了它的形式。折线图很简单,不过通过这幅图我们要掌握的是,第一:不单独计算mean、sd来实现折线图添加误差线。第二:利用循环实现批量作图,减轻工作负担。
图片来源
(Reference:Phased diploid genome assemblies and pan-genomes provide insights into the genetic history of apple domestication)
image.png读入数据
数据类似于基因在不同时期或者处理的表达量变化,每个组至少有3个重复,先对数据进行处理,转化为ggplot格式。
setwd("E:/生物信息学/批量折线图")
lc <- read.csv("批量折线图.csv", header=T, row.names=1)
lc$group <- as.factor(lc$group)
library(reshape2)
lc_plot <- melt(lc)
更多精彩内容请关注我的公众号《KS科研分享与服务》!
网友评论