R Histograms
In this article, you’ll learn to use hist() function to create histograms in R programming with the help of numerous examples.
Histogram can be created using the hist()
function in R programming language. This function takes in a vector of values for which the histogram is plotted.
Let us use the builtin dataset airquality
which has Daily air quality measurements in New York, May to September 1973.
R documentation.
> str(airquality)
'data.frame': 153 obs. of 6 variables:
$ Ozone : int 41 36 12 18 NA 28 23 19 8 NA ...
$ Solar.R: int 190 118 149 313 NA NA 299 99 19 194 ...
$ Wind : num 7.4 8 12.6 11.5 14.3 14.9 8.6 13.8 20.1 8.6 ...
$ Temp : int 67 72 74 62 56 66 65 59 61 69 ...
$ Month : int 5 5 5 5 5 5 5 5 5 5 ...
$ Day : int 1 2 3 4 5 6 7 8 9 10 ...
We will use the temperature parameter which has 154 observations in degree Fahrenheit.
Example 1: Simple histogram
Temperature < airquality$Temp
hist(Temperature)
We can see above that there are 9 cells with equally spaced breaks. In this case, the height of a cell is equal to the number of observation falling in that cell.
We can pass in additional parameters to control the way our plot looks. You can read about them in the help section ?hist
.
Some of the frequently used ones are, main
to give the title, xlab
and ylab
to provide labels for the axes, xlim
and ylim
to provide range of the axes, col
to define color etc.
Additionally, with the argument freq=FALSE
we can get the probability distribution instead of the frequency.
Example 2: Histogram with added parameters
# histogram with added parameters
hist(Temperature,
main="Maximum daily temperature at La Guardia Airport",
xlab="Temperature in degrees Fahrenheit",
xlim=c(50,100),
col="darkmagenta",
freq=FALSE
)
Note that the y axis is labelled density instead of frequency. In this case, the total area of the histogram is equal to 1.
Return Value of hist()
The hist()
function returns a list with 6 components.
> h < hist(Temperature)
> h
$breaks
[1] 55 60 65 70 75 80 85 90 95 100
$counts
[1] 8 10 15 19 33 34 20 12 2
$density
[1] 0.010457516 0.013071895 0.019607843 0.024836601 0.043137255
[6] 0.044444444 0.026143791 0.015686275 0.002614379
$mids
[1] 57.5 62.5 67.5 72.5 77.5 82.5 87.5 92.5 97.5
$xname
[1] "Temperature"
$equidist
[1] TRUE
attr(,"class")
[1] "histogram"
We see that an object of class histogram
is returned which has:
breaks
places where the breaks occur,counts
the number of observations falling in that cell,density
the density of cells,mids
the midpoints of cells,xname
the x argument name andequidist
a logical value indicating if the breaks are equally spaced or not.
We can use these values for further processing.
For example, in the following example we use the return values to place the counts on top of each cell using the text()
function.
Example 3: Use Histogram return values for labels using text()
h < hist(Temperature,ylim=c(0,40))
text(h$mids,h$counts,labels=h$counts, adj=c(0.5, 0.5))
Defining the Number of Breaks
With the breaks
argument we can specify the number of cells we want in the histogram. However, this number is just a suggestion.
R calculates the best number of cells, keeping this suggestion in mind. Following are two histograms on the same data with different number of cells.
Example 4: Histogram with different breaks
hist(Temperature, breaks=4, main="With breaks=4")
hist(Temperature, breaks=20, main="With breaks=20")
In the above figure we see that the actual number of cells plotted is greater than we had specified.
We can also define breakpoints between the cells as a vector. This makes it possible to plot a histogram with unequal intervals. In such case, the area of the cell is proportional to the number of observations falling inside that cell.
Example 5: Histogram with nonuniform width
hist(Temperature,
main="Maximum daily temperature at La Guardia Airport",
xlab="Temperature in degrees Fahrenheit",
xlim=c(50,100),
col="chocolate",
border="brown",
breaks=c(55,60,70,75,80,100)
)

PREVIOUS
R Bar Plot 
NEXT
R Pie Chart