R语言作图——histogram(直方图)

原创:黄小仙

最近小仙同学很是烦恼,本以为自己已经掌握了ggplot2作图的语法,用read.csv(), ggplot()+geom_point()/boxplot()/violinplot()…就可以画遍天下图表,结果却发现到真正画图的时候,还是会出现不少的小问题
R语言作图——histogram(直方图)_第1张图片

比如小仙最近要画一个直方图,最开始用hist()函数试了一下,看了下形状, 好像因为数据取值范围跨度太大(最大值104,724,最小值30),这个图画出来有点丑,于是决定用ggplot美化一下。
R语言作图——histogram(直方图)_第2张图片调整之后好看是好看了,但是大家有没有看出什么不对的地方,明明bins=10但是只画出8个格子,之后调整bins的值,每次都会比我指定的值少2个格子。而且,图中第一个格子(取值范围0-1250)应该有700多个数据,但是图上显示只有不到300个,问题出在哪里呢?小仙同学百思不得其解。在geom_histogram()函数中,bins就是用来指定分组数目(格子),为什么总是会少两个?
R语言作图——histogram(直方图)_第3张图片

小仙同学考虑到自己能力有限,决定量力而行,另辟蹊径。

于是设置另外一个参数breaks,终于找回了丢失的格子
R语言作图——histogram(直方图)_第4张图片

经过此事,小仙同学深刻认识到了自己的有限水平,哈哈。不过还是分享一下,希望能帮助到大家。实际的数据可真是比书上的例子难处理呢。

按照惯例写一下整个作图的过程

Step1. 绘图数据的准备
首先要把你想要绘图的数据调整成R语言可以识别的格式,建议大家在excel中保存成csv格式。
R语言作图——histogram(直方图)_第5张图片

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

Step3.绘图所需package的安装、调用
library(ggplot2)
#注释:package使用之前需要调用

Step4绘图
p<-ggplot(data, aes(x=data$销量)) +
geom_histogram(breaks=seq(0,10000,1000))+ xlim(0,10000)
p
R语言作图——histogram(直方图)_第6张图片

(公众号:生信了)

你可能感兴趣的:(#,R语言,R语言,生物信息,绘图,直方图)