今天要学做的是一篇发表在《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科研分享与服务》！
网友评论