Graphics basics in R - Plot()

By: Emiley J Printer Friendly Format    

One of the greatest strengths of R, relative to other programming languages, is the ease with which we can create publication-quality graphics. In this tutorial, you'll learn about base graphics in R.

We do not cover the more advanced portions of graphics in R in this tutorial. These include lattice, ggplot2 and ggvis.

Load the included data frame cars with data(cars).

> data(cars)

To fix ideas, we will work with simple data frames. Our main goal is to introduce various plotting functions and their arguments. All the output would look more interesting with larger, more complex data sets.
Pull up the help page for cars.

> ?cars

As you can see in the help page, the cars data set has only two variables: speed and stopping distance. Note that the data is from the 1920s.
Run head() on the cars data.

> head(cars)  speed dist1     4    22     4   103     7    44     7   225     8   166     9   10

Before plotting, it is always a good idea to get a sense of the data. Key R commands for doing so include, dim(), names(), head(), tail() and summary().

Run the plot() command on the cars data frame.

> plot(cars)

As always, R tries very hard to give you something sensible given the information that you have provided to it. First, R notes that the data frame you have given it has just two columns, so it assumes that you want to plot one column versus the other.

Second, since we do not provide labels for either axis, R uses the names of the columns. Third, it creates axis tick marks at nice round numbers and labels them accordingly. Fourth, it uses the other defaults supplied in plot().
We will now spend some time exploring plot, but many of the topics covered here will apply to most other R graphics functions. Note that 'plot' is short for scatterplot.

Look up the help page for plot().

> ?plot

The help page for plot() highlights the different arguments that the function can take. The two most important are x and y, the variables that will be plotted. For the next set of questions, include the argument names in your answers. That is, do not type plot(cars$speed, cars$dist), although that will work. Instead, use plot(x = cars$speed, y = cars$dist).

Use plot() command to show speed on the x-axis and dist on the y-axis from the cars data frame. Use the form of the plot command in which vectors are explicitly passed in as arguments for x and y.

> plot(x=cars$speed, y=cars$dist)

Note that this produces a slightly different answer than plot(cars). In this case, R is not sure what you want to use as the labels on the axes, so it just uses the arguments which you pass in, data frame name and dollar signs included.

Note that there are other ways to call the plot command, i.e., using the "formula" interface. For example, we get a similar plot to the above with plot(dist ~ speed, cars). However, we will wait till later in the lesson before using the formula interface.

Use plot() command to show dist on the x-axis and speed on the y-axis from the cars data frame. This is the opposite of what we did above.

> plot(x=cars$dist, y=cars$speed)

It probably makes more sense for speed to go on the x-axis since stopping distance is a function of speed more than the other way around. So, for the rest of the questions in this portion of the lesson, always assign the arguments accordingly.

Recreate the plot with the label of the x-axis set to "Speed".

> ?plot> plot(x=cars$speed, y=cars$dist, xlab="Speed")

Recreate the plot with the label of the y-axis set to "Stopping Distance".

> plot(x=cars$speed, y=cars$dist,  ylab="Stopping Distance")

Recreate the plot with "Speed" and "Stopping Distance" as axis labels.

> plot(x=cars$speed, y=cars$dist, xlab="Speed", ylab="Stopping Distance")

The reason that plots(cars) worked at the beginning of the lesson was that R was smart enough to know that the first element (i.e., the first column) in cars should be assigned to the x argument and the second element to the y argument. To save on typing, the next set of answers will all be of the form, plot(cars, ...) with various arguments added.

For each question, we will only want one additional argument at a time. Of course, you can pass in more than one argument when doing a real project.
Plot cars with a main title of "My Plot". Note that the argument for the main title is "main" not "title".

> plot(cars, main="My Plot")

Plot cars with a sub title of "My Plot Subtitle".

> ?plot> plot(cars, sub="My Plot Subtitle")

The plot help page (?plot) only covers a small number of the many arguments that can be passed in to plot() and to other graphical functions. To begin to explore the many other options, look at ?par. Let's look at some of the more commonly used ones. Continue using plot(cars, ...) as the base answer to these questions.

Plot cars so that the plotted points are colored red. (Use col = 2 to achieve this effect.)

> plot(cars, col=2)

Plot cars while limiting the x-axis to 10 through 15.  (Use xlim = c(10, 15) to achieve this effect.)

> plot(cars, xlim=c(10,15))

You can also change the shape of the symbols in the plot. The help page for points (?points) provides thedetails. Plot cars using triangles.  (Use pch = 2 to achieve this effect.)

> plot(cars, pch=2)

Arguments like "col" and "pch" may not seem very intuitive. And that is because they aren't! So, many/most people use more modern packages, like ggplot2, for creating their graphics in R.

It is, however, useful to have an introduction to base graphics because many of the idioms in lattice andggplot2 are modeled on them.

Let's now look at some other functions in base graphics that may be useful, starting with boxplots. Load the mtcars data frame.

> data(mtcars)

Anytime that you load up a new data frame, you should explore it before using it. In the middle of a swirl lesson, just type play(). This temporarily suspends the lesson (without losing the work you have already done) and allows you to issue commands like dim(mtcars) and head(mtcars). Once you are done examining the data, just type nxt() and the lesson will pick up where it left off.

Look up the help page for boxplot().

> ?boxplot

Instead of adding data columns directly as input arguments, as we did with plot(), it is often handy to pass in the entire data frame. This is what the "data" argument in boxplot() allows.

boxplot(), like many R functions, also takes a "formula" argument, generally an expression with a tilde ("~") which indicates the relationship between the input variables. This allows you to enter something like mpg ~ cyl to plot the relationship between cyl (number of cylinders) on the x-axis and mpg (miles per gallon) on the y-axis.

Use boxplot() with formula = mpg ~ cyl and data = mtcars to create a box plot.

> boxplot(formula = mpg ~ cyl, data = mtcars)

The plot shows that mpg is much lower for cars with more cylinders. Note that we can use the same set of arguments that we explored with plot() above to add axis labels, titles and so on.

When looking at a single variable, histograms are a useful tool. hist() is the associated R function. Like plot(), hist() is best used by just passing in a single vector.

Use hist() with the vector mtcars$mpg to create a histogram.

> hist(mtcars$mpg)

In this tutorial, you learned how to work with base graphics in R. The best place to go from here is to study the ggplot2 package. If you want to explore other elements of base graphics, then this web page (http://www.ling.upenn.edu/~joseff/rstudy/week4.html) provides a useful overview.





Most Viewed Articles (in Data Science )

Latest Articles (in Data Science)

Comment on this tutorial