Skip to main content
Statistics LibreTexts

4.2: 1-Dimensional Plots

  • Page ID
    3562
  • Our firm has just seven workers. How to analyze the bigger data? Let us first imagine that our hypothetical company prospers and hired one thousand new workers! We add them to our seven data points, with their salaries drawn randomly from interquartile range of the original sample (Figure \(\PageIndex{1}\)):

    CodeBox (R) \(\PageIndex{1}\): Boxplots

    salary <- c(21, 19, 27, 11, 102, 25, 21)
    new.1000 <- sample((median(salary) - IQR(salary)) : + (median(salary) + IQR(salary)), 1000, replace=TRUE)
    salary2 <- c(salary, new.1000)
    boxplot(salary2, log="y")
    

    In a code above we also see an example of data generation. Function sample() draws values randomly from a distribution or interval. Here we used replace=TRUE, since we needed to pick a lot of values from a much smaller sample. (The argument replace=FALSE might be needed for imitation of a card game, where each card may only be drawn from a deck once.) Please keep in mind that sampling is random and therefore each iteration will give slightly different results.

    Screen Shot 2019-01-09 at 3.09.45 PM.png

    Figure \(\PageIndex{1}\) The boxplot.

    Let us look at the plot. This is the boxplot (“box-and-whiskers” plot). Kathryn’s salary is the highest dot. It is so high, in fact, that we had to add the parameter log="y" to better visualize the rest of the values. The box (main rectangle) itself is bound by second and fourth quartiles, so that its height equals IQR. Thick line in the middle is a median. By default, the “whiskers” extend to the most extreme data point which is no more than 1.5 times the interquartile range from the box. Values that lay farther away are drawn as separate points and are considered outliers. The scheme (Figure \(\PageIndex{2}\)) might help in understanding boxplots.

    Screen Shot 2019-01-09 at 3.19.12 PM.png

    Figure \(\PageIndex{2}\) The structure of the boxplot (“box-and-whiskers” plot).

    Numbers which make the boxplot might be returned with fivenum() command. Boxplot representation was created by a famous American mathematician John W. Tukey as a quick, powerful and consistent way of reflecting main distribution-independent characteristics of the sample. In R, boxplot() is vectorized so we can draw several boxplots at once (Figure \(\PageIndex{3}\)):

    CodeBox (R) \(\PageIndex{2}\): Boxplots

    boxplot(scale(trees))
    

    (Parameters of trees were measured in different units, therefore we scale()’d them.)

    Histogram is another graphical representation of the sample where range is divided into intervals (bins), and consecutive bars are drawn with their height proportional to the count of values in each bin (Figure \(\PageIndex{4}\)):

    CodeBox (R) \(\PageIndex{3}\): Histograms

    hist(salary2, breaks=20)
    

    (By default, the command hist() would have divided the range into 10 bins, but here we needed 20 and therefore set them manually. Histogram is sometimes a rather cryptic way to display the data. Commands Histp() and Histr() from the asmisc.r will plot histograms together with percentages on the top of each bar, or overlaid with normal curve (or density—see below), respectively. Please try them yourself.)

    A numerical analog of a histogram is the function cut():

    CodeBox (R) \(\PageIndex{4}\)

    table(cut(salary2, 20))
    

    There are other graphical functions, conceptually similar to histograms. The first is stem-and-leaf plot. stem() is a kind of pseudograph, text histogram. Let us see how it treats the original vector salary:

    CodeBox (R) \(\PageIndex{5}\): stem-and-leaf plot

    stem(salary, scale=2)
    

    Screen Shot 2019-01-09 at 3.22.21 PM.png

    Figure \(\PageIndex{3}\) Three boxplots, each of them represents one column of the data.

    The bar | symbol is a “stem” of the graph. The numbers in front of it are leading digits of the raw values. As you remember, our original data ranged from 11 to 102—therefore we got leading digits from 1 to 10. Each number to the left comes from the next digit of a datum. When we have several values with identical leading digit, like 11 and 19, we place their last digits in a sequence, as “leafs”, to the left of the “stem”. As you see, there are two values between 10 and 20, five values between 20 and 30, etc. Aside from a histogram-like appearance, this function performs an efficient ordering.

    Screen Shot 2019-01-09 at 3.23.49 PM.png

    Figure \(\PageIndex{4}\) Histogram of the 1007 hypothetical employees’ salaries.

    Another univariate instrument requires more sophisticated calculations. It is a graph of distribution density, density plot (Figure \(\PageIndex{5}\)):

    CodeBox (R) \(\PageIndex{6}\): Density Plots

    plot(density(salary2, adjust=2))
    rug(salary2)
    

    (We used an additional graphic function rug() which supplies an existing plot with a “ruler” which marks areas of highest data density.)

    Here the histogram is smoothed, turned into a continuous function. The degree to which it is “rounded” depends on the parameter adjust. Aside from boxplots and a variety of histograms and alike, R and external packages provide many more instruments for univariate plotting.

    Screen Shot 2019-01-09 at 3.25.20 PM.png

    Figure \(\PageIndex{5}\) Distribution density of the 1007 hypothetical employees’ salaries.

    One of simplest is the stripchart. To make stripchart more interesting, we complicated it below using its ability to show individual data points:

    CodeBox (R) \(\PageIndex{7}\): stripchart

    trees.s <- data.frame(scale(trees), class=cut(trees$Girth, 3, labels=c("thin", "medium", "thick")))
    stripchart(trees.s[, 1:3], method="jitter", cex=2, pch=21, col=1:3, bg=as.numeric(trees.s$class))
    legend("right", legend=levels(trees.s$class), pch=19, pt.cex=2, col=1:3)
    

    (By default, stripchart is horizontal. We used method="jitter" to avoid overplotting, and also scaled all characters to make their distributions comparable. One of stripchart features is that col argument colorizes columns whereas bg argument (which works only for pch from 21 to 25) colorizes rows. We split trees into 3 classes of thickness, and applied these classes as dots background. Note that if data points are shown with multiple colors and/or multiple point types, the legend is always necessary.‘)

    Screen Shot 2019-01-09 at 3.28.12 PM.png

    Figure \(\PageIndex{6}\) Stripchart for modified trees data.

    Beeswarm plot requires the external package. It is similar to stripchart but has several advanced methods to disperse points, plus an ability to control the type of individual points (Figure \(\PageIndex{7}\)):

    CodeBox \(\PageIndex{8}\) (R): Beeswarm plots

    library(beeswarm)
    beeswarm(trees.s[, 1:3], cex=2, col=1:3, pwpch=rep(as.numeric(trees.s$class), 3))
    bxplot(trees.s[, 1:3], add=TRUE)
    legend("top", pch=1:3, legend=levels(trees.s$class))
    

    Screen Shot 2019-01-09 at 3.29.42 PM.png

    Figure \(\PageIndex{7}\) Beeswarm plot with boxplot lines.

    (Here with bxplot() command we added boxplot lines to a beehive graph in order to visualize quartiles and medians. To overlay, we used an argument add=TRUE.)

    And one more useful 1-dimensional plot. It is a similar to both boxplot and density plot (Figure \(\PageIndex{8}\)):

    CodeBox (R) \(\PageIndex{9}\): Bean plots

    library(beanplot)
    beanplot(trees.s[, 1:3], col=list(1, 2, 3), border=1:3, beanlines="median")
    

    Screen Shot 2019-01-09 at 3.30.25 PM.png

    Figure \(\PageIndex{8}\) Bean plot with overall line and median lines (default lines are means).