Soy un novato en R, así que perdona mi ignorancia. Hice una gráfica de barras pseudo-apilada en la que dibujé 4 conjuntos de barras una encima de la otra utilizando geom_bar. Hay 4 categorías de estado de salud (vivo, muerto, infectado y muerto) para tres especies de robles (QUAG, QUKE, QUCH).

Mi código es el siguiente:

x <- as.data.frame(list(variable=c("QUAG", "QUKE", "QUCH"), alive = c(627,208,109),  infected = c(102,27,0), dead = c(133,112,12), sod.dead=c(49,8,0)))

x.plot = ggplot(x, aes(variable, alive)) + geom_bar(fill="gray85") + 
  geom_bar(aes(variable,dead), fill="gray65") +
  geom_bar(aes(variable, infected), fill="gray38") +
  geom_bar(aes(variable, sod.dead), fill="black")+
  opts(panel.background = theme_rect(fill='gray100'))
x.plot

Ahora quiero hacer una leyenda que muestre qué tono de gris se relaciona con el estado del árbol, es decir, "gray65" es "árboles muertos", etc. Lo he intentado durante la última hora y no puedo hacerlo funcionar.