Выделение значений выше порога с помощью решетки в R xyplot

Такой график легко построить с помощью ggplot2 (в этом примере используется примерный набор данных mtcars):

library(ggplot2)
ggplot(mtcars, aes(x = wt, y = mpg, color = mpg > 30)) + geom_point() + 
    scale_color_manual(values = c('red', 'blue'), breaks = c(TRUE, FALSE))

Я использовал отметку даты и времени для форматирования столбца даты. Затем я использовал предложение @PaulHiemstra с ggplot2.

образец д.ф.

Date    sseye   max month   season
2009/10/6   Shaded Nest-boxes   29.9508 11  Non-breeding season
2011/06/7   Sun Exposed Nest-boxes  38.1586 11  Non-breeding season
2013/11/10  Shaded Nest-boxes   29.2015 11  Non-breeding season
2014/05/9   Sun Exposed Nest-boxes  39.4666 11  Non-breeding season

library(scales) # for units
library(ggplot2)
library(gridExtra)
library(plyr)
library(ffbase)


Nestemp2$Date <- as.Date(Nestemp2$Date, format="%Y/%m/%d")
plotnest=ggplot(Nestemp2, aes(x=Date, y=max, color = max > 35)) + 
  geom_point() + scale_color_manual(values = c('black', 'red')) +  

Затем я хотел добавить вертикальные и горизонтальные линии в определенные даты и значение

geom_hline(yintercept=35, colour="red", linetype="dashed") + 
  geom_hline(yintercept=0, colour="black")+
  geom_vline(xintercept=as.numeric(as.Date("2010-04-01")), vjust=0.5, linetype=2) + 
  annotate("text", x=as.Date("2010-04-10"), y=25, label=" Not Breeding", angle=90, size=3,  vjust=-1, hjust=1) + 
  geom_vline(xintercept=as.numeric(as.Date("2010-10-01"), vjust=0.5), linetype=2) + 
  annotate("text", x=as.Date("2010-09-30"), y=25, label=" Breeding", angle=90, size=3,  vjust=-1, hjust=1) + 

Мой столбец df состоит из нескольких факторов для местоположения гнезда (солнце, тень и эталонное местоположение), поэтому я использовал обтекание граней, чтобы сделать панельный график, состоящий из фактора ~ssey (три уровня (солнце, тень и эталон). Я не могу опубликовать картинка, но выглядит мило.

ylab("Temperature (c)") +
  facet_wrap( ~ sseye, ncol=1) +                  #this is what puts them in a grid
  scale_y_continuous(limits=c(0,60)) +    

  theme_bw() + theme(
  strip.text.x = element_text(size=16),           #sets the size of the title to each grid section as 10
  strip.background = element_blank(),             #removes the grey background in the title of each grid section
  axis.text.x = theme_text(size=16),               #sets the x axis tick text size to 8
  axis.text.y = theme_text(size=16),               #sets the y axis tick text size to 8
  axis.title.x=theme_blank(),                     #removes the x axis title
  axis.title.y=theme_text(size=16, vjust = .3),               #sets the y axis title to size 10
  axis.line=theme_segment(colour="black"),        #sets the axis lines
  panel.grid.minor = theme_blank(),               #removes minor grid lines
  panel.grid.major = theme_blank(),               #removes major grid lines
  panel.border=theme_blank(),                     #removes the border around the graph
  panel.background=theme_blank(),                 #removes the grey background of the plot
  legend.justification=c(10,10), 
  legend.position=c(10,10), # Position legend in top right
  legend.title=element_blank(),
  legend.text = theme_blank(),
  legend.title = theme_blank(),                   #removes the legend title
  legend.key = theme_blank())   

plotnest