R语言作图——Line plot with error

原创：黄小仙

为了画今天的这个图，小仙决定凭空想象一台可以实时监控基因表达水平的设备，成功得到了这么一组数据。

想要画的图是这样子滴。

Step1. 绘图数据的准备

首先要把你想要绘图的数据调整成R语言可以识别的格式，建议大家在excel中保存成csv格式。

Step2. 绘图数据的读取

data<-read.csv(“your file path”, header = T)
#注释：header=T表示数据中的第一行是列名，如果没有列名就用header=F

Step3. 绘图所需package的安装、调用

library(ggplot2)                
# 注释：package使用之前需要调用

Step4. 绘图

ggplot(data,aes(x=Time, y=Expression,linetype=type)) + 
  geom_line() +
  geom_ribbon(aes(ymin= Expression.min,ymax=Expression.max,fill=type),alpha=0.5)+
  xlab("Time") + 
  ylab("Expression") +
  ylim(0.8,1.1)

简单画出曲线直接用

ggplot(data,aes(x=Time, y=Expression,linetype=type)) + 
  geom_line() 

今天画的图也比较简单，小仙就多啰嗦几句。有段时间小仙觉得R语言作图这个系列差不多要结束了，因为常见的图表都快画了一遍了，有点江郎才尽的感觉，不知道接下来还能画什么。好在经过一段时间的学习之后，又有了新的感悟，突然又觉得图是画不完的，即使是line plot也有好多花样。目前基本的图形小仙都能画出来，但距离随心所欲画出想要的效果，还是差十万八千里呐。学无止境，大家都加油呀！

（公众号：生信了）