原创 黄小仙
上次分享了小提琴曲线(violin plot)的作图方法,今天小仙同学给大家介绍一下如何用R画出漂亮的密度图(density plot)。
Step1. 绘图数据的准备
首先还是要把你想要绘图的数据调整成R语言可以识别的格式，把excel保存成csv格式。
数据的格式如下图：一列表示一种变量，第一行是列名
Step2. 绘图数据的读取
data<-read.csv(“your file path”, header = T)
#注释：header=T表示数据中的第一行是列名，如果没有列名就用 header=F
Step3.绘图所需package的安装、调用
install.package(“ggplot2”)
#注释：ggplot2是目前公认绘图很强的一个安装包，如果已经安装，这句就不需要啦
library(ggplot2)
#注释：package使用之前需要调用(每次用前都要调用哦)
Step4绘图
p<-ggplot(data, aes(x = weight))
#注释：x轴表示weight，y轴表示频率就不需要指定啦
p + geom_density(color = “black”, fill = “gray”)
#注释：color是线条的颜色参数，fill表示填充颜色
Step5.美化
p + geom_density(aes(color = sex))
#注释:按照性别不同组改变线条颜色
p + geom_density(aes(fill = sex), alpha=0.4)
#注释:按照性别不同组改变填充颜色，alpha表示调整透明度
到这里你已经可以画出比较高B格的density plot了，不过有些同学可能有
特殊的要求，比如把两组的平均值画出来，那本小仙再小秀一下?
1.准备平均值数据并读取
mean<-read.csv(“your file path”)
2.画图
p+ geom_density(aes(color = sex), alpha=0.4)+
geom_vline(data = mean, aes(xintercept = mean, color=sex),linetype=“dashed”)
#注释:前面是和之前一样作图语句，在后面直接加上这一句就好啦，xintercept表示x轴截断线，dashed表示是虚线
好啦，今天的分享就到这里了，希望可以帮助到大家。有什么新的作图要求，可以留言给我哦！
导出高清图的方法在这里：
R语言作图技巧——导出高清图
R语言作图系列还有：
R语言作图——Histogram
R语言作图——Violin plot
R语言作图——Facet violin plot
R语言作图——Line plot
R语言作图——Density plot
R语言作图——Bubble matrix
R语言作图——Radar chart
R语言作图——Pie chart
（公众号：生信了）