R - ordering in boxplot

Question

I am trying to produce a series of box plots in R that is grouped by 2 factors. I've managed to make the plot, but I cannot get the boxes to order in the correct direction.

My data farm I am using looks like this:

Nitrogen    Species    Treatment
2           G          L
3           R          M
4           G          H
4           B          L
2           B          M
1           G          H

I tried:

boxplot(mydata$Nitrogen~mydata$Species*mydata$Treatment)

this ordered the boxes alphabetically (first three were the "High" treatments, then within those three they were ordered by species name alphabetically).

I want the box plot ordered Low>Medium>High then within each of those groups G>R>B for the species.

So i tried using a factor in the formula:

f = ordered(interaction(mydata$Treatment, mydata$Species), 
            levels = c("L.G","L.R","L.B","M.G","M.R","M.B","H.G","H.R","H.B")

then:

boxplot(mydata$Nitrogen~f)

however the boxes are still shoeing up in the same order. The labels are now different, but the boxes have not moved.

I have pulled out each set of data and plotted them all together individually:

lg = mydata[mydata$Treatment="L" & mydata$Species="G", "Nitrogen"]
mg = mydata[mydata$Treatment="M" & mydata$Species="G", "Nitrogen"]
hg = mydata[mydata$Treatment="H" & mydata$Species="G", "Nitrogen"]
etc ..

boxplot(lg, lr, lb, mg, mr, mb, hg, hr, hb)

This gives what i want, but I would prefer to do this in a more elegant way, so I don't have to pull each one out individually for larger data sets.

---

Loadable data:

mydata <-
structure(list(Nitrogen = c(2L, 3L, 4L, 4L, 2L, 1L), Species = structure(c(2L, 
3L, 2L, 1L, 1L, 2L), .Label = c("B", "G", "R"), class = "factor"), 
    Treatment = structure(c(2L, 3L, 1L, 2L, 3L, 1L), .Label = c("H", 
    "L", "M"), class = "factor")), .Names = c("Nitrogen", "Species", 
"Treatment"), class = "data.frame", row.names = c(NA, -6L))

Answer 1

The following commands will create the ordering you need by rebuilding the Treatment and Species factors, with explicit manual ordering of the levels:

mydata$Treatment = factor(mydata$Treatment,c("L","M","H"))

mydata$Species = factor(mydata$Species,c("G","R","B"))

---

edit 1 : oops I had set it to HML instead of LMH. fixing.

edit 2 : what factor(X,Y) does:

If you run factor(X,Y) on an existing factor, it uses the ordering of the values in Y to enumerate the values present in the factor X. Here's some examples with your data.

> mydata$Treatment
[1] L M H L M H
Levels: H L M
> as.integer(mydata$Treatment)
[1] 2 3 1 2 3 1
> factor(mydata$Treatment,c("L","M","H"))
[1] L M H L M H                               <-- not changed
Levels: L M H                                 <-- changed
> as.integer(factor(mydata$Treatment,c("L","M","H")))
[1] 1 2 3 1 2 3                               <-- changed

It does NOT change what the factor looks like at first glance, but it does change how the data is stored.

What's important here is that many plot functions will plot the lowest enumeration leftmost, followed by the next, etc.

If you create factors simply using factor(X) then usually the enumeration is based upon the alphabetical order of the factor levels, (eg "H","L","M"). If your labels have a conventional ordering different from alphabetical (ie "H","M","L"), this can make your graphs seems strange.

At first glance, it may seem like the problem is due to the ordering of data in the data frame - ie if only we could place all "H" at the top and "L" at the bottom, then it would work. It doesn't. But if you want your labels to appear in the same order as the first occurrence in the data, you can use this form:

 mydata$Treatment = factor(mydata$Treatment, unique(mydata$Treatment))

Answer 2

This earlier StackOverflow question shows how to reorder a boxplot based on a numerical value; what you need here is probably just a switch from factor to the related type ordered . But it is hard say as we do not have your data and you didn't provide a reproducible example.

Edit Using the dataset you posted in variable md and relying on the solution I pointed to earlier, we get

R> md$Species <- ordered(md$Species, levels=c("G", "R", "B"))
R> md$Treatment <- ordered(md$Treatment, levels=c("L", "M", "H"))
R> with(md, boxplot(Nitrogen ~ Species * Treatment))

which creates the chart you were looking to create.

This is also equivalent to the other solution presented here.