The box-whisker plot (or a boxplot) is a quick and easy way to visualize complex data where you have multiple samples.

A box plot is a good way to get an overall picture of the data set in a compact manner.

## The boxplot() function

You can use the `boxplot()`

function to create box-whisker plots.

It has many options and arguments to control many things, such as the making it horizontal, adding labels, titles and colors.

### Syntax

The syntax for the `boxplot()`

function is:

boxplot(x,names,xlab,ylab,border,col,notch,horizontal,add,…)

### Parameters

Parameter | Description |

x | A vector of values from which the boxplots are to be produced |

names | Group labels to be printed under each boxplot |

xlab | The label for the x axis |

ylab | The label for the y axis |

border | A vector of colors for the outlines of the boxplots |

col | The foreground color of symbols as well as lines |

notch | if TRUE, a notch is drawn in each side of the boxes |

horizontal | Set it to TRUE to draw the box-plot horizontally |

add | Set it to TRUE to add boxplot to current plot |

… | other graphical parameters |

## Create a Box-Whisker Plot

To get started with plot, you need a set of data to work with. Let’s consider the built-in ToothGrowth data set as an example data set.

Here are the first six observations of the data set.

```
# First six observations of the 'ToothGrowth' data set
head(ToothGrowth)
len supp dose
1 4.2 VC 0.5
2 11.5 VC 0.5
3 7.3 VC 0.5
4 5.8 VC 0.5
5 6.4 VC 0.5
6 10.0 VC 0.5
```

**ToothGrowth data set**

ToothGrowth data set contains observations on effect of vitamin C on tooth growth in 60 guinea pigs, where each animal received one of three dose levels of vitamin C (0.5, 1, and 2 mg/day) by one of two delivery methods, orange juice (coded as OJ) or ascorbic acid (coded as VC).

To create a box plot just specify any variable of the data set in `boxplot()`

function.

```
boxplot(ToothGrowth$len)
```

## Horizontal Box Plot

You can also draw the box-plot horizontally by setting the horizontal argument to TRUE.

```
boxplot(ToothGrowth$len,
horizontal = TRUE)
```

## Notched Box Plot

The notched box plot allows you to assess whether the medians are different. If the notches do not overlap, there is **strong evidence** (95% confidence) their medians differ.

You add notches to a box plot by setting the notch argument to TRUE.

```
# Add notches to a box plot
boxplot(ToothGrowth$len,
notch = TRUE)
```

## Side-by-Side Box Plots

Often your data set contains a numeric variable (quantitative variable) and a factor (categorical variable). It can be quite tedious to find whether the numeric variable changes according to the level of the factor.

Information of that nature can be gained by plotting box plots **side by side**.

In R, you can do this by using the boxplot() function with a formula:

boxplot(x ~ f)

Here, x is the numeric variable and f is the factor.

```
# Creating one box plot for each factor level (dose)
boxplot(len ~ dose, data = ToothGrowth)
```

## Grouped Box Plot

A grouped box plot is used when you have a numerical variable, several groups and subgroups.

You can create a grouped box plot by putting interaction of two categorical variables on x-axis and a numeric variable on y-axis.

The interaction of two variables is indicated by separating their names with an asterisk `*`

```
# Box plot of length based on interaction of two variables (supplement and dose)
boxplot(len ~ supp*dose, data = ToothGrowth)
```

## Change Group Names

To change names for group of boxes, use names argument.

```
boxplot(len ~ dose, data = ToothGrowth,
names=c("0.5 mg","1 mg","2 mg"))
```

## Change Colors

Use col argument to change the fill colors used for the boxes.

```
boxplot(len ~ dose, data = ToothGrowth,
col = "dodgerblue1")
```

You can change the colors of individual boxes by passing a vector of colors to the col argument.

```
boxplot(len ~ dose, data = ToothGrowth,
col = c("orange1", "dodgerblue1", "olivedrab2"))
```

By using the border argument, you can even change the color used for the border of the boxes.

```
boxplot(len ~ dose, data = ToothGrowth,
col="lightblue1",
border="dodgerblue3")
```

## Adding Titles and Axis Labels

You can add your own title and axis labels easily by specifying following arguments.

Argument | Description |

main | Main plot title |

xlab | x‐axis label |

ylab | y‐axis label |

```
boxplot(len ~ dose, data = ToothGrowth,
main="Tooth Growth in Guinea Pigs",
xlab="Vitamin C dose (mg/day)",
ylab="Length of odontoblasts")
```

## Add Means to a Box Plot

The horizontal line in the middle of a box plot is the median, not the mean.

The median alone will not help you understand if the data is normally distributed. So, you need to add mean markers on your box plot.

```
boxplot(len ~ dose, data=ToothGrowth,
col="dodgerblue1")
meanval <- by(ToothGrowth$len, ToothGrowth$dose, mean)
points(meanval, col="white", pch=18)
```