Goals

• Review some basic graphical formats and when they are useful

• Make graphs in R, such as histograms, bar charts, box plots, and scatter plots

• Critique basic graphs to improve readability and accurate communication

Get the data for all labs here

Start a script!

For this and every lab or project, begin by starting a new script, create a good header section and table of contents, save the script file with an informative name, and set your working directory. Aim to make the script useful as a future reference for doing things in R - this will come in handy for projects and assessments!

Learning the tools

{ggplot2} package

The function ggplot() allows us to graph most kinds of data relatively simply. Its syntax is slightly odd but very flexible. We’ll practice specific commands for several types of plots below.

To begin, remember to install ggplot2 if need to and then load ggplot2 with the library() function:

# Neat trick for checking whether a package is installed
# and installing it if it is not, loading it if it is...

if(!require("ggplot2")) install.packages("ggplot2")

# The "hard" way
# install.packages("ggplot2", dep = T) # run if needed...

library(ggplot2)

Graphing in {ggplot2} is slightly different than with the base R plot functions (like plot(), boxplot(), hist(), etc.) . To make a graph with ggplot(), you need to specify at least two elements in your command.

1. First part: Use the function ggplot() itself, to specify which data frame you want to use and also which variables are to be plotted.

2. Second part: Tell R what kind of graph to make using a geom() function.

The odd part is that these two parts are put together with a + sign. It’s simplest to see this with an example. We’ll draw a histogram with ggplot() in the next section.

Histograms

A histogram represents the frequency distribution of a numerical variable in a sample.

Let’s see how to make a basic histogram using the age data from the Titanic data set. Make sure you have loaded the data (using read.csv) into a data frame called titanicData.

# Example code to load the data
# NB you will probably have to supply your own file path
titanic <- read.csv("data/titanic.csv")

# Code to make a simple histogram of age:

ggplot(titanic, aes(x=age)) + geom_histogram()

# There are a couple of informative warnings, but we can ignore them
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.
## Warning: Removed 680 rows containing non-finite values (stat_bin).

ggplot()

Notice that there are two functions called here, put together in a single command with a plus sign. The first function is ggplot(), and it has two input arguments. Listed first is titanicData; this is the name of the data frame containing the variables that we want to graph. The second input to ggplot is an aes() function. In this case, the aes() function tells R that we want Age to be the x-variable (i.e. the variable that is displayed along the x-axis). (The “aes” stands for “aesthetics”,” but if you’re like me it won’t necessarily help you remember it any better)

geom_histogram()

The second function in this command is geom_histogram(). This is the part that tells R that the “geometry” of our plot should be a histogram.

This is not the most beautiful graph in the world, but it conveys the information. Later, we’ll see a couple of options that can make a ggplot graph look a little better.

Bar graphs

A bar graph plots the frequency distribution of a categorical variable.

In ggplot(), the syntax for a bar graph is very similar to that for a histogram. For example, here is a bar graph for the categorical variable Sex in the titanic data set. Aside from specifying a different variable for x, we use a different geom function here, geom_bar.

ggplot(titanic, aes(x=sex)) + geom_bar(stat="count")

Boxplots

A boxplot is a convenient way of showing the frequency distribution of a numerical variable in multiple groups. Here’s the code to draw a boxplot for age in the titanic data set, separately for each sex:

ggplot(titanic, aes(x=sex, y=age)) + geom_boxplot()

Notice that the y variable here is age, and x is the categorical variable sex that winds up on the x-axis. See the result below, and look at where the variables are. The other new feature here is the new geom function, geom_boxplot().

Here the thick bar in the middle of each boxplot is the median of that group. The upper and lower bounds of the box extend from the first to the third quartile. The vertical lines are called whiskers, and they cover most of the range of the data (except when data points are pretty far from the median (see text), when they are plotted as individual dots, as on the male boxplot).

Scatterplots

The last graphical style that we will cover here is the scatter plot, which shows the relationship between two numerical variables.

We will use a dataset showing the relationship between the ornamentation of father guppies and the sexual attractiveness of their sons, using the file guppy-attractiveness.csv.

# Example code to load the data
# NB you will probably have to supply your own file path
guppy <-read.csv("data/guppy-attractiveness.csv", header =T)

head(guppy)

  fatherOrnamentation sonAttractiveness
1                0.35             -0.32
2                0.03             -0.03
3                0.14              0.11
4                0.10              0.28
5                0.22              0.31
6                0.23              0.18`

To make a scatter plot of the variables fatherOrnamentation and sonAttractiveness with ggplot, you need to specify the x and y variables, and use geom_point():

ggplot(guppy,  
      aes(x = fatherOrnamentation, y = sonAttractiveness)) +
  geom_point()

Better looking graphics with options

(sexing up the sexual selection scatterplot)

The code we have tried here for {ggplot2} does not begin to scratch the surface of what ggplot and R are capable of. Not only are there far more choices about the kinds of plots available, but there are many, many options for customizing the look and feel of each graph. You can choose the font, the font size, the colors, the style of the axes labels, etc., and you can customize the legends and axes legends nearly as much as you want.

Let’s dig a little deeper into just a few options that you can add to any of the graphs to improve them. For example, you can change the text of the x-axis label or the y-axis label by using xlab or ylab. Let’s do that for the scatterplot, to make the labels a little nicer to read for humans.  

ggplot(guppy,   
    aes(x = fatherOrnamentation, y = sonAttractiveness)) +
    geom_point() +
    xlab("Father ornamentation") + 
    ylab("Son attractiveness")

The can see the labels included in quotes inside the xlab and ylab arguments now appear on your graph.

It can also be nice to remove the default gray background, to make what some feel is a cleaner graph. Try adding + theme_classic() to the end of one of your lines of code making a graph, to see whether you prefer the result to the default design. Also, you could add a title by adding the labs function and setting the main argument to a character string + labs(title = "your title").

ggplot(guppy,   
    aes(x = fatherOrnamentation, y = sonAttractiveness)) +
    geom_point() +
    xlab("Father ornamentation") + 
    ylab("Son attractiveness") +
    theme_classic() +
    labs(title = "Ooooo, now that is classic")
    

Getting help

The help pages in R (using help() with the name of the function you would like help with) are the main source of help, but the amount of detail might be off-putting for beginners. For example, to explore the options for ggplot(), enter the following into the R Console.

help(ggplot)

This will cause the contents of the manual page for this function to appear in the Help window in RStudio. These manual pages are often frustratingly technical. What many of us do instead is simply google the name of the function-there are a great number of resources online about R.

There are also many introductory books available. A good one is Dalgaard (2008) Introductory Statistics with R, 2nd ed.

Challenge questions

Make a script in RStudio that collects all your R code required to answer the following questions. Include answers to the qualitative questions using comments.

1. For each of the following pairs of graphs, identify features that communicate better on one version than the other.

Survivorship as a function of sex for passengers of the RMS Titanic

Ear length in male humans as a function of age

2. Let’s use the data from “countries.csv” to practice making some graphs.

Read the data from the file “countries.csv” in the Data folder.

Use ggplot() (Make sure that you have run library(ggplot2)).

Make a histogram to show the frequency distribution of values for measles_immunization_oneyearolds, a numerical variable. (This variable gives the percentage of 1-year-olds that have been vaccinated against measles.) Describe the pattern that you see.

Make a bar graph to show the numbers of countries in each of the continents. (The categorical variable continent indicates the continent to which countries belong.)

Draw a scatterplot that shows the relationship between the two numerical variables life_expectancy_at_birth_male and life_expectancy_at_birth_female.

3. The ecological footprint is a widely-used measure, developed at UBC, of the impact a person has on the planet. It measures the area of land (in hectares) required to generate the food, shelter, and other resources used by a typical person and required to dispose of that person’s wastes. Larger values of the ecological footprint indicate that the typical person from that country uses more resources.

The countries data set has two variables for many countries showing the ecological footprint of an average person in each country. ecological_footprint_2000 and ecological_footprint_2012 show the ecological footprints for the years 2000 and 2012, respectively.

Plot the relationship between the ecological footprint of 2000 and of 2012.

Describe the relationship between the footprints for the two years. Does the value of ecological footprint of 2000 seem to predict anything about its value in 2012?

From this graph, does the ecological footprint tend to go up or down in the years between 2000 and 2012? Did the countries with high or low ecological footprint change the most over this time? (Hint: you can add a one-to-one line to your graph by adding + geom_abline(intercept = 0, slope = 1) to your ggplot command. This will make it easier to see when your points are above or below the “line of equivalence”).

4. Use the countries data again. Plot the relationship between continent and female life expectancy at birth. Describe the patterns that you see.

5. Muchhala (2006) measured the length of the tongues of eleven different species of South American bats, as well as the length of their palates (to get an indication of the size of their mouths). All of these bats use their tongues to feed on nectar from flowers. Data from the article are given in the file “BatTongues.csv”. In this file, both Tongue Length and Palette Length are given in millimeters.

Import the data and inspect it using summary(). You can call the data set whatever you like, but in one of the later steps we’ll assume it is called bat_tongues. Each value for tongue length and palate length is a species mean, calculated from a sample of individuals per species.

Draw a scatter plot to show the association between palate length and tongue length, with tongue length as the response variable. Describe the association: is it positive or negative? Is it strong or weak?

All of the data points that went into this graph have been double checked and verified. With that in mind, what conclusion can you draw from the outlier (the point with an unusual value) on the scatter plot?

Let’s figure out which species is the outlier. To do this, we’ll use a useful function called filter() from the package {dplyr}. Use library to load dplyr to your R session.

The function filter() gives us the row (or rows) of a data frame that has a certain property. Looking at the graph, we can tell that the point we are interested in has a very long tongue length, at least over 80 mm long! The following command will pull out the rows of the data frame bat_tongues that have tongue_length greater than 80 mm:

filter(bat_tongues,tongue_length>80)

The unusual species is Anoura fistulata See a photo here. This species has an outrageously long tongue, which it uses to collect nectar from a particular flower (can you guess what feature of the flower has led to the evolution of such a long tongue?). See the article by Muchhala (2006) to learn more about the biology of this strange bat.

6. Pick one of the plots you made using R today. What could be improved about this graph to make it a more effective presentation of the data?