¿Cómo descartar los factores no utilizados en el diagrama de caja de ggplot facetado?

A continuación hay un código de ejemplo que utilizo para hacer algunas gráficas de caja:

stest <- read.table(text="    site  year    conc
    south   2001    5.3
    south   2001    4.67
    south   2001    4.98
    south   2002    5.76
    south   2002    5.93
    north   2001    4.64
    north   2001    6.32
    north   2003    11.5
    north   2003    6.3
    north   2004    9.6
    north   2004    56.11
    north   2004    63.55
    north   2004    61.35
    north   2005    67.11
    north   2006    39.17
    north   2006    43.51
    north   2006    76.21
    north   2006    158.89
    north   2006    122.27
", header=TRUE)

require(ggplot2)
ggplot(stest, aes(x=year, y=conc)) +
  geom_boxplot(horizontal=TRUE) +
  facet_wrap(~site, ncol=1) +
  coord_flip() +
  scale_y_log10()

Lo que resulta en esto:

Intenté todo lo que pude pensar pero no puedo hacer una trama donde elfaceta sur solo contiene años donde se muestran los datos (2001 y 2002). ¿Es lo que estoy tratando de hacer posible?

Aquí hay unenlazar (MUERTO) a la captura de pantalla que muestra lo que quiero lograr:

Respuestas a la pregunta(2)

Su respuesta a la pregunta