原创:黄小仙
为了画今天的这个图,小仙决定凭空想象一台可以实时监控基因表达水平的设备,成功得到了这么一组数据。
想要画的图是这样子滴。
Step1. 绘图数据的准备
首先要把你想要绘图的数据调整成R语言可以识别的格式,建议大家在excel中保存成csv格式。
Step2. 绘图数据的读取
代码语言:javascript复制data<-read.csv(“your file path”, header = T)
#注释:header=T表示数据中的第一行是列名,如果没有列名就用header=F
Step3. 绘图所需package的安装、调用
代码语言:javascript复制library(ggplot2)
# 注释:package使用之前需要调用
Step4. 绘图
代码语言:javascript复制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)
简单画出曲线直接用
代码语言:javascript复制ggplot(data,aes(x=Time, y=Expression,linetype=type))
geom_line()
今天画的图也比较简单,小仙就多啰嗦几句。有段时间小仙觉得R语言作图这个系列差不多要结束了,因为常见的图表都快画了一遍了,有点江郎才尽的感觉,不知道接下来还能画什么。好在经过一段时间的学习之后,又有了新的感悟,突然又觉得图是画不完的,即使是line plot也有好多花样。目前基本的图形小仙都能画出来,但距离随心所欲画出想要的效果,还是差十万八千里呐。学无止境,大家都加油呀!
(公众号:生信了)