¿Cómo ordenar correctamente los diagramas de caja de facetas por mediana?

Estoy usando el conjunto de datos de "diamantes" que viene con R. Al intentar clasificar el factor de "color" con respecto a su mediana de precios, no funcionará.

Esto es lo que conseguí:

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

Y me da eso (no está ordenado en absoluto):

¿Hay algo que estoy haciendo mal o me falta?

Respuestas a la pregunta(2)

Su respuesta a la pregunta