Problema
Cuando intentamos representar una serie continua, en nuestro ejemplo WS (Winter Solstice), ggplot2 une el último dato de WS (invierno) en marzo con el primero en diciembre, obteniendo el siguiente resultado indeseado.
library(ggplot2)
getSeason <- function(DATES) {
#found here https://stackoverflow.com/questions/9500114/find-which-season-a-particular-date-belongs-to
WS <- as.Date("2012-12-15", format = "%Y-%m-%d") # Winter Solstice
SE <- as.Date("2012-3-15", format = "%Y-%m-%d") # Spring Equinox
SS <- as.Date("2012-6-15", format = "%Y-%m-%d") # Summer Solstice
FE <- as.Date("2012-9-15", format = "%Y-%m-%d") # Fall Equinox
# Convert dates from any year to 2012 dates
d <- as.Date(strftime(DATES, format="2012-%m-%d"))
ifelse (d >= WS | d < SE, "Winter",
ifelse (d >= SE & d < SS, "Spring",
ifelse (d >= SS & d < FE, "Summer", "Fall")))
}
zz <- sample(1:10000,365)/1000
dag <- seq(as.Date("2014-01-01"), as.Date("2014-12-31"), by = "day")
seas <- getSeason(dag)
test <- data.frame(zz,dag,seas)
ggplot(data=test, aes(x=dag,ymax=zz,ymin=0,fill=seas))+
geom_ribbon()
Solución
Dividimos los datos con dplyr en dos subconjuntos, igual o mayor que el 2014-12-15 y por debajo de dicha fecha. Así convertimos la serie WS en dos secciones discretas, evitando que una los dos extremos de WS.
library(dplyr)
ggplot() +
geom_ribbon(data = filter(test, dag >= "2014-12-15") ,
aes(x = dag, ymax = zz, ymin = 0, fill = seas)) +
geom_ribbon(data = filter(test, dag < "2014-12-15") ,
aes(x = dag, ymax = zz, ymin = 0, fill = seas))
Resultados
Referencias
Nube de datos