r ggplot 丝带季节填充错误

Question

我无法正确地为时间序列中的组着色。好像顺序乱了示例代码如下。有什么提示吗？

library(ggplot2)

getSeason <- function(DATES) {
#found here 
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()

Answer 1

我运行这是为了测试你的数据，看起来你的 winter 在数据的两端，这导致 2014-03-14 Winter 连接到 2014-12-15 Winter其他季节。

ggplot(test, aes(dag)) + 
   geom_ribbon(aes(ymin = zz - 1, ymax = zz + 1), fill = "grey70") +
   geom_line(aes(y = zz, group = seas)) + facet_grid(. ~ seas)

你也可以在这里看到，

ggplot(test, aes(x= dag, y= zz)) +
geom_line(aes(colour= seas)) +
geom_ribbon(aes(ymax= zz, ymin= zz, fill= seas), alpha = 0.2)

一个可能的解决方案是像这样对数据进行子集化，

ggplot() + 
 geom_line(data=subset(test, dag < as.Date("2014-12-15") ), aes(dag,zz,colour=seas))  + 
 geom_line(data=subset(test, dag > as.Date("2014-12-15") ), aes(dag,zz,colour=seas))

Answer 2

这是一个通用的解决方案，需要 lubridate 包。它也不需要难以阅读的 ifelse() 语法，而是使用 cut() 将日期向量拆分为与每个季节相对应的因子级别。

函数 return 给定日期的季节特征向量

getSeason <- function(DATES) {
  require(lubridate)

  # Get day of year of solstices
  solstices <- as.Date(c('3-15','6-15','9-15','12-15'),format = '%m-%d')
  solstice_day <- yday(solstices)

  # get year and day of year of dates
  dates_year <- year(DATES)
  dates_day <- yday(DATES)

  # Split dates vector into seasons
  season_label <- as.character(cut(dates_day, 
                                   breaks = c(0,solstice_day,367), 
                                   labels = c('winter','spring','summer','fall','winter_next')))

  # Deal with the second winter
  dates_year[season_label == 'winter_next'] <- dates_year[season_label == 'winter_next'] + 1
  season_label[season_label == 'winter_next'] <- 'winter'

  # Return seasons pasted with year
  paste0(season_label, dates_year)

}

绘制色带的代码

library(ggplot2)
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() +
  scale_fill_manual(values = c('winter2014'='red','spring2014'='blue','summer2014'='green','fall2014'='yellow','winter2015'='red'))

结果

Answer 3

您可以通过将您的数据分成两个子集，在 WS 之上和之下，并绘制两个 geom_ribbon 来解决它。因此，您将连续的 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))

Answer 4

您可以在两个冬天之间的某处添加一个虚拟缺失值：

test[nrow(test), ] <- list(NA, as.Date("2014-6-1"), "Winter")

ggplot(data=test, aes(x=dag,ymax=zz,ymin=0,fill=seas))+
  geom_ribbon()

r ggplot 丝带季节填充错误

r ggplot ribbon season fill wrong

r

date

fill

ggplot2