直方图绘制hist()

素娥

2020/9/27

本文我们要用iris数据集进行直方图的绘制

本文主要使用hist、plot函数进行绘制图片

首先还是观察数据类型

data<-iris
library(knitr)
tt<-head(data,15)
kable(tt)

简单的直方图

相关参数

1.x:要画图的数据

2.breaks:数据分多少组

3.freq:若为T（默认）则绘制频数分布直方图，为F则绘制频率分布直方图

4.probability:与freq相反

我们绘制iris的叶片宽度的直方图

绘制分12组的频数分布直方图

hist(x=data$Sepal.Width,breaks = 12)

绘制分12组的频率分布直方图

x<-hist(x=data$Sepal.Width,breaks = 12,freq = F)

注意，此处的纵坐标是频率/组距

再叠加一条密度曲线。

hist(x=data$Sepal.Width,breaks = 12,freq = F)
lines(density(data$Sepal.Width))

#美化图像(图形颜色，打斜线等)

###相关参数

1.density :阴影线密度

2.angle:阴影线倾角

3.col:填充颜色

4.border:直方图边缘颜色

在第一幅图的基础上进行美化

先加斜线

hist(x=data$Sepal.Width,breaks = 12,density = 10,angle = 45)

在加上颜色和边框

hist(x=data$Sepal.Width,breaks = 12,density = 10,angle = 45,col = "red",
border = "blue")

不加斜线，直接写颜色

hist(x=data$Sepal.Width,breaks = 12,col = "red",
border = "blue")

美化图像(坐标轴,标题,标签)

相关参数

1.main, xlab, ylab :标题，x名称，y名称

2.xlim, ylim :x轴，y轴范围

3.axes：逻辑值，是否画坐标轴

4.labels:是否添加标签

先加上标题，x名称，y名称

hist(x=data$Sepal.Width,breaks = 12,main="直方图",xlab="叶片宽度",
ylab="频数")

设置x,y轴的范围

hist(x=data$Sepal.Width,breaks = 12,main="直方图",xlab="叶片宽度",
ylab="频数",xlim=c(2,4.5),ylim=c(0,40))

加上数据标签

hist(x=data$Sepal.Width,breaks = 12,main="直方图",xlab="叶片宽度",
ylab="频数",xlim=c(2,4.5),ylim=c(0,40),labels=T)

去掉坐标轴

hist(x=data$Sepal.Width,breaks = 12,main="直方图",axes=F,labels=T)

参数plot

是否画出图像

hh<-hist(x=data$Sepal.Width,breaks = 12,plot=F)
hh

## $breaks
## [1] 2.0 2.2 2.4 2.6 2.8 3.0 3.2 3.4 3.6 3.8 4.0 4.2 4.4
##
## $counts
## [1] 4 7 13 23 36 24 18 10 9 3 2 1
##
## $density
## [1] 0.13333333 0.23333333 0.43333333 0.76666667 1.20000000 0.80000000
## [7] 0.60000000 0.33333333 0.30000000 0.10000000 0.06666667 0.03333333
##
## $mids
## [1] 2.1 2.3 2.5 2.7 2.9 3.1 3.3 3.5 3.7 3.9 4.1 4.3
##
## $xname
## [1] "data$Sepal.Width"
##
## $equidist
## [1] TRUE
##
## attr(,"class")
## [1] "histogram"

hh展示的是直方图的相关值

间隔，频数，频率/组距，组中值等

接着用plot函数就可以继续绘制

hh<-hist(x=data$Sepal.Width,breaks = 12,plot=F)
plot(hh,col=rgb(0.7,0.2,0.2,alpha=0.5),border = "white",
main="直方图",xlab="叶片宽度",
ylab="频数",xlim=c(2,4.5),ylim=c(0,40),labels=T,
)

再往上加一条密度图

plot(hh,
freq=F,col=rgb(0.7,0.2,0.2,alpha=0.5),border = "white",
main="直方图",xlab="叶片宽度",
ylab="频数",xlim=c(2,4.5),ylim=c(0,1.5),labels=T,
)
lines(density(data$Sepal.Width),lwd=2,col=rgb(0.1,0.5,0.2))