Como classificar corretamente os boxplots de faceta por mediana?

Estou usando o conjunto de dados 'diamantes' que acompanha o R. Ao tentar classificar o fator 'cor' em relação à sua mediana de preços, não funcionará.

Isto é o que eu tenho:

ggplot(diamonds, aes(x = reorder(color, -price, FUN=median), y = price)) + 
  geom_boxplot() + 
  facet_wrap(~cut) + 
  ylim(0, 5500)

E isso me dá o seguinte (não classificado):

Existe algo que estou fazendo errado ou faltando?

questionAnswers(2)

yourAnswerToTheQuestion